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Abstract 


Symmetry plays a fundamental role in physics. The quasi-degeneracy between 
single-particle orbitals (n,/, j = 1+ 1/2) and (n — 1,/+ 2, 7 = L+ 3/2) indicates 
a hidden symmetry in atomic nuclei, the so-called pseudospin symmetry (PSS). 
Since the introduction of the concept of PSS in atomic nuclei, there have been 
comprehensive efforts to understand its origin. Both splittings of spin doublets 
and pseudospin doublets play critical roles in the evolution of magic numbers 
in exotic nuclei discovered by modern spectroscopic studies with radioactive ion 
beam facilities. Since the PSS was recognized as a relativistic symmetry in 1990s, 
many special features, including the spin symmetry (SS) for anti-nucleon, and 
many new concepts have been introduced. In the present Review, we focus on the 


recent progress on the PSS and SS in various systems and potentials, including 
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extensions of the PSS study from stable to exotic nuclei, from non-confining to 
confining potentials, from local to non-local potentials, from central to tensor 
potentials, from bound to resonant states, from nucleon to anti-nucleon spectra, 
from nucleon to hyperon spectra, and from spherical to deformed nuclei. Open 
issues in this field are also discussed in detail, including the perturbative nature, 
the supersymmetric representation with similarity renormalization group, and the 
puzzle of intruder states. 
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1. Introduction 


The establishment of nuclear shell model is one of the most important 
milestones in nuclear physics. Similar to that of electrons orbiting in an atom, 
protons and neutrons in a nucleus form shell structures. The corresponding so- 
called magic numbers are found to be 2, 8, 20, 28, 50, and 82 for both protons 
and neutrons as well as 126 for neutrons in stable nuclei [1, 2]. The abundance 
of a nucleus with the magic numbers of proton and/or neutron is normally more 
than its neighboring nuclei. The magic numbers manifest themselves as a sudden 
jump in the plots of the two-nucleon separation energies [3], the a-decay half- 
lives, neutron-capture cross sections, and also the isotope shifts as functions of 
nucleon number [4]. They also appear as peaks in the abundance pattern in the 
solar systems in astrophysics. 

In order to understand these magic numbers, starting from some simple models 
such as the square well or harmonic oscillator (HO) potential with analytical 
solutions, nuclear physicists tried to solve the corresponding Schrödinger equa- 
tions. In 1949, independently, Haxel, Jensen, and Suess [1] and Mayer [2] 
introduced the spin-orbit (SO) potential which largely splits the states with high 
orbital angular momentum. In combination with the usual mean-field harmonic 
oscillator, square well, or Woods-Saxon potentials, the strong spin-orbit potential, 
although added by hand, excellently reproduces all traditional magic numbers 
in nuclear physics. Apart from the magic numbers, it also provides wonderful 


descriptions for nuclear ground-state and some low-lying excited-state properties. 
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Therefore, the substantial spin symmetry (SS) breaking between the spin doublets 
(n, 1, j = 1+ 1/2) is one of the most important concepts in nuclear structure. 

The success of the nuclear shell model or spin-orbit potential is unprecedented. 
For light nuclei (A < 25) the rotational features of nuclear spectra can be 
understood in a many-particle spherical shell-model framework in terms of 
the SU(3) coupling scheme of Elliott and Harvey [5-7]. By introducing the 
deformation-dependent oscillator length, Nilsson et al. [8, 9] extended the shell 
model to the deformed cases and built the foundation for describing not only 
the deformed nuclei but also nuclear rotation phenomena. Even nowadays, 
searching for new magic numbers and investigating shell structure evolution for 
unstable nuclei is still one of the key topics for the radioactive ion beam facilities 
worldwide [10]. 

After the successful reproduction of the magic numbers, the shell model 
with strong spin-orbit potential became the strongest candidate of the standard 
nuclear model and almost the entire nuclear physics community was exploring and 
enjoying the new physics brought in by this model. Here we mentioned “almost” 
because there were a few groups who were examining the nuclear shell model in 
a different way. 

These groups were not satisfied with simply reproducing the magic numbers 
or the splittings between the spin doublets. By examining the single-particle 
spectra, Hecht and Adler [11] and Arima, Harvey, and Shimizu [12] found 
the near degeneracy between two single-particle states with quantum numbers 
(n,l, j = 1+ 1/2) and (n — 1,1 + 2,7 = 1+ 3/2). They introduced the so-called 
pseudospin symmetry (PSS) and defined the pseudospin doublets as (ñ = n—1,/ = 


l+ 1, j =1+ 1/2) to explain this near degeneracy. This is illustrated in Fig. 1. 
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Figure 1: (Color online) Schematic nuclear single-particle spectrum. On one hand, the strong spin- 
orbit splittings between the spin doublets (n, l, j = L+ 1/2) lead to the traditional magic numbers. 


On the other hand, pairs of single-particle states in braces, (n,/, j = 1+ 1/2) and(n-1,1+2,j = 


1+3/2), are quasi-degenerate. They are defined as the pseudospin doublets (ñ = n—1,/ = 1+1, j = 


/ + 1/2), and the pseudospin symmetry is introduced for such near degeneracy. 


The pseudospin symmetry remains an important physical concept in axially 
deformed [13—16] and even triaxially deformed [17, 18] nuclei. Based on this 
concept, a simple but useful pseudo-SU(3) model was proposed, and this model 
was generalized to be the (pseudo-)symplectic model [15, 16, 19-22], (also see 
Ref. [23] and references therein). The pseudospin symmetries have been also 
extensively used in the odd-mass nuclei in the context of the interacting Boson- 
Fermion model [24-27]. 

Although the concept of pseudospin symmetry was based on the observation 
of empirical single-particle spectra, it remained to be a purely theoretical subject 
for nearly 20 years before the discovery of the nuclear superdeformed identical 
rotation bands [28]. From then on, a lot of phenomena in nuclear structure have 
been successfully interpreted directly or implicitly by the pseudospin symmetry, 
including nuclear superdeformed configurations [29, 30], identical bands [31-33], 
quantized alignment [34], and pseudospin partner bands [35, 36]. The pseudospin 
symmetry may also manifest itself in magnetic moments and transitions [37-39] 
and y-vibrational states in nuclei [40] as well as in nucleon-nucleus and nucleon- 
nucleon scatterings [41—44]. In addition, the relevance of pseudospin symmetry 
in the structure of halo nuclei [45] and superheavy nuclei [46, 47] was pointed 
out. 

In the 21st century, it has been intensively shown that the traditional magic 
numbers can change in nuclei far away from the stability line [10]. It is found 
that both splittings of spin and pseudospin doublets play critical roles in the shell 
structure evolution. For example, the N = 28 shell closure disappears due to 
the quenching of the spin-orbit splitting for the vl f spin doublets [48-51], and 


the Z = 64 subshell closure is closely related to the restoration of pseudospin 


symmetry for the 72) and 21f pseudospin doublets [52-54]. Therefore, it will 
be quite interesting and challenging to understand shell closure and pseudospin 
symmetry on the same footing, in particular for superheavy and exotic nuclei near 
the limit of nucleus existence. 

Since the recognition of pseudospin symmetry in atomic nuclei, there have 
been comprehensive efforts to understand its origin. Apart from the formal 
relabelling of quantum numbers, various explicit transformations from the normal 
scheme to the pseudospin scheme have been proposed [55-58]. In 1982, Bohr, 
Hamamoto, and Mottelson [55] discussed the pseudospin symmetry in rotating 
nuclear potentials, and found that such a symmetry is very helpful to understand 
qualitatively the feature of quasi-particle motion in rotating potentials. Based 
on the single-particle Hamiltonian of the harmonic oscillator shell model, they 
tried to understand the origin of pseudospin symmetry in terms of the spin- 
orbit potential introduced by hand, and also the orbit-orbit term, which has been 
artificially added in the Nilsson model. It turns out that the origin of pseudospin 
symmetry is connected with a special ratio between the strengths of the spin-orbit 
and orbit-orbit interactions. This idea was followed by the groups at Louisiana 
State University, University of California, and National Autonomous University 
of Mexico (UNAM), and they tried to understand the spin-orbit and orbit-orbit 
potentials in order to explain the pseudospin symmetry [30, 56-58]. 

The relation between the pseudospin symmetry and the relativistic mean-field 
(RMF) theory [59-64] was first noted in Ref. [30], where the relativistic mean- 
field theory was used to explain approximately such a special ratio between the 
strengths of the spin-orbit and orbit-orbit interactions. 


In order to see the connection with the relativistic mean-field or the covariant 


density functional theory (CDFT), it will be quite illuminating to examine the 
Dirac equation governing the motion of nucleons. The corresponding single- 
particle wave functions are given in the form of the Dirac spinor, which has both 
the upper and lower components. For the spherical case, by looking into the Dirac 
spinor, it is interesting to note that the upper and lower components have the same 
total angular momentum j but not the orbital angular momentum /. Their orbital 
angular momenta / differ by one unit. In the normal labelling of the single-nucleon 
states, the / of the dominant upper component is used. Equivalently, one can also 
use the quantum number of the lower component to label the single-nucleon states. 
In 1997, Ginocchio [65] revealed that the pseudospin symmetry is essentially a 
relativistic symmetry of the Dirac Hamiltonian and the pseudo-orbital angular 
momentum / is nothing but the orbital angular momentum of the lower component 
of the Dirac spinor. He also showed that the pseudospin symmetry in nuclei is 
exactly conserved when the scalar potential S(r) and the vector potential V(r) 
have the same size but opposite sign, i.e., X(r) = S(r) + V(r) = 0. This discovery 
not only reveals the origin of pseudospin symmetry but also demonstrates an 
unexpected success of the relativistic mean-field theory. It should be noted the 
pseudospin symmetry is a special case of the Bell-Reugg symmetries [66], as 
pointed out in Section 2 of Ref. [67]. 

One can also go one step further to reduce the Dirac equation into the second- 
order differential equation for either the upper or lower component. Then for 
the upper and lower components there will be the spin-orbit and pseudospin- 
orbit (PSO) potentials, respectively governing the relevant partner splittings. If 
either the spin-orbit or pseudospin-orbit potential vanishes, it will lead to the 


corresponding spin or pseudospin symmetry. The derivative for the sum of the 
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scalar and vector potentials, i.e., dX(r)/dr, determines the pseudospin symmetry. 
The pseudospin symmetry is exact under the condition dX(r)/dr = 0 [68]. This 
condition means that the pseudospin symmetry becomes much better for exotic 
nuclei with a highly diffused potential [69]. Approximately, the pseudospin 
symmetry is connected with the competition between the pseudo-centrifugal 
barrier (PCB) and the pseudospin-orbit potential. However, in either limit, Ł(r) = 
0 or dX(r)/dr = Q, there are no longer bound states, thus the pseudospin symmetry 
is always broken in realistic nuclei. In this sense, the pseudospin symmetry is 
viewed as a dynamical symmetry [70-72] or it is of the non-perturbative nature 
[73-76]. 

Following discussions for spherical nuclei, the study of pseudospin symmetry 
within the relativistic framework was quickly extended to deformed ones [77, 78]. 
As the pseudospin symmetry is a relativistic symmetry, the wave functions of the 
pseudospin partners satisfy certain relations. These relations have been tested in 
both spherical and deformed nuclei [77, 79-82]. 

Although the doubt on the connection between the pseudospin symmetry 
and the condition X(r) = 0 or dX(r)/dr = O exists [83-85], following the 
pseudospin symmetry limit, a lot of discussions about the pseudospin symmetry 
in single-particle spectra have been made by exactly or approximately solving the 
Dirac equation with various potentials, for examples, the one-dimensional Woods- 
Saxon potential [86], the two-dimensional Smorodinsky-Winternitz potential 
[87], the spherical harmonic oscillator [88—96], anharmonic oscillator [97], 
Coulomb [76, 98-101], Deng-Fan [102], diatomic molecular [103, 104], Eckart 
[105, 106], Hellmann [107], Hulthén [108-111], Manning-Rosen [112-114], 
Mie-type [115-117], Morse [118-123], Poéschl-Teller [124-133], Rosen-Morse 
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[134-137], Tietz-Hua [138], Woods-Saxon [139-143], and Yukawa [144-147] 
potentials, as well as the deformed harmonic oscillator [148-153], anharmonic 
oscillator [154, 155], Hartmann [156, 157], Hylleraas [158], Kratzer [159], 
Makarov [160], Manning-Rosen [161], and ring-shaped [162-165] potentials. 
Self-consistently, the pseudospin symmetry in spherical [45, 73, 79, 166-177] 
and deformed [77, 78, 81, 82, 178-180] nuclei have been investigated within 
relativistic mean-field and relativistic Hartree-Fock (RHF) [181—184] theories. 
One of interesting topics is the tensor effects on the pseudospin symmetry or spin 
symmetry, which has been investigated in some of the above-mentioned works 
[45, 53, 89, 92, 173] and also in Refs. [185, 186]. 

For the many-body problems in quantum mechanics, the basis expansion is 
one of the standard methods, e.g., the solutions of the Schrödinger equation with a 
harmonic oscillator potential as a widely used basis. For the Dirac equation, there 
exist not only the positive-energy states in the Fermi sea but also the negative- 
energy states in the Dirac sea, where the negative-energy states correspond to 
the anti-particle states. When the solutions of the Dirac equation are used as a 
complete basis, e.g., the Dirac Woods-Saxon basis [187], the states with both 
positive and negative energies must be included [184, 187-194]. 

When Zhou, Meng, and Ring developed the relativistic mean-field theory 
in a Dirac Woods-Saxon basis [187], they examined carefully the negative- 
energy states in the Dirac sea and found that the pseudospin symmetry of those 
negative-energy states, or equivalently, the spin symmetry in the anti-nucleon 
spectra is very well conserved [195]. Furthermore, they have shown that the 
spin symmetry in the anti-nucleon spectra is much better developed than the 


pseudospin symmetry in normal nuclear single-particle spectra. It should be noted 
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that, by applying the charge conjugate transformation, the spin symmetry for anti- 
nucleon states have been formally conjectured in Ref. [196]. The spin symmetry 
in the anti-nucleon spectra was also tested by investigating relations between 
Dirac wave functions of spin doublets with the relativistic mean-field theory [197]. 
Later, this spin symmetry was studied with the relativistic Hartree-Fock theory 
and the contribution from the Fock term was analyzed [198]. It has been pointed 
out in Ref. [195] that an open problem related to the experimental study of the 
spin symmetry in the anti-nucleon spectra is the polarization effect caused by 
the annihilation of anti-nucleons in a normal nucleus. Detailed calculations of 
the anti-baryon annihilation rates in the nuclear environment showed that the 
in-medium annihilation rates are strongly suppressed by a significant reduction 
of the reaction Q values, leading to relatively long-lived anti-baryon-nucleus 
systems [199]. Recently, the spin symmetry in the anti-A spectra of hypernuclei 
was studied quantitatively [200-202], which may be free from the problem 
of annihilation. This kind of study would be of great interests for possible 
experimental tests. 

In recent years, there has been an increasing interest in the exploration of 
continuum and resonant states, especially in the studies of exotic nuclei with 
unusual N/Z ratios. In exotic nuclei, the neutron (or proton) Fermi surface is close 
to the particle continuum; thus, the contribution of the continuum is important 
[190-194, 203-218]. Many methods or models developed for the studies of 
resonances [219] have been adopted to locate the position and to calculate the 
width of a nuclear resonant state, e.g., the analytical continuation in coupling 
constant (ACCC) method [220-225], the real stabilization method [226—231], the 
complex scaling method (CSM) [232-235], the coupled channels method [236— 
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238], and some others [239, 240]. 

The study of symmetries in resonant states is certainly interesting, e.g., the 
pseudospin symmetry [241-245] and spin symmetry [246] in single-particle 
resonant states. Recently, Lu, Zhao, and Zhou [247] gave a rigorous verification 
of the pseudospin symmetry in single-particle resonant states. They have shown 
that the pseudospin symmetry in single-particle resonant states in nuclei is exactly 
conserved under the same condition for the pseudospin symmetry in bound states, 
i.e., &(r) = 0 or dX(r)/dr = 0 [247-249]. The exact conservation and breaking 
mechanism of the pseudospin symmetry in single-particle resonances for spherical 
square-well potentials have been investigated, in which the pseudospin symmetry- 
breaking part can be separated from other parts in the Jost functions. By 
examining zeros of Jost functions corresponding to the lower components of radial 
Dirac wave functions, general properties of pseudospin symmetry splittings of the 
energies and widths are examined. As noted in Ref. [247], it is straightforward to 
extend the study of the pseudospin symmetry in resonant states in the Fermi sea to 
that in the negative-energy states in the Dirac sea or spin symmetry in anti-particle 
continuum spectra. 

Works are also in progress for understanding the origin of pseudospin 
symmetry and its breaking mechanism in a perturbative way. On one hand, the 
perturbation theory was used in Refs. [250, 251] to investigate the symmetries of 
the Dirac Hamiltonian and their breaking in realistic nuclei. This provides a clear 
way for investigating the perturbative nature of pseudospin symmetry. An illu- 
minating example is that the energy splittings of the pseudospin doublets can be 
regarded as a result of perturbation of the Hamiltonian with a relativistic harmonic 


oscillator (RHO) potential, where the pseudospin doublets are degenerate [250]. 
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On the other hand, supersymmetric (SUSY) quantum mechanics [252, 253] 
was used to investigate the symmetries of the Dirac Hamiltonian [254-256], 
also see Refs. [143, 257-261]. In particular, by employing both exact and 
broken supersymmetries, the phenomenon that all states with 7 > O have 
their own pseudospin partners except for the so-called intruder states can be 
interpreted naturally within a unified scheme. A pseudospin symmetry-breaking 
potential without a singularity can also be obtained with the supersymmetric 
technique [255], in contrast to the singularities appearing in the reduction of the 
Dirac equation to a Schrédinger-like equation for the lower component of the 
Dirac spinor. However, by reducing the Dirac equation to a Schrédinger-like 
equation for the upper component, the corresponding effective Hamiltonian is 
not Hermitian, since the upper component wave functions alone, as the solutions 
of the Schrédinger-like equation, are not orthogonal to each other. In order to 
fulfill the orthonormality, an additional differential relation between the lower and 
upper components must be taken into account. By doing so, effectively, the upper 
components alone are orthogonal with respect to a different metric [255]. Such 
fact that the corresponding effective Hamiltonian is not Hermitian prevents us 
from being able to perform quantitative perturbation calculations. 

Recent works by Guo and coauthors [262-264] bridged the perturbation 
calculations and the supersymmetric descriptions by using the similarity renor- 
malization group (SRG) [265-269] to transform the Dirac Hamiltonian into a 
diagonal form. The effective Hamiltonian expanded in a 1/M series in the 
Schrédinger-like equation is Hermitian. This makes the perturbation calculations 
possible. Therefore, one can understand the origin of pseudospin symmetry and its 


breaking by combining supersymmetric quantum mechanics, perturbation theory, 
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and the similarity renormalization group [270, 271]. 

Another open issue in the study of pseudospin symmetry is the special status 
of nodeless intruder states which do not have their own pseudospin partners. 
The nodal structure of radial Dirac wave functions of pseudospin doublets was 
studied in Ref. [272], which was helpful particularly for understanding the 
reason why between a pair of pseudospin doublets the pseudospin-down state 
with j. = Ī— 1/2 has one more radial node than the pseudospin-up state with 
j> = Ī + 1/2. However, in this case, there exist no bound states in the Fermi sea 
at the pseudospin symmetry limit. In contrast, as pointed out in Ref. [88], there 
can exist bound states in the Fermi sea at the pseudospin symmetry limit if the 
potential A(r) = V(r) — S(r) is confining. Quite recently, the nodal structure of 
radial Dirac wave functions for this case was demonstrated in an analytical way 
by Alberto, de Castro, and Malheiro [273]. It is interesting to note that in such 
a case all states with / > 0 have their own pseudospin partners, but instead some 
states with / > 0 lose their own spin partners. 

In this Review, we will mainly focus on the progress in the studies of 
pseudospin symmetry and spin symmetry hidden in atomic nuclei and the related 
topics in the past decade. Section 3 will be devoted to highlighting the progress 
from several different aspects, and Section 4 will be devoted to discussing the 
selected open issues. Some of the topics covered in a former review [67] will not 
be repeated here, such as the Bell-Reugg symmetries, the pseudospin symmetry 
in the magnetic dipole, electric quadrupole, and Gamow-Teller transitions, the 
pseudospin symmetry in the nucleon-nucleon and nucleon-nucleus scattering, as 
well as the spin symmetry in hadrons. 


The paper will be organized as follows. The typical Dirac equation widely 
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used in nuclear physics will be presented together with its Schrédinger-like 
equations in Section 2. In the same Section, different analytical solutions for 
Dirac equation at the pseudospin symmetry limit and the pseudospin symmetry 
breaking in realistic nuclei will be reviewed briefly. Recent progress on the 
pseudospin symmetry, ranging from stable to exotic nuclei, non-confining to 
confining potentials, local to non-local potentials, central to tensor potentials, 
bound to resonant states, nucleon to anti-nucleon spectra, nucleon to hyperon 
spectra, and spherical to deformed nuclei, will be presented in Section 3. 
Section 4 will be devoted to discussing the open issues in this field, including the 
perturbative nature, puzzle of intruder states, and supersymmetric representation 
for pseudospin symmetry. Finally, summary and perspectives will be given in 


Section 5. 


2. General Features 


2.1. Dirac and Schrédinger-like equations 
2.1.1. Dirac equations 
In the relativistic framework, the motion of nucleons is described by the Dirac 


equation. The corresponding eigenfunction equation for nucleons reads 
Ayr) = eW(r), (a) 


where e is the single-particle energy including the rest mass of nucleon M. 
Originating from the minimal coupling of the scalar and vector mesons to the 
nucleons in the covariant density functional theory [59-64], the single-particle 


Dirac Hamiltonian H is written as 
H=a-p+f£[M+S(r)]+ V(r). (2) 
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In this expression, œ and £ are the Dirac matrices, while S(r) and V(r) are the 
scalar and vector potentials, respectively. In addition, we set h = c = 1 in this 
paper. 

From mathematical point of view, the conclusions concerning the symmetry 
limits discussed below remain valid if either the scalar or the vector potential is 


modified by an arbitrary constant, i.e., 
Sir) > Sr) +¢s, Vir) > V@) + cy, (3) 


because one can simply adjust the mass and energy by the same constant so that 


the Dirac equation remains unchanged [67]: 
M—->M-c;, €—>€+Cy. (4) 


When the spherical symmetry is imposed, the single-particle eigenstates are 
specified by a set of quantum numbers œ = (a, Ma) = (Na, la, ja, Ma), and the 


single-particle wave functions can be factorized as 


ioe), 


r| FD- 


we (a), (5) 


with the spherical harmonics spinor DY ie) for the angular and spin parts [274]. 


The corresponding normalization condition reads 


i YTT) r = T [GU + FA] dr = 1. (6) 


Note that in this paper, to label the single-particle eigenstates we use the main 
quantum number n equal to the number of the internal nodes plus one for the 
dominant component of the Dirac spinor. Namely, the single-particle spectra start 


from the n = 1 states. 
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For the lower component of the Dirac spinor (5), one has 
GRY" (BH) =-Y! À), (7) 
with 
f=2j-1. (8) 
Thus, the single-particle wave functions can also be expressed as 


iGa(r)%, (P) | 


7 (9) 
FAY on, (f) 


1 
Wot) = = 


r 


In such a way, the pseudo-orbital angular momentum / is found to be the orbital 
angular momentum of the lower component of the Dirac spinor [65]. 


The corresponding radial Dirac equation reads 


M+X(r) -TF Il Ga(r) Gal) 

d K r = €a $ (10) 
— +2 -M+A(r) || Far) Fa(r) 

dr r 


where X(r) = S(r) + V(r) and A(r) = V(r) — S(r) denote the combinations of 
the scalar and vector potentials, and x is a good quantum number defined as x = 
+(j + 1/2) for the j = / + 1/2 orbitals. 

The SS and PSS of the Dirac Hamiltonian in Eq. (2) or (10) can be studied by 
the Bell-Reugg condition [66], see also Section 2 of Ref. [67]. Alternatively, these 
symmetries can be investigated by reducing the Dirac equation to the Schrédinger- 


like equations. 


2.1.2. Schrédinger-like equations 
Focusing on the spherical case, one can derive the Schrédinger-like equation 


for the upper component G(r) of the Dirac spinor by substituting 


1 d k 
m) = oe i sJow ie 
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in Eq. (10), and obtain 


Me E a 


1 æ 1 dM, d 1 k(k+1) 1 dM,K 
=—— + == G=6G, 
M,dr? M? dr dr M, Pr M2 


with the energy-dependent effective mass M,(r) = M — A(r) + e. For brevity 
we omit the subscripts if there is no confusion. Similarly, one can derive the 


Schrédinger-like equation for the lower component F (r) by using 


1 d K 
SS] ee 13 
oy) aa 245) ), We) 
and obtain 

1 æ 1 dM_d 1 K(k- 1) 1 dM_K 
ae e a a a GY 
| M_dP | We dr a+ | + )+ JE r? M? dr “It le 
(14) 


with the energy-dependent effective mass M_(r) = -M — X(r) + e. In Refs. [275— 
278], it has been shown that each of these two Schrédinger-like equations, 
together with its charge conjugated one, are fully equivalent to Eq. (10). 

For Eq. (12), in analogy with the Schrödinger equations, X(r) is the central 
potential in which particles move; the term proportional to (Z + 1) = x(k + 1) 
corresponds to the centrifugal barrier (CB); and the last term corresponds to the 
SO potential, which leads to the substantial SO splittings in single-particle spectra. 


Namely, 


1 xK(k+1) 1 dM,(r) k 


Ya =a Pp mE NOD a a F 


(15) 


It is well known that there is no SO splitting if the Vso vanishes. In other words, 


dM,(r) _ dA(r) | 
dr dr 


0 (16) 
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is the SS limit. 

If one uses the Schrédinger-like equation (14) for the lower component 
instead, although usually A(r) does not stand for the potential in which particles 
move, all terms except one, —(1/M?)(dM_/dr)(k/r), are identical for the pseu- 
dospin doublets a and b with Ka(Ka — 1) = Kp(Kp — 1), 1.€., Ka + Kp = 1. As pointed 
out in Refs. [68, 78], if this term vanishes, i.e., 


dM(r) dEr) _ 
"r gg Oa 


each pair of pseudospin doublets will be degenerate and the PSS will be exactly 
conserved. This is called the PSS limit, which is more general and includes the 
limit X(r) = O discussed in Ref. [65]. From the physical point of view, X(r) = 0 
is never fulfilled in realistic nuclei as in which there exist no bound states for 
nucleons [272], but dX(r)/dr ~ 0 can be approximately satisfied in exotic nuclei 
with highly diffuse potentials [69]. 

Analogically, such a term is regarded as the PSO potential, while the term 
proportional to (7 + 1) = x(x — 1) is regarded as the PCB, i.e., 


1 xKk-1) 7 1 dM_(r) k 
Min F ant VeRO "ey ae r 


Veca(7) = 


2.2. Analytical solutions at PSS limit 

Within the pseudospin symmetry limit shown in Eq. (17), the potential X(r) 
is simply a constant Xo, and then Eq. (14) for the lower component of the Dirac 
spinor can be reduced to 


d _ KK = 1) 


dr? r? 


+(e- M -oe + M- A(r))| FV) =0. (19) 


During the past decade, it is a very active field to investigate the exact or 


approximate analytical solutions of this equation within certain special forms 
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of potential A(r), by using the Nikiforov-Uvarov (NU) method [279], the 
supersymmetric quantum mechanics [252], the asymptotic iteration method 
[280], the exact quantization rule [281], and so on. For the spherical case, 
extensive investigations have been made for the spherical harmonic oscillator 
[88-90, 92, 95, 96], anharmonic oscillator [97], Coulomb [76, 99-101], Deng- 
Fan [102], diatomic molecular [103, 104], Eckart [105, 106], Hellmann [107], 
Hulthén [108-111], Manning-Rosen [112-114], Mie-type [115-117], Morse 
[118-123], Poschl-Teller [124—133], Rosen-Morse [134-137], Tietz-Hua [138], 
Woods-Saxon [139, 141, 143], and Yukawa [144-147] potentials, etc. Note that 
some of these potentials are good approximations to model the atom-atom or 
nucleon-nucleon interactions, but not nuclear mean-field potentials. 

In this Section, we will take the relativistic harmonic oscillator and relativistic 
Morse potentials as examples, and introduce the corresponding analytical solu- 
tions of equation (19) at the pseudospin symmetry limit. In the second part, the 
Nikiforov-Uvarov method [279] and the Pekeris approximation [282] for the non- 


vanishing (pseudo-)centrifugal barrier will be discussed as well. 


2.2.1. Relativistic harmonic oscillator potential 

In analogy with the Schrödinger equations with the harmonic oscillator 
potentials, the Dirac equations with the RHO potentials have received extensive 
attention in different fields of mathematics, physics, and chemistry. These 
equations have analytical solutions in many cases, for example, the corresponding 
equations at the PSS limit [88—90, 92, 95, 96]. 


By taking one of the simplest cases as an example, i.e., 


1 
A(r) = 5Mw*r and %)=0, (20) 


22 


the Schrédinger-like equation (19) for the lower component of the Dirac spinor at 
the PSS limit becomes [88, 89] 


a i+ 1) = +(€-M\e+M- TMa) F= 21) 
as K(k — 1) = (1+ 1). Note that some notations here are changed from the original 
papers for the self-consistency through the present Review, and similar changes 
have been done for the whole paper. 


By further introducing [89] 


[M -M ” 2 -MP 
y= Aor = BP and eee as (22) 
a 


the above equation is rewritten as 


d d +1) 


4p a 
Fag +25 


-5-4} F0 =0. (23) 
An asymptotic analysis suggests searching for the solutions of the type of 
FK) = Be? Dw), 24) 


where w(Ẹ) is a function to be determined and B is a normalization factor. Inserting 
this expression into Eq. (23), the equation for w(¥) reads 


es » 3 d 1 3 A 
fa Piles —j|—-=<(/4+ >+ 


5) =0. 2 
dy 2a e | A (a 


The solutions of this equation, which guarantee that lim;_,.. F,(¥) = 0, are the 


generalized Laguerre polynomials of degree ñ, ie (¥), where 


4... Ik 3.2 >, 1 
n=-5(I+3+5] and alae (26) 
Finally, one can get the eigenenergies [88, 89] 
E&x — M = 4 
(Eir + M) >M ` fon +1+ 5) ; (27) 
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which are discrete since ñ is an integer equal to or greater than zero. The 


corresponding eigenfunctions for the lower component of the Dirac spinor read 


Flr) = Be% ary ar). (28) 


Here ñ is the number of the internal nodes of F(r), denoted as np in the following 
Sections. 

The solutions of the Dirac equation with the RHO potential and their special 
features will be re-visited in Section 3.2 for the PSS in confining potentials, 
Section 3.4 for the PSS in tensor potentials, Section 3.6 for the SS in anti-nucleon 
spectra, Section 4.1 for the perturbative nature of PSS, and Section 4.2 for the 


puzzle of the intruder states. 


2.2.2. Relativistic Morse potential 


Another widely discussed potential is the relativistic Morse potential [283], 
AG\aD leer oe" |, (29) 


for atomic systems. In this expression, D > 0 is the dissociation energy, ro is the 
equilibrium internuclear distance, and a > 0 is a parameter controlling the width 
of potential well. In Refs. [118—121], the analytical solutions of the relativistic 
Morse potential at the PSS limit were investigated by using the NU method [279], 
the asymptotic iteration method [280], the exact quantization rule [281], and the 
confluent hypergeometric functions, respectively. 

In the following, we will briefly introduce one of the widely used methods, the 
Nikiforov-Uvarov method [279], and the Pekeris approximation [282] for the non- 
vanishing (pseudo-)centrifugal barrier, then discuss the solutions of the relativistic 


Morse potential [118-121]. 


24 


The Nikiforov-Uvarov method [279] is based on solving the hypergeometric- 
type second-order differential equations by means of the special orthogonal 
functions. Detailed derivations about the NU method and specific examples for 
standard Schrödinger equations with the harmonic oscillator, Coulomb, Kratzer, 
Morse, and Hulthén potentials can be found in Ref. [284]. 


The main equation which is closely associated with the NU method reads 


wor y+ O 


o(s) o? (s) 


y(s)=0, (30) 


where o(s) and G(s) are polynomials at most second-degree and 7T(s) is a first- 


degree polynomial. By letting W(s) = ¢(s)y(s) and 


p(s) _ a(s) l (31) 
p(s) o(s) 

Eq. (30) can be finally reduced into an equation of hypergeometric type, 
o(s)y"(s) + t(s)y'(s) + Ay(s) = 0, (32) 


with t(s) = T(s) + 27(s). Here both t(s) and z(s) are polynomials of degree at 
most one and J is a constant. 


The polynomial 7(s) satisfies a quadratic equation, 
n (s) + [#(s) — o'(s)la(s) + [6(s) — ko(s)] = 0, (33) 
with 
kK=A-7(s). (34) 
The corresponding solutions read 


_ o(s)~75) | (= — 7s) 


2 
a(s) = 5 5 — G(s) + ko(s). (35) 
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Since 7(s) is a polynomial of degree at most one, the expression under the square 
root has to be the square of a polynomial. In such a way, the constant k can be 
determined. In addition, the derivative of t(s) thus obtained must be negative for 
bound states. This is the main essential condition for any choice of particular 
solutions. 

To generalize the solutions of Eq. (32), it is shown that the nth-order derivative 


of y(s), v,(s) = d”y(s)/ds", is also a hypergeometric-type function, which satisfies 
(8) (S) + Tr(8)V,(5) + MnVa(s) = 0, (36) 
with 


a’ (s). (37) 


T,(S) =T(s)+no"(s) and un =A+ntT'(s)+ 


n(n — 1) 
2 


When u, = 0, Eq. (36) has a particular solution of the form y(s) = y,(s), which 


is a polynomial of degree n, and Eq. (37) becomes 


-1 
A, = -nT (s) — Dons with n=0,1,2,... (38) 


The whole set of eigenvalues for the second-order differential equation (30) can 
be obtained by comparing 4 in Eq. (34) and 4, in Eq. (38), i.e., A = Ay. 

Finally, for the corresponding eigenfunctions W(s) = @(s)y,(s), (s) is 
obtained by solving Eq. (31), and the polynomial solutions y,(s) are given by 


the Rodrigues relation, 


B, d" n 
nls) = ads [o"(s)p(s)] , (39) 


where B,, is a normalization constant and the weight function p(s) must satisfy the 


following condition: 


[o(s)p(s)!’ = T(s)p(s). (40) 
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For the relativistic Morse potential shown in Eq. (29), by assuming new 
variables œ = arọ and x = (r — rp)/To, Eq. (19) becomes 


a. _ KK 1) 
dx* (14+) 


+74(M + Xp — €)D(e* — 2e-) + r5(€ — M — Xe + M)| F(x) 


(41) 


First of all, this equation can be solved exactly by the NU method if the 
pseudo-centrifugal barrier is absent, i.e., for the p1/2 orbitals (x = 1). In this case, 
by making a new change of independent variable s = e~°*, Eq. (41) is rewritten as 


@F(s) 1dF(s 1 
ae Os Sess? este] F(s) =0, (42) 


with £; = ro(e — M -X \(e + M)/a’, & = 2r,D(M + Zp — €)/a’, and s; = 


1r,D(M + Xo — €)/a*. With reference to Eq. (30), it indicates 
T(ss=1, o(s)=s, G(s)= £35" — ES +E). (43) 


Then, one has (s) = + -~—€38? + (k + £2)s — £1, and ką = —e&) + 2 4E1863 to make 


the expression under the square root be the square of a polynomial of the first 


degree. In such a case, there are four possible forms of z(s): 


n(s) = +i(-Ve3s+ yE) for k_ =-e) -2 Ve es, 
na(s) = +i(-Ve3s— Ve,) for k, = -~-e +2 Veie3. (44) 
One of these four possible forms must be chosen to obtain the bound state 
solutions. The most suitable form can be a(s) = —i(.fé35 + y£1) for k_, as the 
derivative of t(s) thus obtained, t’(s) = —2i 4/83, is negative [118]. 
Finally, a particular solution of Eq. (42) which is a polynomial of degree n can 


be calculated from Eq. (38), 


A=A, > —&—2 VE1E3 — ive3 = 2ni VE; . (45) 
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Namely, the eigenvalues of Eq. (41) for the pı;2 orbitals can be obtained by 


(2 VD(M + Xo - 6) +2 Ve, + Me, -M - mm) +(1+2n)’a =0. (46) 


For the non-vanishing pseudo-centrifugal barrier, the Pekeris approximation 
[282] is widely used to expand the PCB about r = rọ in a series of powers of 
x = (r-1o)/To, 


K(k — 1) y 
Vecg(r) = FE = (+x)? 


= ¥(1 -—2x+ 3x --:--), (47) 


with y = x(k — 1)/rj. Here one should approximate this potential by the 


exponential forms in the Morse potential, then the approximate PCB reads 


Vecs (r) = y(Do + Die + Dre***) 


2d 4 e 
Dy + Di(1 -axs -+ D(1 -20x4 = =): 


ay 
(48) 
Comparing with these two expressions, the coefficients are 
3 3 4 6 1 3 
Do=1--+—, Di=--—, D =--+—. (49) 
a@ a a @ a @ 
In such a way, the original PCB is approximated as 
1 [oy 2 
= eS (Do + Ds + Dzs ) 3 (50) 
S a 


acting on F(s) in Eq. (42). All procedures for solving Eq. (42) remain, but simply 
substituting €; > £1 — r¢yDo/a”, €2 > £2 + yD,/a’, and £3 > €3 — rgyD2/a°. 
The final solutions of eigenvalues can be obtained by 
-2DM_ + yD i 
| +2 en + M)M_— 7s} +(1 +2n)a2 = 0, (51) 
vV-DM_ - yD, 


28 


and here M- = -M — Xo + €x. The readers are referred to Ref. [121] for 
the corresponding eigenfunctions, which can be expressed in terms of confluent 
hypergeometric function ; F. 

The same results were obtained in Refs. [119] and [120] by using the 
asymptotic iteration method [280] and the exact quantization rule [281] together 
with the Pekeris approximation, respectively. Note that the validity of the Pekeris 
approximation deserves more careful examinations. 

As a common example shown in Refs. [118-121], the parameters of the 
relativistic Morse potential in Eq. (29) are taken as D = 5.0 fm“!, rọ = 2.40873 fm, 
and a = 0.988879 fm! with the mass M = 10.0 fm™!. The corresponding 
coefficients in Eq. (49) deduced with the Pekeris approximation read Do = 
0.26928, Dı = 0.62178, and D> = 0.10893. For different choices of the constant 
Xo, one can obtain different numerical solutions for the eigenvalues of Eq. (51). 

It is instructive to choose Xp» = 0 for the comparison with the results shown 
in the coming Sections. In Table 1, the bound-state eigenenergies 6, of the 
Dirac particle in the PSS Morse potential are shown for several n and x states. 
It is confirmed that the single-particle energies of pseudospin doublets are exactly 
degenerate at the PSS limit. However, it is found that the single-particle spectra 
are bound from the top, and for each given x the single-particle energies decrease 
when the radial quantum number n increases. This indicates these bound states 
are indeed of the characteristics of the states belonging to the Dirac sea, but not 
those in the Fermi sea, even though the eigenvalues are positive and close to the 
threshold of the continuum states in the Fermi sea. It can be seen in a clearer way 
by showing the corresponding single-particle potentials explicitly in Fig. 2. Due 


to the particular shape of the Morse potential, the potential A(r) — M extends to 
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Figure 2: (Color online) Relativistic Morse potential with D = 5.0 fm™!, ro = 2.40873 fm, a = 
0.988879 fm™!, and M = 10.0 fm~! corresponding to the PSS limit £(r) = Xo = 0. 
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Table 1: The bound-state eigenenergies €x of the Dirac particle in the pseudospin-symmetry 


Morse potential with the Pekeris approximation and 2%» = 0. Energy units are in fm™. 


data are taken from Ref. [119]. 


n 


1 


N NY N N 


(L j) 
51/2 
P3/2 
ds 
fir 
51/2 
P3/2 
ds 
fin 


Enk 
9.9935101 
9.9838165 
9.9737712 
9.9656754 
9.9929544 
9.9807043 
9.9652868 
9.9484873 


(j) 
dz 
fsi2 
87/2 
hoj2 
d3)2 
fn 
87/2 


ho)2 


Enk 
9.9935101 
9.9838165 
9.9737712 
9.9656754 
9.9929544 
9.9807043 
9.9652868 
9.9484873 


1 


The 


the positive-energy region and forms a pocket for the bound states, which belong 


to the Dirac sea. 


Furthermore, when Xo is chosen as —M instead of 0, the general pattern of the 


single-particle spectra does not change [118—121], since it is a trivial modification 


by a constant as shown in Eqs. (3) and (4). In this case, it is clear that all solutions 


of Eq. (51) are of negative energy and there are no bound states with positive 


energy. 


Therefore, more precisely, the phenomenon discussed in Section 2.2.2 corre- 


sponds to the spin symmetry of the anti-particle spectra, which will be discussed 


in more detail in Section 3.6. 


31 


2.3. PSS breaking in realistic nuclei 


For the isolated atomic nuclei, both the scalar S(r) and vector V(r) potentials 
in the single-particle Dirac Hamiltonian H in Eq. (2) vanish at large distance from 
the center, i.e., lim, (r) = lim,. A(r) = 0. This kind of potentials will be 
specified as the non-confining potentials in Section 3.2. Within these potentials, 
it can be proven that there are no single-particle bound states in the Fermi sea 
at the pseudospin symmetry limit shown in Eq. (17) [272]. In other words, the 
pseudospin symmetry must be broken in realistic nuclei. 

Therefore, much more meaningful and important tasks are to investigate to 
which extent the pseudospin symmetry is approximately conserved in realis- 
tic nuclei and what the symmetry-breaking mechanism is. This approximate 
pseudospin symmetry can be examined by the quasi-degenerate single-particle 
energies, as well as by the relation between the pseudospin-orbit potential and 
pseudo-centrifugal barrier pointed out by Meng et al. [68, 69] and the relations 
between the single-particle wave functions pointed out by Ginocchio and Madland 
[79, 80]. 

As the pseudospin symmetry is shown to be a relativistic symmetry [65], let 
us start with the introduction of the covariant density functional theory [59-64], 
which is one of the most appropriate microscopic and self-consistent approaches 


for studying the properties of pseudospin symmetry in realistic nuclei. 


2.3.1. Covariant density functional theory 

The CDFT [59-64] can be traced back to the successful RMF models 
introduced by Walecka and Serot [285]. The most popular RMF models are 
based on the finite-range meson-exchange representation, in which the nucleus 


is described as a system of Dirac nucleons that interact with each other via 
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the exchange of mesons and photons. The nucleons and mesons are described 
by the Dirac and Klein-Gordon equations, respectively. Together with the 
electromagnetic field, the isoscalar-scalar o meson, the isoscalar-vector w meson, 
and the isovector-vector o meson build the minimal set of meson fields that 
is necessary for a description of bulk and single-particle nuclear properties. 
Moreover, a quantitative treatment of nuclear matter and finite nuclei needs a 
medium dependence of effective mean-field interactions, which can be introduced 
either by including nonlinear meson self-interaction terms in the Lagrangian 
[286-290] or by assuming explicit density dependence for the meson-nucleon 
couplings [290-293]. 

The Lagrangian density of the RMF theory with nonlinear meson self- 


interactions [287—290] can be written by using the conventions in Ref. [61] as 


-T 
3 Aaly 


: 1 
L =f | iy"ðu — M - 80T — 8oy Wy — 8Y"? -Bu — ey" 


1 1 lsa 1 
+ A -— U, (T) — g w + U,(@,) — a Ry + U,(B,) - re 


pv? 


(52) 


where M and m; (g;) (i = o@,w,p) are the masses (coupling constants) of the 


nucleon and mesons, respectively, and 
W = Ho" - Ow", RP = Op’ — OF, FP = OKAY - OA", (53) 


are the field tensors of the vector mesons and electromagnetic field. We adopt the 
arrows to indicate vectors in isospin space and bold type for the space vectors. 
Greek indices u and v run over 0, 1, 2,3 or t, x, y, z, while Roman indices i, j, etc. 
denote the spatial components. The nonlinear self-coupling terms U,(~), U,,(w,), 


and U,((,,) for the a, w, and p mesons in the Lagrangian density (52) respectively 
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have the following forms: 
1 1 1 
U (T) = z"o + z8r + 7e 
1 1 
U.(Wy) = so Oy + gE O) : 
> 1 = 1 2? 
U Ba) = 5P But GBP PY (54) 


The system Hamiltonian density can be obtained via the general Legendre 


transformation, 
Of . 
KH = — p- L, 55 
id; $ (55) 


where ¢; represent the nucleon-, meson-, and photon-field operators. The ground- 


state trial wave function is taken as a Slater determinant, 
A 
Io) =| [et +), (56) 
a=1 


where |—) is the physical vacuum and the single-particle states œ are confined 
to those with positive energies in the Fermi sea, i.e., the no-sea approximation. 
Combining these two expressions together, one has the energy density functional 


for the whole system, 
Epor = (®o| Z |®o) . (57) 


In the RMF theory, only the direct contributions of the meson and Coulomb fields, 
i.e., the so-called Hartree terms, are taken into account. 

For the systems with time-reversal symmetry, the space-like components of 
the vector fields vanish. Furthermore, one can assume that the nucleon single- 
particle states do not mix isospin, i.e., the single-particle states are eigenstates of 


T3, therefore only the third component of 8, survives. 


34 


By the variation principle, the Dirac equation for nucleons reads [cf. Eqs. (1) 


and (2)] 
[a-p+B(M + S(r)) + VO)| Wo) = Coal), (58) 
and the Klein-Gordon equations for mesons and photons read 
-V’o + Uj(0) = -8ops » 
—V*wo + U, (wo) = SuPv> 
-V° p + U00) = SPV’ 
-V°Ao = epc. (59) 
The scalar and vector potentials in Eq. (58) are, respectively, 


S) = gor), 


1-7; 


Vr) = guWo(t) + Sptspp(r) + eAo(r). (60) 


The scalar density ps, the baryonic density py, the isovector density P, and the 
charge density pc in the Klein-Gordon equations (59) are, respectively, 


A 
ps (r) = X Paal), 


a=1 


A 
pytt) = X Wi@Welr), 


a=1 


A 
pp Œ) = X Eyal), 
a=1 


Z 
pelt) = X HO). (61) 


p=1 
In the density-dependent RMF approach [290-293], the nonlinear meson self- 


couplings in the Lagrangian density are replaced by the density dependence of 
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the coupling strengths g,(p), g.(e), and g,(p), and an additional term, i.e., the 
rearrangement term, will appear in the Dirac equation (58). 

More recently, this framework has been re-interpreted by the relativistic Kohn- 
Sham density functional theory, and the functionals have been developed based 
on the zero-range point-coupling interaction [294—298], in which the meson 
exchange in each channel (isoscalar-scalar, isoscalar-vector, isovector-scalar, 
and isovector-vector) is replaced by the corresponding local four-point contact 
interaction between nucleons. The point-coupling model has attracted more and 
more attention owing to the following advantages. First, it avoids the possible 
physical constrains introduced by explicit usage of the Klein-Gordon equation 
to describe mean meson fields, especially the fictitious o meson. Second, it is 
possible to study the role of naturalness [299, 300] in effective theories for nuclear- 
structure-related problems. Third, it is relatively easy to include the Fock terms 
[301], and provides more opportunities to investigate its relationship to the non- 
relativistic approaches [302]. 

In order to describe open-shell nuclei, the pairing correlation and the coupling 
to continuum must be taken into account properly, which are in particular 
crucial for the descriptions of drip line nuclei. The extension of the RMF 
theory to take into account both bound states and (discretized) continuum via 
Bogoliubov transformation in a microscopic and self-consistent way has been 
done in Refs. [205, 207], the so-called relativistic continuum Hartree-Bogoliubov 


(RCHB) theory. 


2.3.2. PSS in single-particle energies 
First of all, let us start with the neutron-rich doubly magic nucleus '°?Sn as an 


example, which shows that the self-consistent CDFT or RMF theory can nicely 
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Figure 3: (Color online) Self-consistent single-particle potentials for neutrons in '**Sn calculated 


by the RMF theory with the effective interaction PK1 [290]. 


reproduce its ground-state properties including the single-particle spectra [303]. 

In Fig. 3, the potentials X(r) and A(r) for neutrons calculated by the RMF 
theory with the effective interaction PK1 [290] are shown. The depths of potentials 
are X(r) ~ 70 MeV and A(r) ~ 700 MeV, respectively. 

The single-particle energies of the neutron bound states thus obtained are 
shown in Fig. 4, where E = e — M excluding the rest mass of nucleon. In 
order to show the SO and PSO splittings and to see their energy dependence more 
clearly, the reduced SO splittings AEso = (E;.—E;,)/(2/+ 1) and the reduced PSO 
splittings AEpso = (E; - E; )/ (2Ï+ 1) versus their average single-particle energies 
Ew = (Ej. + E;,)/2 are plotted in the left and right panels of Fig. 5, respectively. 
In this paper, j< (J>) denotes the states with j = /—- 1/2 (j = L+ 1/2) for the spin 
doublets and the states with j = 7- 1/2 (j = Ï+ 1/2) for the pseudospin doublets. 
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E (MeV) 


Figure 4: (Color online) Neutron single-particle spectrum in !??Sn calculated by the RMF theory 
with PK1. For each pair of spin doublets, the left levels are those with j< = /— 1/2 and the right 


ones with j> = l + 1/2. Potential X(r) is shown as the solid line. 
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Figure 5: (Color online) Reduced SO splittings (Ej. — E; )/(2l + 1) (the left panel) and reduced 
PSO splittings (Ej. — Ej,)/ (27 + 1) (the right panel) versus their average single-particle energies 


Ey = (Ej, + Ej,)/2 in neutron single-particle spectrum of 132Sn calculated by the RMF theory 
with PK1. 
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A dramatic energy dependence can be seen in the reduced PSO splittings, 
whereas the reduced SO splittings are less energy dependent. While the reduced 
PSO splitting for the 1p pseudospin doublets is 1.291 MeV, that for the 2p 
doublets is 0.138 MeV, roughly smaller than the former one by a factor of 10. 
Thus, the PSS becomes better near the Fermi surface. This is in agreement with 
the experimental observation. 

Around the Fermi surface, E ~ —8 MeV in this case, AEpsg = 0.138 MeV 
for the 2) pseudospin doublets, compared to AEso = 0.332 MeV for the 2d spin 
doublets. Note that these two cases in the comparison share a common state 2d3/9. 
Further approaching the single-particle threshold, on one hand, the reduced SO 
splitting for the 1h doublets is almost the same as those for the 1g and 1 f doublets 
below the Fermi surface; on the other hand, the PSO splittings become smaller 
and even reversed, e.g., AEpso = —0.068 MeV for the 1g doublets. 

To understand why the energy splitting between the pseudospin partners 
changes with the single-particle energies, the PSO potential Vpso and the pseudo- 
centrifugal barrier Vpcg in Eq. (18) should be examined carefully [68, 69]. Their 
contribution to the single-particle energy E can be evaluated by the integrals with 


the lower component of the Dirac spinor, i.e., 


T LROVOFOdr (62) 
f Fe@)F@)dr 
It was pointed out in Refs. [68, 69] that 
|Veso(r)| < |Vecs(7)| (63) 


is the condition under which the PSS is conserved approximately. Unfortunately, it 


is difficult to plot and compare these potentials, as both of them have a singularity 
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at ro where M_(r) = 0, i.e., E = X(r)|,-<,,. AS one is only interested in the relative 


magnitude of the PCB and the PSO potential, the effective PCB, 

Veen (") = M-K(K -= 1/7, (64) 
and the effective PSO potential, 

Vigor) = (-dM_/dr)(k/r), (65) 


are introduced for comparison [68, 69]. They correspond respectively to the PCB 
and the PSO potential multiplied by a common factor M2. 


5.0 
4.0 
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Centrifugal & PSO Potential '°°Zr [fm°] 


Figure 6: Comparison of the effective PCB (dashed and dot-dashed lines) and the effective PSO 
potential (solid line) in arbitrary scale for d3/2 (the upper panel) and s1/2 (the lower panel) states 
in "Zr. The dashed lines are for 1d3 j2 and 21/2 states, and the dot-dashed lines are for 2d3/2 and 
35/2 states. The inserted boxes show the behaviors of the effective PCB and the effective PSO 


potential near the nuclear surface. Taken from Ref. [69]. 


The effective PSO potential in Eq. (65) depends on the angular momentum 


and parity, but does not depend on the single-particle energy. On the other hand, 
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the effective PCB in Eq. (64) depends on the energy. They are given in Fig. 6 for 
the sı;2 and d3/. orbitals of 1202r in arbitrary scale, and their behavior near the 
nuclear surface are shown in the inserts. Note that the solid lines in the figure are 
the effective PSO potentials, which are enlarged in the inserts. It is found that the 
PSS is conserved much better for the less bound pseudospin partners, because the 


effective PCB is smaller for the more deeply bound states. 
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Figure 7: Comparison of the effective PCB (dashed lines) and the effective PSO potentials (solid 


lines) multiplied by F 2(r) in arbitrary scale for 51/2 (left panels) and d3/2 (right panels) states in 


1202r. Taken from Ref. [68]. 


To see clearly their contributions to the single-particle energy, the effective 
PCB and the effective PSO potential multiplied by the squares of the lower 
component wave function F(r) are given in Fig. 7 for the 25)/2, 351;2, 1d3/2, 
and 2d3/ states of "Zr in arbitrary scale. It is clear that the contributions of 
the effective PCB are much bigger than those of the effective PSO potential, and 


generally they differ by two orders of magnitude. In a semi-quantitative sense, 


this indicates the condition in Eq. (63) is satisfied. 


4] 


However, it was also pointed out in Ref. [167] that, from a more quantitative 
point of view, one should directly compare the PSO potential and PCB, instead 
of using the effective ones, because the common factor M? multiplied depends 
on r. As a result, the magnitude of the PSO potential is drastically modified 
around rọ, and the inequality |M2(r)Vpso(7)| < |M2(r)Vece(r)| differ from 
\Vpso(‘)| < |Vpcg(r)| shown in Eq. (63). Furthermore, although Vpso(r) and 
Vpcg(r) have a singularity, it has been proven that the principal values of the 
integrals, P f F*(r)V,(r)F(r)dr, are still finite due to the nodal structure of F(r) 
[72, 167]. This makes a direct comparison possible. Several examples will be 


shown in Sections 3.3, 3.4, 3.7, and 4.1. 


2.3.3. PSS in single-particle wave functions 
There are intensive discussions on the approximate energy degeneracy be- 
tween the pseudospin doublets since the introduction of PSS in 1969, but less 
on their wave functions until the relativistic origin of PSS was revealed [79, 80]. 
Within the PSS limit shown in Eq. (17), the Schrédinger-like equation for the 
lower component F'(r) of the Dirac spinor is expressed as 


d _ K(k-1) 


dr? r2 


+(e- M -oe + M- A(r)|F(r)=0, (66) 


as seen in Eq. (19). It is clear that this equation is identical for the pseudospin 
doublets a and b with Ka(Ka — 1) = Kp(Kp — 1), i.e., Ka + Kp = 1. Therefore, the 
eigenfunctions F,(r) and F;,(r) are exactly the same up to a normalization factor. 
Moreover, there holds Eq. (13), 
1 d K 
G(r) = ———— |-— + -| F (r), 67 
”) | dr ”) (on) 


and now Xp is just a common constant. Therefore, for the single-particle wave 
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functions of the state a, the normalization condition in Eq. (6) reads 


2 2 = 1 4 Ka 
fione: feitir 


2 
+ ro dr=1. 


(68) 
Integrating by parts and using Eq. (66), one will end up with 
{ (2e, — A(r) — Xo) F2(r)dr = & — M- Xo. (69) 
In the same way, for its pseudospin partner b, one has 
1 (26, — A(r) — Xo) F? (rdr = & — M- Xo. (70) 


AS €, = & at the PSS limit, the normalization factors for F,(r) and F;,(r) are the 
same. Therefore, for a pair of pseudospin doublets, their lower components of the 


Dirac spinor are identical (up to a phase) at the PSS limit [79], 
F,(r) = Fy). (71) 


As a step further, together with Eq. (11), the first-order differential relation for 


their upper components of the Dirac spinor can be obtained [80], 


d Ka [d k 
(5 + <) G,(r) = (5 + 2) G(r), (72) 
or written as 
d I d i+1 
Gaa a i 


with j (j>) labelling the j = Ï- 1/2 (j = Ï+ 1/2) orbital. 
The wave-function relations (71) and (72) between the pseudospin doublets 


can also be derived from the pseudospin SU(2) generator [304] 


7 S 0 
Š = ; (74) 
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Figure 8: (Color online) Neutron single-particle wave functions of the 1p, 2p, and 1 f pseudospin 


doublets in !?Sn calculated by the RMF theory with the effective interaction PC-PK1 [297]. 


with s = 0/2 and § = (ø - p)s(o - p). The details can be found in Ref. [80]. 

To test the relation shown in Eq. (71), the single-particle wave functions for 
the neutrons in "Sn are calculated by the self-consistent point-coupling RMF 
theory with the effective interaction PC-PK1 [297]. In panels (a), (b), and (c) of 
Fig. 8 are shown the wave functions of the 1), 2p, and 1f pseudospin doublets, 
respectively. For each pair of pseudospin doublets, their upper components G 
have different number of nodes and radial shape, however, their lower components 
F are almost identical except on the nuclear surface. By comparing these three 
panels, it is found that the relation in Eq. (71) is better satisfied for smaller /. 

To test the differential relation of the upper components shown in Eq. (72), 
the corresponding results obtained by using the first-order differential operators 
are plotted in Fig. 9. It is found that, with the x-dependent first-order differential 
operators, one obtains a remarkable similarity in the differential wave functions 


except near the nuclear surface. By comparing three panels, it is also found that 
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Figure 9: (Color online) Differential relation of the upper components G in Eq. (72) for the 1), 


2p, and 1f pseudospin doublets in '*?Sn calculated by the RMF theory with PC-PK1. 


the relation in Eq. (72) is better satisfied with smaller /. 

In Section 2, the general features for the Dirac equation and its corresponding 
Schrédinger-like equations were discussed. The analytical solutions for Dirac 
equation at the pseudospin symmetry limit were shown by taking the relativistic 
harmonic oscillator and relativistic Morse potentials as examples. The pseudospin 
symmetry and its breaking in realistic nuclei were discussed in a general 
framework of the covariant density functional theory. The evaluations of the 


pseudospin symmetry in the single-particle energies and wave functions were 
reviewed. 
3. PSS and SS in Various Systems and Potentials 


3.1. From stable nuclei to exotic nuclei 


The concept of pseudospin symmetry [11, 12] was introduced originally based 


on the observation of the empirical single-particle spectra in stable nuclei. Since 
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then, intensive discussions of PSS were mainly concentrated around the nuclear 
-stability valley. During the past decades, more and more highly unstable nuclei 
with extreme N/Z ratios have been accessible with the radioactive ion beam 
facilities. The physics connected to the extreme neutron richness in these nuclei 
and the low density in the tails of their matter distributions have attracted much 
attention, and new exciting discoveries have been made by exploring hitherto 
inaccessible regions in the nuclear chart. One of the examples is the investigations 
of PSS from stable to exotic nuclei [68, 69]. 

From the theoretical point of view, for open-shell exotic nuclei, the RCHB 
theory [61] is able to take the pairing correlation and the coupling to continuum 
into account properly. Furthermore, as pointed out in Ref. [206], in the RCHB 
theory, the particle levels for the bound states in the canonical basis are the same 
as those by solving the Dirac equation with the corresponding scalar and vector 
potentials. Therefore, the Schrédinger-like equations (12) and (14) remain the 
same in the canonical basis even after the pairing interaction has been taken into 
account. 

One observes from Eq. (14) that formally the only term which breaks the PSS 
is the PSO potential Vpso (18), which is proportional to dX(r)/dr. Therefore, it 
is expected that the PSS is better conserved when |dX(r)/dr| becomes small [68]. 
This conjecture can be verified with the exotic nuclei, whose potentials can be 
much more diffuse than the stable ones [69]. 

In the left panel of Fig. 10, the potentials Ł(r) for neutrons in Sn isotopes 
calculated by the self-consistent RCHB theory with the effective interaction 
NLSH [287] are shown. One can see a gradual change in the diffuseness of X(r) 


from the neutron-deficient nucleus '°°Sn to the extremely neutron-rich nucleus 
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Figure 10: Left panel: Neutron potentials X(r) for Sn isotopes calculated by the RCHB theory [61] 
with the effective interaction NLSH [287]. Taken from Ref. [305]. Right panel: The corresponding 
neutron single-particle energies in the canonical basis as a function of the mass number, where the 
pseudospin doublets are marked by boxes and the Fermi surface is shown by a dashed line. Taken 


from Ref. [69]. 


"Sn. The corresponding evolution of the single-particle energies in the canonical 
basis can be seen in the right panel of Fig. 10. 

The reduced PSO splittings AEpso versus their average single-particle energies 
are plotted in Fig. 11. It is seen that the PSO splittings in Sn isotopes have a 
monotonous decreasing behavior with increasing isospin. In particular, for the 1p 
doublets, AEpso in '”°Sn is only half of that in °°Sn. Furthermore, a monotonous 
decreasing behavior of AEpso with increasing single-particle energies E maintains 
from the proton drip line to the neutron drip line. 

From these studies, the pseudospin symmetry remains a good approximation 
for both stable and exotic nuclei. A better pseudospin symmetry can be expected 
for the orbitals near the threshold, in particular for nuclei near the particle drip 


line. 
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Figure 11: Reduced PSO splittings for neutrons in Sn isotopes versus their average single-particle 


energies [7+ DE jet IE l (21+ 1) calculated by RCHB theory with NLSH. Taken from Ref. [69]. 


3.2. From non-confining potentials to confining potentials 

It is observed that the main quantum numbers n of a pair of pseudospin 
doublets differ by one, e.g., (251/2, 1d3/2), (2p3/2, 1fs/2), etc. This phenomenon 
motivated Leviatan and Ginocchio [272] for the analytical proof on the nodal 
structure of the Dirac spinor. For the so-called non-confining potentials, which 
mean S(r), V(r) — 0 for r > ov, it is proven that the number of internal nodes of 


the upper G and lower F components of the Dirac spinor, ng and nr, obeys [272] 
nr=ng for «x<0O, Nr =not1 for Kx>Q. (75) 


It is also proven that there exist no bound states in the Fermi sea at the pseudospin 
symmetry limit (17) within the non-confining potentials. 
In contrast to the non-confining potentials, Chen et al. [88] showed there 


exist bound states in the Fermi sea at the pseudospin symmetry limit (17) when 
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the potential A(r) is confining. A typical example is the relativistic harmonic 
oscillator potential [88—90, 92]. Recently, the corresponding nodal structure of the 
Dirac spinor within the confining potentials was derived analytically by Alberto, 
de Castro, and Malheiro [273]. 

In this Section, we will highlight the key steps of these analytical proofs. 
We will then discuss the single-particle spectra and wave functions given by 
the relativistic harmonic oscillator potential, which conserves the pseudospin 


symmetry exactly. 


3.2.1. Nodal structure for non-confining potentials 

First of all, let us focus on the so-called non-confining potentials, which 
generally occur in isolated atomic nuclei. Their scalar and vector potentials satisfy 
rS(r),rV(r) —> 0 for r —> 0, and S(r), V(r) — 0 for r — œ [272]. A typical 


example for such potentials calculated by the self-consistent RMF theory is shown 


in Fig. 3. 
From the radial Dirac equations (10), it is seen that the radial wave functions 
follow 
PG(r) 
= ae Oe) (M? - &)G(r), 
@ F(r) 
Fa YM OM_O)FO) ~ M? -= eFC), (16) 


at large r. As for a bound state, the wave functions G(r) and F (r) should vanish 
exponentially, G(r) œ F(r) œ e~” with A = VM? — e. Therefore, -M < € < M 
is the condition for the bound states. 

Focusing on the single-particle bound states in the Fermi sea, € > 0, the 


effective mass M,(r) is always positive, and the effective mass M-(r) is positive at 
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the origin and becomes negative at large r, changing its sign at ro. The asymptotic 


behaviors of their radial wave functions at r > 0 read 


F M_(O 
Gar“, Faar", lim a =O <0 for «<0, 
G M,(0 
Gdne«r*!, Far)“, lim ce ee A >0 for x>0O. (77) 


390 F(r) 2K +1 

In order to study the properties of the radial wave functions, it is helpful 

to introduce G,(r) = r“G,(r) and F,(r) = r“F,(r), and then Eqs. (10) can be 
simplified as [272] 


Ae a = PM, (Fr) and 


dF Ar) _ 
d dr 


-r ™M_(r)GKr). (78) 
In the open interval (0, co), the nodes of G and F coincide with the nodes of G 
and F, respectively. 

Equations (78) lead to a number of observations [272]. First, it is impossible 
for G and F, or G and F, to vanish simultaneously at the same point because, 
if they did, then all other higher-order derivatives would vanish at that point and 
hence the functions themselves would vanish everywhere. Moreover, a node of F 
corresponds to an extremum of G, and a node of G corresponds to an extremum 
of F. Since M_(r) changes sign at ro, F can have an additional extremum at this 
point, which does not correspond to a node of G. It follows that the nodes of F 
and G alternate, i.e., between every pair of adjacent nodes of F (G) there is one 
node of G (F). 

Furthermore, for bound states, both G(r) and F (r) vanish at r = œ and their 
extrema are concave towards the r-axis. Therefore, the nodes of G and F can 
occur only where M,(r)M_(r) > 0, i.e., both M,(r) > 0 and M_(r) > O in the 


present cases. This is consistent with the non-relativistic case that the nodes of the 
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radial wave function can occur only in the region of classically allowed motion, 
where the kinetic energy is positive. 
One can now use the above results to obtain a relation between the radial nodes 


of G(r) and F(r), together with 
GNF) = (GNF (NY = M (FF) — MGA). (79) 


For the single-particle bound states in the Fermi sea within the non-confining 
potentials, (GFY = (e + M)F? — (e - M)G* > O at large r. At small r, 
(GFY = —M_(0)G* < 0 for x < 0, while (GF) = M,(0)F? > 0 fork > 0. 


Since GF vanishes at both r = 0 and r = ov, one confirms that 
lim G,rF.d(r)<0O for «<0, lim G.r)F.(r)>0 for «x>0, (80) 
and 
lim G.Ar)F (r) <0. (81) 


Furthermore, since both M,(r) and M_(r) are positive at nodes of G and F, one 


has 

IGF) = -M-r1)G (r) < 0, (82) 
at F(r,) = 0, and 

[GFT en = M72) F°(r2) > 0, (83) 


at G(r2) = 0. This indicates G,(r)F;,(r) is a decreasing function at the nodes of F, 
and an increasing function at the nodes of G. 
Combining all the properties discussed above, one can conclude that: (i) For 


K > 0, G,(r)F,(r) has an odd number of zeroes, and the first and the last zeroes 
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belong to the nodes of F. As a result, F has one more node than G. (ii) For x < 0, 
G,(r)F,(r) has an even number of zeroes, and the first and the last zeroes belong 
to the nodes of G and F, respectively. As a result, G and F have the same number 
of nodes. Note that this also includes the case that neither G nor F has internal 
nodes. Therefore, the number of internal nodes of G and F, ng and nr, obeys 


[272] 


nr=ng for «<0, nr=ng+1 for Kx>Q. (84) 


Figure 12: Products of the upper and lower components G(r)F (r) of the neutron 31/2, 2d3/2, and 


1ds/2 states in '°?Sn calculated by the RMF theory with the effective interaction PC-PK1 [297]. 


In Fig. 12, the products of the upper and lower components G(r)F(r) are 
shown for the 35/2, 2d3/2, and lds); states. From panels (a) and (b), one can 
confirm all the properties discussed above. If there exist internal nodes, the last 
one must belong to F. For the x > 0 states, F has one node more than G, while 
for the x < 0 states, F has the same number of nodes as G. Meanwhile, G(r)F(r) 
is a decreasing function at the nodes of F and an increasing function at the nodes 


of G. 
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It should be noted that there is a kind of special case that np = ng = 0 for all 
the lowest j = / + 1/2 orbitals, e.g., the 1d5;2 state shown in panel (c) of Fig. 12. 
These orbitals are special because they have no pseudospin partners. We will 
come back to this point with the discussions on the puzzle of intruder states in 
Section 4.2. 

Before ending this part, it is important to point out that, within the non- 
confining potentials, there exist no bound states in the Fermi sea at the PSS limit 
[272]. The reason is following: On one hand, for the bound states, M_(r) must 
be negative at large r because M,(r) is always positive. The exact PSS limit (17) 
means M_ is simply a constant, which means M_ is a negative constant here. On 
the other hand, from Eqs. (80) and (81), we know that G(r)F(r) goes to zero for 
both small and large r, but G(r)F(r) cannot be identically zero for all r. This 
means [G(r)F(r)]’ in Eq. (79) must be negative for some range of r and positive 
for some other range of r. However, this is impossible if M,(r) is always positive 
and M_(r) is a negative constant. This contradiction demonstrates that within the 
non-confining potentials there exist no bound states in the Fermi sea at the PSS 
limit. 

3.2.2. Nodal structure for confining potentials 


Chen ef al. pointed out in Ref. [88] that there can exist bound states in the 


Fermi sea at the PSS limit (17), dx(r)/dr = 0, if the potential A(r) satisfies 


lim Air) > M+e>2M+2Xo. (85) 


roo 


As discussed in the previous part, for a bound state, its wave functions G(r) and 
F(r) should go to zero exponentially at large r. From Eqs. (76), one finds that this 
requirement can be fulfilled either by M,(r) > 0 and M_(r) < 0 or by M,(r) < 0 
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and M_(r) > 0 at large r. The bound states in the non-confining potentials belong 
to the former case. On the other hand, the exact PSS limit (17) indicates M_ is 
simply a constant, in particular, it should be positive for all possible solutions 
belonging to the Fermi sea. These solutions become the bound states as long as 


M,(r) < 0 at large r, which corresponds to the condition given in Eq. (85). 
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Figure 13: (Color online) Relativistic harmonic oscillator potential for V(r) = -S (r) = Mor 


with M = 10.0 fm™! and w = 2.0 fm™!. 


A typical example for such kind of potential is the relativistic harmonic 


oscillator potential [88—90, 92], 
l 2,2 
A(r) = zvu r, (86) 


as illustrated in Fig. 13. This is called the confining potential [273] as A(r) — œ 


for r — oo, and it is different from the non-confining potential as shown in Fig. 3. 
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Recently, the nodal structure of the wave functions obtained within the 
confining potential was demonstrated in an analytical way by Alberto, de Castro, 
and Malheiro [273]. Assuming the potential A(r) goes to infinity at large r as 
a power law, lim, A(r) = Cr*, with C > 0 anda > 0, the wave functions 
G(r) and F(r) go to zero exponentially as G(r) œ F(r) œ e*” "Tn this case, 
A = 2 VCle— M —Xp)/(a + 2) indicating that € > M + Xp is the condition for the 
bound states. 

For these single-particle bound states, the effective mass M_(r) is a positive 
constant, and the effective mass M,(r) is positive at the origin and becomes 
negative at large r, changing its sign at rọ. Therefore, the asymptotic behaviors of 
their radial wave functions at small r still have the properties shown in Eqs. (77). 
In addition, Eqs. (79), (80), (82), and (83) also remain valid. It follows that nodes 
of G and F can occur only where both M,(r) > 0 and M_(r) > 0, as well as 
G(r)F(r) is a decreasing function at the nodes of F and an increasing function at 
the nodes of G. 

Compared with the non-confining potentials, the difference happens in Eq. (81). 
Within the confining potentials, M_(r) > 0 and M,(r) < 0 at large r, and thus the 


right hand side of Eq. (79) is negative. This indicates 
lim G,(r)F,(r) > 0, (87) 


and thus if there exist internal nodes in the wave function the last one must belong 


to G. Therefore, the number of internal nodes of G and F, ng and np, obeys [273] 


nr =ng-1 for «<0, Nnr=ng for x>O. (88) 
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3.2.3. Exact PSS in confining potentials 
As pointed out above, the exact PSS can be fulfilled for the single-particle 
bound states in the Fermi sea within the confining potentials. The first example 


was shown in Ref. [88]. 
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Figure 14: Single-particle spectrum in units of fm~! for V(r) = —S(r) = Mu*r*/4 (the left 
panel) and V(r) = -S (r) = Mw?r?/4 — Vo/[1 + e°-®/] (the right panel) with M = 10.0 fm}, 
w = 2.0 fm, Vo = 5.0 fm™!, R = 0.3 fm, and a = 0.05 fm. The pseudospin singlets are shown 
with the dashed line, and the pseudospin doublets with the solid lines. Taken from Ref. [88] and 


modified to present notations. 


In the left panel of Fig. 14, the single-particle spectrum for V(r) = —S (r) = 
Mw’r*/4 with M = 10.0 fm! and w = 2.0 fm™! is shown. Note that different 
from the figures in Ref. [88], here the main quantum number n of each state equals 
the number of its internal nodes for the dominant component of the Dirac spinor 
plus one. It can be seen that there are the exact pseudospin degeneracies for the 
pseudospin doublets, such as (21/2, 1d3/2), (2p3/2, 1fs/2), etc. This spectrum has 
higher degeneracy than that of PSS due to the speciality of the harmonic oscillator 
potential. To eliminate the extra degeneracies, an additional term of Woods-Saxon 


potential is introduced, V(r) = —S(r) = Mw*r?/4 — Vo/[1 + e°-®/] with Vo = 
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5.0 fm™!, R = 0.3 fm, and a = 0.05 fm. The single-particle spectrum thus obtained 
is shown in the right panel of Fig. 14, where the pseudospin singlets are shown 
with the dashed lines, and the doublets with the solid lines. By comparing these 
two panels, we see this new term reduces all redundant degeneracies of harmonic 
oscillator potential and keeps the pseudospin symmetry. 

The radial wave functions for the pseudospin singlets, 1pı;2 and 2p,/2, are 
shown in Fig. 15(a) and (b), where the dashed lines are for the lower components 
Fx of the Dirac spinor and the solid lines for the upper components G,,. Those 
for the pseudospin doublets, (2.1/2, 1d3/2), (3.81/2, 2d3/2), (2p3/2, 1fs/2), and (3p3/2, 
2f5/2), are shown in Fig. 15(c)-(e). It can be seen that for every pair of doublets, 
their lower components are identical as proven before Eq. (71). Meanwhile, their 
upper components look very different since they obey the first-order differential 
relation shown in Eq. (72). 

Furthermore, Fig. 15 confirms that if there exist internal nodes in the radial 
wave functions, the last one belongs to the upper components G. It is also seen 
that the upper and lower components are now in phase at large r, rather than out 
of phase for the cases of the non-confining potentials shown in Figs. 8 and 12. 

It is interesting and important to point out that within the confining potentials 
all of the n = 1 eigenstates for the j> = l+ 1/2 (x < 0) orbitals disappear in the 
Fermi sea, such as 1581/2, 1p3/2, 1ds/2, etc., as seen in Fig. 14. This is because the 
nodal structure (88) of the wave functions obtained within the confining potentials 
is different from the conventional nodal structure (84) appearing in the isolated 
atomic nuclei. For the x < 0 orbitals, it has been proven that np = ng — 1, and 
since np cannot be negative, all of such eigenstates have at least one internal node 


in their upper components of the Dirac spinor. In other words, the single-particle 


57 


Figure 15: Unnormalized radial wave functions for V(r) = —S(r) = Mw*r?/4 — Vo/[1 + e’-®/4] 
with M = 10.0 fm™!, w = 2.0 fm™!, Vo = 5.0 fm™!, R = 0.3 fm, and a = 0.05 fm. Pseudospin 
singlet states: (a) 1pi/2 and (b) 2pi/2. Pseudospin doublet states: (c) (251/2, 1d3/2), (d) (381/2, 
2d3/2), (e€) (2p3/2, 1fs2), and (f) (3p3/2, 2f5/2). Taken from Ref. [88] and modified to present 


notations. 
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spectra for x < 0 orbitals start from the main quantum number n = 2. As a result, 
there are no intruder states in the single-particle spectrum at all, and every / > 0 
state has its own pseudospin partner. 

In summary, the nodal structure of the Dirac spinor can be derived analytically. 
For the non-confining and confining potentials, the corresponding relations are 
shown in Eq. (84) and Eq. (88), respectively. This also explains the reason why at 
the pseudospin symmetry limit there are no bound states in the Fermi sea within 
the non-confining potentials, but there exist bound states within the confining 
potentials. Meanwhile, it is interesting to see that all states with / > 0 have their 
own pseudospin partners in the confining potentials, as shown in Fig. 14, but this 
is not the case for the non-confining potentials appearing in the isolated atomic 
nuclei, as shown in Fig. 4. We consider this as a puzzle of intruder states and will 


discuss in more detail in Section 4.2. 


3.3. From local potentials to non-local potentials 


So far, most of studies on the pseudospin symmetry mainly focus on the 
single-particle Hamiltonian with only local potentials. Ldopez-Quelle et al. 
[171] performed one of the first investigations of the pseudospin symmetry 
with non-local potentials in the framework of the relativistic Hartree-Fock 
theory [181, 182]. Recently, the relativistic Hartree-Fock theory achieved 
lots of success in describing nuclear ground-state and excited-state properties, 
by introducing the density-dependent meson-nucleon couplings by Long and 
coauthors [53, 183, 306]. Along this direction, more detailed investigations of 
the pseudospin symmetry in the relativistic Hartree-Fock theory were performed 
in Refs. [45, 175], as well as the spin symmetry in the anti-nucleon spectra in 


Ref. [198]. 
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In this Section, we will briefly introduce the relativistic Hartree-Fock theory, 
and then mainly focus on the special features of pseudospin symmetry due to 
the non-local potentials. The spin symmetry in the anti-nucleon spectra will be 


discussed in Section 3.6. 


3.3.1. Relativistic Hartree-Fock theory 

Similar to the RMF theory, the starting point of the RHF theory [181, 182] 
is a Lagrangian density (52), in which the nucleus is described as a system 
of Dirac nucleons that interact with each other via the exchange of mesons and 
photons. With the general Legendre transformation (55) as well as the mean-field 
and no-sea approximations (56), one has the energy density functional for the 


whole system, 
Errr = (®o| 2 |®o) . (89) 


However, different from the RMF theory, here both the direct and exchange 
contributions of the meson and Coulomb fields, the so-called Hartree and Fock 
terms, are included. 

In such a way, the effective nucleon-nucleon tensor interactions can be 
naturally taken into account, which are practically absent at the Hartree level. 


The corresponding parts in the Lagrangian density read [53, 306] 
ZT = Seg WO FR 
pe ` Vind vPu y >, 


meat it, (90) 


NI — 


s 1 
EA = ET . tw + 5 Out OR - 
My 


corresponding to the p-N tensor and z-N pseudovector couplings. These tensor 
interactions are crucial for improving the descriptions of the nuclear shell 


structures and their evolutions [50, 51, 53, 54, 306]. 
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In recent studies, it has been also shown that the meson exchange terms 
play very important roles in the nucleon effective mass splitting [183], sym- 
metry energies [307], halo structure [45], deformation [308], and spin-isospin 
resonances [309, 310]. Meanwhile, the Coulomb exchange term [311] is crucial 
for understanding the isospin symmetry-breaking corrections to the superallowed 
B decays [312] in order to test the unitarity of the Cabibbo-Kobayashi-Maskawa 
matrix [313]. 

In the RHF theory, the eigenfunction equations for nucleons involve non- 
local potentials due to the finite-range meson and photon exchanges. For the 
spherical case, the corresponding radial Dirac equations are expressed as the 


coupled integro-differential equations [182], 


M+) -24i G(r) Y(r) G(r) 
d ‘ dr r + =E 5 (91) 
ae AP(r) || FO) X(r) F(r) 


where È? and A? contain the contributions from the Hartree terms as shown in 
Eq. (10), while X and Y functions represent the results of the non-local Fock 


potentials acting on the wave functions F and G, respectively. 


3.3.2. PSS in non-local potentials 

As shown in Refs. [45, 171, 175], in many cases the single-particle spectra 
obtained in the self-consistent RHF theory show a good PSS. It was also found 
that the stability of neutron halo structures, e.g., in the drip-line Ce isotopes, is 
closely related to the PSS conservation of the single-proton spectrum. In this PSS 
conservation the p-N tensor coupling in Eq. (90) plays an essential role via the 
Fock terms [45]. However, in these cases, it is more sophisticated to trace the 


origin of the PSS and the exact PSS limit is no longer as simple as d£?” (r)/dr = 0 
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(17) because of the non-local potentials [171]. 
One possible way for tracing the origin of the PSS with non-local potentials 


was pointed out by Long et al. [175] by introducing the localized equivalent 


potentials, 
_ _ GMX) FO)X(7) z 
XO = ap P O t po rO EXON + XKOFO, 
GOY) F(r)Y(r) _ 
YO = FO * Bigg FG = PaO) + HFC, 
(92) 
so that the radial Dirac equations (91) can be rewritten as 
d K 
E == r) F(r)- [X(r) - E] G(r) = 0, 
r r 
£ + ; + xan] G(r) + [A(r)- E] Fr) = 0, (93) 


where E = e— M, È =X? + Yg, and A = A? + Xp — 2M. 


The corresponding Schrédinger-like equation for the lower component F(r) 


reads 
qa ig ec + Veso + V2) F = -(A -EÈ -—E\)F, (94) 
with 
1 dz 
-(X,-Y 2 
V = (Xc — Yr) Soar’ 
V =Y i XY. Yr + YX’ - E)+ XF È -E 
2 = Yro p; olac r G( )+XrÈ- E), 
K 1 d& 
Vpso zra, TT) ; 
k(l -K 
Va =P, (95) 
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where Vpcg and Vpso correspond to the PCB and PSO potential, respectively. 
While the potential V entirely originates from the Fock contributions, the 
Hartree and Fock contributions to V, and the PSO potential can be approximately 


separated as, 


1 d=? 1 dY, 
Do = 2% E y yp Ke 
Y= ap Ya en) = ar” 
5 K 1 g x x| i do 
= -— —— pe ey 
Vrso r&—-E dr’ PSO" r|Z-E dr (Xe + Erd) (26 


denoted with superscripts “D” and “E”, respectively. 

Finally, to have a better understanding of the PSO splitting, especially the 
effects of non-local Fock terms, it will be more transparent to rewrite Eq. (94) as 
[175] 


1 æ 


1 D E D 
oro o ee eee, On) 


where the operators V? and VF are respectively 
d 
VeVi PV and VF =VF— ont V. (98) 
r r 


From Eq. (97), one can estimate the contributions of the potentials Vpcg, V”, and 


VE to the single-particle energy E. For example, the PCB contribution can be 


evaluated by 
> Vee 19 f S 2 
F F*dr. 
> IE ar| f dr (99) 


Although there exists a singularity in 1/(A? — E), it has been proven that the 
principal values of these integrals are finite due to the nodal structure of F(r) 
(72, 167]. 

In Table 2, the results calculated by the self-consistent RHF theory with the 


effective interaction PKO1 [183] are shown for the neutron 1 and 2 pseudospin 
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Table 2: Single-particle energies E in '**Sn and the contributions from each term in the left-hand 
side of Eq. (97) given by the RHF theory with the effective interaction PKO1 [183], in comparison 
with those by the RMF theory with the effective interaction PKDD [290]. All units are in MeV. 
The data are taken from Ref. [175]. 

State E F” EP Vecer qyP yE 


RHF v2sı2 -31.41 18.11 -75.35 9.30 -2.99 19.51 
PKOI vld -34.90 14.87 -79.01 9.54 0.44 19.26 
v3sı2 —8.33 34.25 -72.00 11.11 0.09 18.22 
v2d3,2 -8.66 31.93 -73.96 11.32 3.89 18.17 
RMF v2sızn —34.81 21.86 -64.65 11.04 -3.07 — 
PKDD vld; -38.87 18.17 -68.08 11.41 —0.37 — 
v3sı2  —8.15 40.13 -61.97 13.02 0.67 — 
V2d3,. -8.44 37.65 -63.75 13.36 4.30 — 
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Figure 16: (Color online) Radial wave functions G and F (left panels), and the non-local terms 
X and Y (right panels) for the v1 and v2 doublets in '*Sn given by the RHF theory with the 
effective interaction PKO1 [183]. Taken from Ref. [175]. 
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doublets in '°’Sn. For comparison, the corresponding results calculated by the 
RMF theory with the effective interaction PKDD [290] are shown in the lower 
part. In addition, their radial wave functions G(r) and F(r) as well as the non- 
local terms X(r) and Y(r) given by the RHF theory are shown in Fig. 16. 

First of all, the RHF and RMF theories share the same properties on the F”, 
EP, Vpcg, and VP terms. The differences between the pseudospin partners in 
the PCB are very small, whereas the differences in the F”, £P, and VP terms are 
substantial individually. However, these three terms cancel largely one another 


and the PSS is conserved very well. 


T $ T 
8 F DDRHF i v2s,. T gp vid, | 
aL ii -31.41 MeV | -34.90 MeV 
| a n 
0 7 Y EO $ 
4b Y — vio FI (VE) + i 7 
go ViFFV(V-E) | i 4 
f J p Į | | al | 
T T aa Vr ne A 
v3s i v 2d. 


112 " 3/2 
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3 t h 
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Figure 17: (Color online) Functions Vf..F*/(A? — E) and VE F F’/(A? — E) for the v1p and v2p 
doublets in !?Sn given by the Fock terms of the RHF theory with PKO1. The singular points at 
r ~ 6 fm for v1p and at r ~ 7.5 fm for v2p are due to the denominator (2 — E) in the PSO potential, 


while the other local peaks are due to the nodes of the upper component G. Taken from Ref. [175]. 


For the special features in the RHF theory, it is found that the gross differences 
in the Fock terms V® are almost negligible. The main reason is that there exist 


significant cancellations between Vf. and VF in Eq. (98), especially in the 
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inner part of the nucleus. The functions V5,,F?/(A? — E) and ViFF'/(A? — E) 
for the neutron 1p and 2p doublets are illustrated in Fig. 17. In other words, 
although the Fock terms bring substantial contributions to the PSO potential, these 
contributions are canceled by the other exchange term V”, which stems mainly 
from the non-locality of the state-dependent exchange potentials. Furthermore, 
it was pointed out that these cancellations are not accidental but because of the 
similar radial dependence between the Fock-related terms (X, Y) and the wave 
functions (F, G) [175], as shown in Fig. 16. 

In summary, the Fock terms in the relativistic Hartree-Fock theory bring 
significant contributions to the pseudospin-orbit potential and make it comparable 
to the pseudo-centrifugal barrier. However, these Fock terms in the pseudospin- 
orbit potential are counteracted by other exchange terms due to the non-locality 
of the exchange potentials. The physical mechanism of these cancellations was 
discussed in relation to the similarity between the exchange potentials and the 
Dirac wave functions. Therefore, although the non-local potentials make the 
analysis more much complicated, they do not substantially violate the general 


features of pseudospin symmetry described by only using the local potentials. 


3.4. From central potentials to tensor potentials 


The tensor component is a non-central contribution of the nucleon-nucleon 
interaction, which has been extensively discussed during the past decades. Ab ini- 
tio calculations based on realistic nucleon-nucleon interactions have demonstrated 
the important role played by the bare tensor force in the description of the binding 
energy of nuclei [314, 315]. Its role has also been investigated for nuclear matter 
with realistic potentials [316]. 


With the rich spin-isospin structure of finite nuclei, it is reasonable to expect 
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that some nuclear observables are sensitive to the tensor force. For instance, the 
shell evolution and the modification of magic numbers far from stability are often 
discussed and interpreted in terms of tensor effects [10]. In the framework of 
the shell model, the tensor contribution in shell evolution has been explored and 
its effects have been underlined by Otsuka et al. [317-319]. In the mean-field 
models, the tensor component has been also extensively discussed during the 
past years, including the non-relativistic Skyrme [320-324] and Gogny [325, 326] 
theories. The readers are referred to Ref. [327] for a recent review. 

As the covariant symmetry is conserved in the covariant density functional 
theory, the tensor interactions in this scheme usually indicate the interactions 
with Lorentz-type tensor couplings, i.e., those involving the Dirac matrix o*” 
or ysy”. The Lorentz-type tensor interactions lead to both central and non- 
central potentials, and the later ones correspond to the non-relativistic tensor 
components embedded in the covariant couplings. The Lorentz-type tensor 
effects on the shell evolution have been discussed in the relativistic Hartree-Fock 
theory [53, 183, 306]. Furthermore, the tensor effects in the relativistic and non- 
relativistic theories were systematically compared in Refs. [50, 51]. 

In this Section, we will mainly focus on three different examples for illus- 
trating the tensor effects on the pseudospin symmetry. The first example is 
the relativistic harmonic oscillator with a linear tensor potential, which admits 
analytical solutions [89]. Self-consistently, the w-N tensor coupling is one of 
the widely used ways to include the tensor potential for the covariant density 
functional theory in the Hartree level [328]. The corresponding effects on 
the pseudospin symmetry were investigated in Ref. [173]. Similarly, the w-A 


tensor effect on the spin symmetry in the single-A spectra will be discussed in 


67 


Section 3.7. Finally, we will discuss the tensor effects in the relativistic Hartree- 
Fock theory [53], which is considered as a more complete way to include the 
effective nucleon-nucleon tensor interactions in the scheme of covariant density 


functional theory. 


3.4.1. Linear tensor potential 
When the Lorentz-type tensor couplings of the vector mesons to nucleons are 
taken into account in the CDFT [328], the single-particle Dirac Hamiltonian H in 


Eq. (2) becomes 
H=a-p+f[M + S(r)] + V(r) -iBa- T(r), (100) 


where T(r) is the tensor potential. If the tensor potential is simply a function of 


the radial coordinate r, this Hamiltonian can be simplified as 
H=a-p+f[M+S(r)] + VY) - Ba -fT(r). (101) 


When the spherical symmetry is adopted, the corresponding radial Dirac equation 


shown in Eq. (10) becomes 


d 
Man -=+--T0) |( eo G(r) 
P a i =€ À (102) 
a -M+A(r) F(r) F(r) 


As mentioned in Section 3.2, there exist bound states in the Fermi sea at the 
PSS limit (17), dX(r)/dr = 0, as long as the potential A(r) is confining. An 
example shown in Figs. 14 and 15 corresponds to the case of relativistic harmonic 
oscillator (20) with X(r) = 0 and A(r) = May,r7/2. On top of this PSS limit, 
Lisboa et al. [89] investigated the effects of the tensor potential T(r), which was 


assumed as a linear function of r, 
T(r) = Mayr. (103) 
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The corresponding single-particle energies can be obtained analytically as [89] 


TE 
€ — M? - (2k + 1)Mo = (on +1+ 5) J2M(e- Mw? + 423. (104) 


Here ñ is the number of the internal nodes of F (r), i.e., np. 


16 


Figure 18: Single-particle energies as a function of w2 with X(r) = 0, wı = 2, and M = 10. Here 


ñ = nr is the number of the internal nodes of F(r). Taken from Ref. [89]. 


In Fig. 18, the single-particle energies of 15 state and 1) pseudospin doublets 
are shown as a function of the tensor potential strength w2. It is shown that in 
general the single-particle energies increase with w2, as w; appears in the right 
hand side of Eq. (104). More importantly, the pseudospin-up states with j> = 
Ī + 1/2 increase faster than the pseudospin-down states je = /— 1/2. This is 
because x > 0 (x < 0) for the j> (j<) orbitals appearing in the left hand side of 
Eq. (104) . 

In this simple example, the SO splittings are not taken into account. Therefore, 
the j> orbitals go higher in energy than the j< orbitals once the tensor potential is 


included. 
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3.4.2. Tensor potential of w meson 

For the CDFT in the Hartree level, i.e., the RMF theory, the w-N tensor 
coupling is one of the widely used ways to include the tensor potential self- 
consistently [328, 329], 


fo - 
L =- Quip, 105 
where fo is the w-N tensor coupling strength and Q, is the field tensor defined in 


Eq. (53). In the spherical systems, the tensor potential thus obtained is 


T(r)= Le) x — =V (r), (106) 


which is approximately proportional to the derivative of the vector potential V(r) 
if the contribution from the p meson is negligible. With this tensor potential the 


Schrödinger-like equation for the lower component F (r) reads 


2 — 1 dM_ 
{4 K(k — 1) d Jo -S4r 
2 


T , ‘| 
+u rT =T F(r) =-M, MFG). 
(107) 


In Ref. [173], this equation is further expressed as 


2 apaatia AT __ >% dF) |_ xY \L-Ss 
pro+(r -T -2 27) Fn = + 4T +2) —F(r) 
=(E-X)(E+2M-A)F(r). (108) 


By dividing it by 2M* = E + 2M — A, the energy decomposition of the above 


equation reads [173] 


2 
(=) + (Vr) + (Ver) + Voanin) + (Veso) + (2) = E, (109) 
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(V) = EE (110) 
of F*Fdr ” 
and 
1 T 
Vp = zp (T-T -2-), 
1 YT 
Vy, 
2T OM* E->’ 
1 Y d 
V Darwin =.  OM* E— dr > 
1 Y \L-S 
Ving 2 Sf ee 111 
PSO arl z) T (111) 


which correspond to the tensor, Darwin, and PSO contributions, respectively. For 
the terms with E -È in the denominator, the integral is taken in the principal value 
sense. 

As an example, the single-particle energies of 1 and 1h pseudospin doublets 
in 7°8Pb and the contributions from each term in the left-hand side of Eq. (109) 
are calculated with Woods-Saxon potentials. The results without and with tensor 
potential in Eq. (106) are compared in Table 3. 

First of all, because the tensor interaction is a higher order interaction in the 
Lagrangian scaled by 1/M, the contributions of the Vr and Vyr terms in energy 
are small compared to the others. However, because the SO or PSO interaction is 
a term of the same order, the tensor effects on the SO and PSO splittings can be 
significant. 

Focusing on the pseudospin-orbit splittings AEpso, for the deeply bound 
pseudospin doublets, the contribution of the PSO potential Vpso almost does 
not change with fù, whereas the contribution of terms such as G /2M*) and 


(È), which do not depend explicitly on the tensor potential T, is larger than 
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Table 3: Single-particle energies E of 1p and 1h pseudospin doublets in 7°8Pb and the contributions 


from each term in the left-hand side of Eq. (109) calculated by the Woods-Saxon potentials without 


and with tensor potential. The Woods-Saxon parameters are R = 7 fm, Ag = 650 MeV, Xo = 
—66 MeV, and a = 0.6 fm. All units are in MeV. The data are taken from Ref. [173]. 


fal Be 
0.0 


1.3 


0.0 


1.3 


State 
2581/2 
1d3) 
AE 

281/2 
1d3/2 
AE 

289/2 
lini 
AE 

2g9/2 
lini 


AE 


(p?/2m") 
24.4396 
21.1032 

3.3364 
23.9037 
20.7075 

3.1962 
55.7666 
51.2033 

4.5633 
56.0706 
51.6221 

4.4485 


(Vr) 
0.0000 
0.0000 
0.0000 
0.2351 
0.1459 
0.0892 
0.0000 
0.0000 
0.0000 
0.5622 
0.8997 

—0.3375 


(Vyr) 
0.0000 
0.0000 
0.0000 

—0.0943 
—0.0712 
—0.0232 
0.0000 
0.0000 
0.0000 
0.2855 
0.1529 
0.1326 
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(Vparwin) 
-3.9527 
—0.8106 
—3.1421 
—3.6870 
—0.4170 
—3.2700 
2.8295 
3.1788 
—0.3493 
2.6487 
2.9268 
—0.2781 


(Vpso) 
—0.5852 
0.0966 
—0.6818 
—0.5632 
0.1559 
—0.7191 
—6.4555 
3.3530 
—9.8085 
—7.2110 
4.1391 
—11.3501 


(2) E 
—61.4644 -—41.5627 
—64.4159 —44.0266 

2.9515 2.4639 
—61.9114 —42.1170 
—64.6181 -—44.1170 

2.7067 2.0000 
—53.3816 —1.2410 
—61.0308 -3.2958 

7.6492 2.0548 
-54.4705 —2.1140 
—61.2232 —1.4826 

6.7527 —0.6314 


that of Vpso. This means that the main contribution to the change of AEpgso is 
mainly due to the change of the wave function induced by 7. Furthermore, the 
energy splitting results from a strong cancellation of the A (È), A(Vparwin), and 
A (p?/2M"*) contributions [173]. 

For the loosely bound pseudospin doublets, much of the previous analysis still 
holds, except for the fact that the PSO potential Vpso is much stronger and also 
changes sensibly as f changes, being responsible for most of the pseudospin 
splitting. For the 1h doublets, the AEpso even becomes negative when fo/8w 
greater than 0.95. It is found that the tensor part of Vpso in Eq. (111) is the 
dominant ingredient to understand this evolution [173]. Because the tensor 
potential here is proportional to the derivative of the vector potential, such tensor 
effects are much more profound for the surface states than the deeply bound ones. 

In short, the pseudospin-orbit splittings AEpso = E;. — E;, decrease, i.e., the 
PSS is conserved better, when the w-N tensor interaction is included. Due to the 
surface-peaked shape of the tensor potential, the tensor effects are much more 


profound for the loosely bound states than the deeply bound states. 


3.4.3. Tensor potential of p and x mesons 

As shown in Section 3.3, a more complete way to include the effective 
nucleon-nucleon tensor interactions in the CDFT scheme is the RHF theory [53]. 
The p-N tensor and z-N pseudovector couplings shown in Eq. (90) are crucial 
for improving the descriptions of the nuclear shell structures and their evolutions 
[50, 51, 53, 54, 306]. In particular, these tensor interactions play important roles 
via the exchange terms. 

In Figs. 19 and 20, the single-particle energies in '*Sn and SPb calculated 


by the self-consistent RHF theory with the effective interactions PKA1 [53] and 
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Figure 19: Single-particle energies in '3?Sn calculated by the RHF theory with the effective 
interactions PKA1 [53] and PKO1 [183], and by the RMF theory with the effective interactions 
PK1 [290] and DD-ME2 [293], compared with the experimental data [330]. Taken from Ref. [53]. 
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Figure 20: Same as Fig. 19, but for the single-particle energies in SPb. Taken from Ref. [53]. 
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PKO1 [183] effective interactions are compared with the corresponding results 
calculated by the RMF theory with the effective interactions PK1 [290] and DD- 
ME2 [293]. Note that both p-N and z-N tensor couplings are included in PKA1, 
and only z-N tensor coupling is included in PKO1, while the tensor effects are 
practically absent in the RMF models. 

Compared with the experimental data [330], the RMF results always show 
too large PSO splittings between the 2ds; and 17/2 states in '**Sn as well as the 
splittings between the 2f;;2 and 1/o/2 states in SPb. This is the case not only 
for the neutron side but also for the proton side. As a result, the overestimated 
sub-shell structures, e.g., N or Z = 58 and N or Z = 92, appear commonly in the 
RMF calculations without tensor interactions [331]. 

The artificial sub-shell closure is depressed to some extent by using the RHF 
theory with PKO1 effective interaction. It is found the PSO splittings of the 1f 
and 1g doublets decrease compared to the RMF results. However, because of the 
strong density dependence of z-N coupling strength f,(o) in PKO1, the net tensor 
effects are mild, and the calculated PSO splittings are still substantially larger than 
the experimental ones. 

The PSO splittings of the 1f and 1% doublets decrease further when both 
the z-N and p-N tensor couplings are included, which is the case of PKA1, and 
eventually the artificial sub-shells disappear. Indeed, these tensor effects can be 
seen in all the pseudospin doublets shown in the figures by comparing the PKA1 
results to the others. This is also crucial for the descriptions of the ordering of the 
single-particle levels, e.g., the neutron states 289/2 and 1i;4/2 in *°°Pb. 

In summary, three different models of tensor forces were discussed in this 


Section. From a phenomenological linear tensor potential to more microscopic 
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meson-nucleon tensor couplings, these models show consistent effects on the 
pseudospin symmetry, i.e., AEpso = E;. — Ej, decrease or become even negative 
when the strengths of tensor interactions increase. This mechanism plays an 
important role in the nuclear shell structure and its evolutions. For example, the 
overestimated sub-shell closure at N or Z = 58 and N or Z = 92 can be better 


reproduced. 


3.5. From bound states to resonant states 


Weakly bound or unbound nuclei with unusual N/Z ratios are open quantum 
many-body systems in which the continuum plays an important role [332]. 
In these nuclei, the neutron (or proton) Fermi surface is close to the particle 
threshold, thus the contribution of the continuum is crucial [60, 61, 190, 192-194, 
203-218]. Many approaches developed for resonances [219], e.g., the analytical 
continuation in coupling constant method [220-225], the real stabilization method 
[226-231], the complex scaling method [232-235], the coupled-channel method 
[236-238], and some others [239, 240], have been used to study nuclear single- 
particle resonant states. Based on some of these methods, the pseudospin 
symmetry [241-245] and the spin symmetry [246] in single-particle resonant 
states have been investigated. The pseudospin symmetry and/or spin symmetry 
in nucleon-nucleus and nucleon-nucleon scatterings have been also investigated 
[41-44]. 

In Ref. [247], by examining the asymptotic behaviors of the Dirac wave 
functions, a rigorous verification of the pseudospin symmetry in single-particle 
resonant states was given by Lu, Zhao, and Zhou. It was shown that the 
pseudospin symmetry in single-particle resonant states in nuclei is exactly 


conserved under the same condition for the pseudospin symmetry in bound states, 
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i.e., X(r) = 0 or dX(r)/dr = 0 given in Eq. (17). To understand more deeply the 
pseudospin symmetry in single-particle resonant states, in Ref. [249], the exact 
conservation and breaking mechanism of the pseudospin symmetry in single- 
particle resonant states in square well potentials were extensively studied. A 
threshold effect in the energy splitting and an anomaly in the width splitting of 
pseudospin partners were found when the depth of the potential varies from zero 
to a finite value. 

In this Section, we will introduce the pseudospin symmetry in single particle 
resonant states. For the pseudospin symmetry in nucleon-nucleus and nucleon- 


nucleon scatterings [41—44], see Ref. [67]. 


3.5.1. PSS in single-particle resonant states with ACCC 
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Figure 21: Left panel: Energies and widths for the neutron resonant states 3p3/2, 3p1/2, 2f7/2, 
2f5/2, [ho/2, and 1113/2 in the potential of 2277 Solid circles represent the results of the ACCC 
method, while open circles denote the results of the scattering method. Right panel: Real parts 
of the upper and lower components of radial wave functions for the neutron resonant state 2 f7/2. 
Solid and dashed curves represent the results of the ACCC and the scattering method, respectively. 
Taken from Ref. [221]. 


In Ref. [220], the ACCC method was combined with the covariant density 
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functional theory to study the unbound states in spherical nuclei. The ACCC 
method [333, 334] is based on the following idea: An unbound state can be 
considered as the continuation of a bound one when the strength of the attractive 
potential decreases. That is, an unbound state will become bound if the coupling 
strength of attractive potential is increased and becomes strong enough. As a 
function of the strength, the energy of a bound state is analytically continued to 
the complex plane in order to get the width and energy of a resonant state. 

In Ref. [221], the analyticity of the eigenvalues and eigenfunctions of the Dirac 
equation with respect to the coupling constant was examined. In the left panel of 
Fig. 21, the energies and widths for the neutron resonant states 33/2, 31/2, 2f7/2, 
2 fs/2, Lho, and 1i}3/2 in the potential of "Zr are shown in a planar E-I plot. The 
results from the ACCC method agree satisfactorily with those from the scattering 
calculation. The wave functions for the neutron resonant state 2 f7/. are exhibited 
in the right panel of Fig. 21. The ACCC and the scattering method give not only 
nearly the same resonance parameters, energies and widths, but also very similar 
wave functions. 

In Ref. [241], the PSS for the resonant states in 7°°Pb was investigated by 
solving the Dirac equation with the Woods-Saxon scalar and vector potentials, 
and the ACCC method was used to determine the resonance parameters. The 
PSS breaking is shown in correlation with the nuclear mean field shaped by the 
central depth Xo, the radius R, and the diffuseness a. The energy-level crossings 
appear in several pseudospin partners of resonant states. The width is found 
to be different for the pseudospin doublets even when their energies are fully 
degenerate. This was attributed to the difference in the centrifugal barrier for 


the pseudospin partners [241]. 
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The ACCC method, combined with the CDFT, was used to study the N- 
dependence of PSS in nuclear resonant states in Ref. [242]. The energies 
and widths of single-particle resonant states in Sn isotopes were investigated 
systematically. An N-dependence of PSS is clearly shown in the resonant states 
and is consistent with that observed in the bound states. The splittings of energies 
and widths between pseudospin doublets are correlated with the quantum numbers 
of the single-particle states, as well as the nuclear mass number. Although the N- 
dependence of the pseudospin splitting in energy is rather complicated, the width 


splitting decreases almost monotonically with increasing mass number. 
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Figure 22: The real part of the upper component G(r) (the left panel) and the lower component 
F(r) (the right panel) for the proton resonant state 19/2 (dashed line) compared with the proton 


resonant state 2 f7/2 (solid line) in 120Sn. Taken from Ref. [245]. 


In Refs. [244, 245], the ACCC combined with CDFT was also used to discuss 
the PSS in the single-proton resonant states. Not only the resonant energies 
and widths, but also the wave functions for pseudospin partners in '*°Sn were 
examined. In Fig. 22 are shown the real parts of the upper and lower components 
of the radial wave functions for the 71% pseudospin doublets in !7°Sn. It is seen 


that the lower components agree very well in the region where nuclear potential 
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dominates, except for some disagreement on the surface, with the number of radial 


nodes being the same. 


3.5.2. PSS in single-particle resonant states with CSM 
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Figure 23: (Color online) Pseudospin splittings in energy and width as a function of Morse 
potential parameters, (a) Vo, (b) ro, and (c) œ, for single-particle resonant states. Taken from 


Ref. [243]. 


In Ref. [243], the CSM was applied to study the PSS of resonant states of a 
Dirac particle in a Morse potential, 


V(r) = Voe 7 (2 =e eae), (112) 


By comparing the energies and widths of the pseudospin doublets, the PSS was 


examined and the relationship between the PSS and the parameters of the Morse 
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potential was studied. In Fig. 23, pseudospin splittings in energy and width for 
single-particle resonant states are shown as a function of the Morse potential 
parameters Vo, ro, and œ, respectively. With rọ and a fixed, the energy splitting of 
a pair of pseudospin partners increases with increasing Vo. However, the opposite 
is true for the width splitting. This is because, with increasing Vo, the barrier 
becomes higher and the potential well becomes deeper. The dependence of the 
energy and width splittings on rọ or œ were also made, and the details can be 


found in Ref. [243]. 


3.5.3. PSS in single-particle resonant states with Jost functions 

Although the PSS in resonant states have been studied numerically, it is still 
desirable to have a mathematical verification of the PSS in the symmetry limit. 
In particular, there remains an open question about the width splitting of the 
pseudospin doublets. Recently, a rigorous justification of the PSS in single- 
particle resonant states has been given by examining the asymptotic behavior of 
the nucleon Dirac wave functions [247]. 

The first-order coupled equation (10) can be rewritten as two decoupled 
second-order differential ones, Eqs. (12) and (14). In Refs. [275-278], it has 
been shown that each of these two Schrédinger-like equations, together with its 
charge conjugated one, are fully equivalent to Eq. (10). The one for the small 
component is given in Eq. (14) to which the PSS is directly connected. Note that 
for bound states, there is always a singularity in 1/M_(r) in Eq. (14), whereas such 
a singularity does not exist for resonant states discussed here. 

For the continuum in the Fermi sea, i.e., € > M, there exist two independent 
solutions for Eq. (14). The physically acceptable solution is the one that vanishes 


at the origin. As usual we define the regular solution F(r) as the one that behaves 
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like j;(pr) function as r > 0 [335], 
lim F(r)/ji(pr) = 1, p= Ve- M°. (113) 


At large r the potentials for neutrons vanish and the radial wave functions os- 
cillate. Equation (14) becomes a Ricatti-Bessel equation with angular momentum 


Ī, and the solution can be written as a combination of the Ricatti-Hankel functions, 
i - * 
FO) = 5 (Te Phy (0) - Ihe], > 0, (114) 


where JE (p) is the Jost function for the small component and he pr) the Ricatti- 
Hankel functions. 
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Figure 24: (Color online) Schematic picture of the zeros of the Jost function J? on the complex 
momentum plane. A cut is made on the imaginary axis, from p = iM to infinity and back to 


p = —iM. Taken from Ref. [247]. 


In nuclei the vector V(r) and scalar S(r) potentials share some general 
properties, e.g., they are analytical functions of r, vanish as r — oo, and have 
no singularities. Under such conditions, the Jost function is an analytic function 


of p, and can be analytically continued to a large area in the complex p plane. Here 
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the structure of the p Riemann surface on which the Jost functions are defined is 
more complex than the non-relativistic case. For example, the square root in the 
relativistic energy-momentum relation €? = p? + M? creates branching points at 
p = +iM, thus the corresponding Riemann surface has at least two sheets. 

In Fig. 24, the zeros of the Jost function .J/"(p) on the complex momentum 
plane are schematically shown. Only the first sheet with Re(e) > 0 is presented in 
Fig. 24 which contains positive-energy bound states and resonant states, while the 
other sheet with Re(e) < 0 can be used to investigate negative-energy ones. These 
two sheets are connected by a cut on the imaginary axis, from p = iM to infinity 
and back to p = —iM. Restricted to the first sheet and relatively small |p|, J" (p) 
is a single-valued analytical function of p. The zeros of J/(p) are denoted by 
the full circles for bound states, open circles for resonant states, and crosses for 
other zeros, respectively. The zeros on the positive imaginary axis of the p plane 
represent bound states of the original eigenvalue problem, while the zeros on the 
lower p plane and near the real axis correspond to resonant states. The resonance 
energy Ees and width Ies are determined by the relation E = Ees — il res/2 = 
Vp? + M? — M. By examining the zeros of the Jost function, one can study the 
bound and resonant states on the same footing, and many known properties of 
bound states can be generalized to resonances straightforwardly. 

At the PSS limit, Eq. (14) is reduced to Eq. (19). For the bound states, it 
is an eigenfunction equation that determines the eigenenergies e. While for the 
continuum e€ can be any value larger than or equal to M and one should focus on 
the wave functions and their asymptotic behavior. For pseudospin doublets with 
different quantum numbers K, and x, with Kp, = —k, + 1, the small components 


satisfy the same equation because they have the same pseudo-orbital angular 
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momentum / [65]. In particular, for continuum states, F«,(€,r) = F,,(e,r) for 
any energy e. Because the definition of the Jost function J/(p) only depends 
on the asymptotic behavior of the small component, on the positive real axis, 
STi (p) = Fi (p). This equivalence can be generalized into the complex p plane 
due to the uniqueness of the analytical continuation. Thus the zeros are the 
same for J, Fa (p) and J, A (p): If there exists a resonant state with energy E,es and 
width Tes and the quantum number x,, there must be another one with the same 
energy and width and quantum number «xp. That is to say, the PSS in single- 
particle resonant states in nuclei is exactly conserved when the attractive scalar 
and repulsive vector potentials have the same magnitude but opposite sign. If 
one focuses on the zeros of the Jost functions of pseudospin doublets on the 
positive imaginary axis of the p plane, the well-known PSS for bound states can 
be investigated similarly. 

In scattering theories, one can also determine resonance parameters from the 
change of cross section or phase shift, which give more insights into the resonant 
phenomena. Using the asymptotic behavior of the Ricatti-Bessel functions, one 


obtains from Eqs. (19) and (114), 
| E 
Fr) œ sin pr—- = +6, (p) , r> œ, (115) 
where the phase shift 6’ (p) is related to the Jost function through 
TEP) = TEPO. (116) 


Whenever 6/(p) = (n + 1/2)z, there is a resonant state and its width is determined 
by the tangent of the phase shift function 6‘ (p). In the PSS limit, the coincidence 
between F,,,(r) and F,,(r) means that of (p) = 6f (p) for any value of p. Therefore, 


resonance parameters of pseudospin doublets are the same. 
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As it has been noted in Ref. [247], it is straightforward to extend the study of 
PSS in resonant states in the Fermi sea to that in the negative-energy states in the 
Dirac sea or SS in anti-particle continuum spectra. More investigations along this 


line should be made. 


3.5.4. PSS in single-particle resonant states in square-well potentials 

To extract the energies and widths of resonant states in realistic potentials is 
relatively complex. In particular, to study the PSS and examine the origin and 
the symmetry breaking mechanism, it is better to start from analytically solvable 
models. In Ref. [247], the square-well potentials were taken as examples to 
illustrate the conservation and breaking of the PSS in the single-particle resonant 
states. Although the diffuseness of realistic potentials cannot be included, it is still 
a good starting point to study general properties of the PSS for the resonant as well 
as bound states by using the square-well potentials, because the PSS-breaking 
term in the Jost function is separated from the PSS-conserving term. 


Spherical square-well potentials for Ł(r) and A(r) read 


C, r<R, D, r<R, 
Lr) = A(r) = (117) 
0, r>R, 0, r>R, 


where C and D are depths and R is the width. The Jost function .J/(p) is derived 
as [247] 


p' 


. + ” + C ” K , + 
TEO) =~ SP {J RR) phi (DR) — KJR (DR) - — r | RipRR) - EIER | DRY} 


M 
(118) 
with k = y(e-— C — M) (e — D + M). The PSS in both bound states and resonant 


states can be explained explicitly. If C = 0, the second term in JE (p) vanishes 
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and the first term only depends on the pseudo-orbital angular momentum /, then 
the Jost functions with different x but the same / are identical. The energies and 


widths of resonant pseudospin partners are exactly the same. 
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Figure 25: (Color online) The zeros of the Jost function JF on the complex energy plane in 
square-well potentials (117) with C = 0 (solid symbols) and C = —66 MeV (half-filled symbols) 
for pseudospin partners 11/2 (diamond) and j13/2 (square). The results with Woods-Saxon-like 


scalar and vector potentials are also shown as open symbols. Taken from Ref. [247]. 


For the large component G(r), there is a similar expression for the asymptotic 


behavior, 

G(r) = = [Ioh (wr) -IT orken, r>w. (119) 
At the origin, 

limG()/j(pr)=1, p=Ve-M?. (120) 


The J °(p) is derived as 


: D 
TEP) = -Era (JR phi PR) — KI (ERYH (PR) - —=— |KER) + = jkR)| HDR}. 


(121) 
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It looks similar to that of the small component J, r (p), with the exception that 
the potential parameter C is substituted by D and the pseudo-orbital angular 
momentum / is substituted by L. In the case of D — 0, this form of Jost function 
can be used to investigate the spin symmetry of single-particle levels. 

Although the solutions of J (p) = 0 cannot be derived analytically, the secant 
method can be used for searching the roots, because the Jost function is analytical 
near its zeros. In Fig. 25 are shown solutions in the complex energy plane for 
PSS doublets with 7 = 6, i.e., hi2 with ka = —6 and ji3/. with «xp = 7, for 
square-well potentials with D = 650 MeV and R = 7 fm. In the PSS limit, i.e., 
C = 0, all the roots locate in the lower half plane and there are no bound states. 
In Fig. 25 three pairs of pseudospin resonant partners are shown by full diamonds 
and squares. The conservation of the PSS for single-particle resonant states is 
clearly seen. When C = —66 MeV, there is one bound state only for /j,/2. Three 
pairs of pseudospin partners of resonant states are shown by half-filled diamonds 
and squares. One finds the breaking of the PSS both in the bound states and in the 
resonant states. For pseudospin doublets with other values of /, similar behaviors 
are observed concerning the exact conservation and the breaking of the PSS. The 
resonances in Woods-Saxon potentials, W(r) = Wo/(1 + exp[(r — R)/a]) (W = V 
or S ) were also studied and the potential parameters are the following: the depths 
Vo-So = 650 MeV and Vo+So = —66 MeV, the diffuseness parameter a = 0.6 fm, 
and R = 7 fm [241]. Resonance parameters are obtained with the real stabilization 
method [227]. The results are shown as open diamonds and squares for /41/2 and 
jiz, respectively. It is found that splittings of energy and width both become 
smaller compared with the results with the square-well potentials. The reason is 


that the derivative of X(r) is smaller due to a non-zero diffuseness parameter. 
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Figure 26: (Color online) Energy (the left panel) and width (the right panel) splitting between the 
d pseudospin doublets as a function of the potential depth C. Taken from Ref. [249]. 


By examining the Jost function, one can trace continuously the PSS partners 
from the PSS limit to the case with a finite potential depth. This has been done 
in Ref. [249] in which were found a threshold effect in the energy splitting and 
an anomaly in the width splitting of pseudospin partners when C varies from zero 
to a finite value. As the depth of the single-particle potential becomes larger, the 
PSS is broken more and a threshold effect in the energy splitting appears which 
can be seen in the left panel of Fig. 26: The energy splitting first increases then 
decreases until the pseudospin doublets encounter the threshold where one of the 
levels becomes a bound state and the splitting takes a minimum value; When the 
potential becomes even deeper, the splitting increases again. When the depth of 
the single-particle potential increases from zero, there also appears an anomaly 
in the width splitting of pseudospin partners which is shown in the right panel of 
Fig. 26: It first decreases from zero to a maximum value with negative sign, then 
increases and becomes zero again; after the inversion of the width splitting, the 


splitting increases and reaches a positive maximum value, then it becomes smaller 
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and eventually reaches zero. 

In summary, the pseudospin and spin symmetries in single-particle resonant 
states in the Woods-Saxon potentials and self-consistent relativistic mean-field 
potentials have been investigated by using the analytical continuation in coupling 
constant method and the complex scaling method. Not only the pseudospin-orbit 
splittings in energy and width, but also the wave functions of pseudospin partners 
have been systematically examined. A rigorous verification of the pseudospin 
symmetry in single-particle resonant states was given by studying the Jost function 
of the lower components of the Dirac wave functions. In an extensive study of the 
pseudospin symmetry in single-particle resonant states in square well potentials, 
a threshold effect in the energy splitting and an anomaly in the width splitting of 
pseudospin partners were found when the depth of the single-particle potential 
varies from zero to a finite value. There are still many open problems concerning 
the pseudospin and spin symmetries in single-particle resonant states, as listed in 


Refs. [248, 249]. 


3.6. From nucleon spectra to anti-nucleon spectra 


In the spectrum of a Dirac Hamiltonian, one finds single-particle states with 
positive energies in the Fermi sea as well as those with negative energies in the 
Dirac sea. The latter are interpreted as anti-particles under a charge conjugation. 
In the relativistic mean-field theory, the no-sea approximation is made: Those 
negative-energy states are assumed to be empty [285]. Therefore, less attention 
was paid to the negative-energy states. However, if the eigenstates of the Dirac 
Hamiltonian are used as basis states in a Hilbert space, the negative-energy states 
in the Dirac sea must be included for the completeness. When Zhou, Meng, and 


Ring developed the relativistic mean-field theory in a Dirac Woods-Saxon basis 
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[187], they included and examined in detail the negative-energy states and found 
that the pseudospin symmetry of those negative-energy states in the Dirac sea, or 
equivalently, the spin symmetry in the anti-nucleon spectra is very well conserved 
[195]. In Ref. [195], it was shown that the spin symmetry in the anti-nucleon 
spectra is much better developed than the pseudospin symmetry in normal nuclear 
single-particle spectra. The spin symmetry in the anti-nucleon spectra of a nucleus 
was later tested by investigating the relations between the Dirac wave functions of 
the spin doublets and examining these relations within the relativistic mean-field 
theory [197]. Since then, the study of the spin symmetry has often been made in 
combination with that of the pseudospin symmetry in various local potentials, see, 
e.g., Refs. [74, 92]. The spin symmetry in the Dirac negative-energy spectrum 
and its origin were also investigated within the relativistic Hartree-Fock theory 
in which the potentials are non-local [198]. Note that it was found that the 
equality of the vector and scalar potentials in the Dirac Hamiltonian results in 
the spin symmetry in Refs. [66, 336], where the authors suggested applications 
to meson spectra. This symmetry was revealed to be valid for mesons with one 
heavy quark [337]. It should also be noted that, by applying the charge conjugate 
transformation, the spin symmetry for the anti-nucleon states have been formally 
conjectured in Ref. [196]. 

In this Section, we will introduce the spin symmetry in the anti-nucleon 
spectra in nuclei, the test of the spin symmetry by examining the wave functions, 


and the spin symmetry in single-(anti-)particle spectra in various potentials. 


3.6.1. SS in single-anti-nucleon spectra 
In this Section, a nucleon state is explicitly labelled with “N” and an anti- 


nucleon state with “A” for convenience. The Dirac equation for nucleons reads 
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[æ - p + Vur) + B(M + Sn(r))] Wn, s) = ExWn(r), (122) 
where Vy(r) = V(r) and Sx(r) = S(r). For a spherical system, the Dirac spinor 
Wn has the form (cf. Eq. (9)) 


1 
Wnt) = = 


r 


(123) 


NI = 


TROLAR) 
$ ; j=l 
FY) 

Charge conjugation leaves the scalar potential Sn(r) invariant, while it 
changes the sign of the vector potential Vy(r). That is, for anti-nucleons (labelled 
by “A’), Va(r) = —Vy(r) = —V(r) and S a(r) = Sn(r) = S(r). Charge conjugation 


of Eq. (123) gives the Dirac spinor for an anti-nucleon, 


1 
Wa(r) = - 


r 


(124) 


blaine) ied 
TROLAR O } = 
with K = —K. 

For particles there are positive- and negative-energy solutions; the same is 
true for anti-particles. For positive-energy states of the Dirac equations, the 
normal quantum numbers follow the upper components which are dominant. A 
particle state is labelled by {n/m}, while its pseudo-quantum numbers are {ñm}. 
Following Ref. [272], ñ = n + 1 for x > 0; ñ =n for x < 0 [cf. Eq. (84)]. An anti- 
particle state is labelled by {ñm} and its pseudo-quantum numbers are {nlm}. In 


analogy to Ref. [272], the following relation holds for anti-nucleon states, 
n=n+1 for k>0; n=n for k<0O. (125) 


With «(K + 1) = R(K — 1) = I(l + 1) and «(K — 1) = (č + 1) = U1 + 1) in mind, one 


derives the Schrédinger-like equations for the upper and the lower components 
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Table 4: Relation between symmetries and external fields. Taken from Ref. [195]. 


Particle Anti-particle 
dA/dr =0 Spin symmetry Pseudospin symmetry 
dx/dr=0  Pseudospin symmetry Spin symmetry 


[cf. Eqs. (12) and (14)], 


1 æ 1 dM, d 1 Wl+1 1 dM, x +eyG(r), 
aca ae gran [Mt go + ae Sesh ow = 
i ad - = gee —eG(r), 
(126) 
and 
1 æ 1 dM_d 1 Wi+1 1 dM_k +F (r), 
(Reet ae ALAE pe 4 Jo- ee 
M_dr Me dr dr M_ r M- dr r —e,F(r), 
(127) 


where M,(r) = M — A(r) + e and M_(r) = -M — X(r) + e with € = +e for particle 
states or —e€, for anti-particle states. Both equations are fully equivalent to the 
exact Dirac equation with the full spectrum of particle and anti-particle states. 
The relation between SS or PSS and the external fields is given in Table 4. If 
dA/dr = 0, there is an exact SS in the particle spectrum and exact PSS in the anti- 
particle spectrum, because states with the same / but different x are degenerate 
in Eq. (126), where / is the orbital angular momentum of particle states and 
pseudo-orbital angular momentum of anti-particle states. When dA/dr # 0, the 
symmetries are broken, but if dA/dr is so small that the SO term is much smaller 


than the centrifugal barrier, there will be approximate symmetries. 
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Figure 27: Anti-neutron potential M — V(r) + S(r) and spectrum of !ĉO. For each pair of spin 
doublets, the left level is with č < 0 and the right one with k > 0. The inset gives neutron potential 


M + V(r) + S(r) and spectrum. Taken from Ref. [195] and modified to present notations. 


Similarly, when dX/dr = 0 in Eq. (127), there is an exact PSS in the particle 
spectra. On the other hand, for anti-particle states, the SS is exactly conserved 
because now / is the orbital angular momentum. If d£/dr + 0 but small, there are 
approximate PSS in particle spectra and approximate SS in anti-particle spectra. 
This implies that the SS in the anti-particle spectra has the same origin as the 
PSS in particle spectra. However, as revealed in Ref. [195], there is an essential 
difference in the degree to which the symmetry is broken: The factor 1/M? = 
1/(e—X— M) is much smaller for anti-nucleon states than that for nucleon states. 
The bound anti-particle energies €, are in the region between M—A(O) S ea SM, 
approximately 0.3 GeV < ea < 1 GeV. On the other hand the bound particle states 
are in the region of M — |X(0)| < en S M, i.e., for realistic nuclei close to 1 GeV. 
Then |M_(e,)| > 2|M — S (0)| and |M_(ey)| < |X(O)|. Thus the factor in front of 


the K term for anti-particle states is smaller than for particle states by more than 
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a factor (2|M — §(0)|/|Z(0)|)? ~ 400. The SS for anti-particle states is therefore 
much less broken than the PSS for particle states [195], as shown in Fig. 27 by 


taken the neutron and anti-neutron spectra in '!°O as examples. 


3.6.2. SS in single-anti-nucleon wave functions 
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Figure 28: (Color online) Radial wave functions (the left panel) and the differential relation 
(129) of the lower components (the right panel) for the anti-neutron spin doublets 1pj1/2 (€ = 
271.91 MeV) and 1p3/2 (€ = 271.55 MeV) in 40Ca. Taken from Ref. [197] and modified to 


present notations. 


Since the spin-orbit term in Eq. (127) is so small, for a pair of spin partners 
in the anti-nucleon spectrum, the dominant components of radial wave functions 


should be very similar to each other in nuclei and be the same at the SS limit, i.e., 
F,(r)=Fj(r) if dł/dr=0, (128) 
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with j< (j>) labelling the j = Ï- 1/2 (j = Ï+ 1/2) orbital. This relation has been 
tested in Refs. [195, 197]: The dominant components F(r) are nearly identical 
for the two spin partners; on the other hand, their small components G(r) show 
dramatic deviations from each other. However, the small components should 


satisfy a certain relation at the SS limit [197], 


d i+ d 2 
(5 ' =) C= (5 7 “G09. (129) 


Here the radial quantum number ñ is omitted for brevity. 

In Ref. [197], it was examined to what extent the relations given in Eqs. (128) 
and (129) are fulfilled in nuclei. The RMF calculations were performed for *°Ca, 
Zr, !4Sn, and SPb with the effective interaction NL3 [289]. A good SS is 
found in both the anti-proton and anti-neutron spectra. The radial wave functions 
for the anti-neutron 1p doublets in *°Ca are shown in the left panel of Fig. 28. 
The energies of these two states are 271.91 and 271.55 MeV, respectively. One 
can see that the upper components F(r) of the eigenfunctions for the spin doublets 
are almost identical with each other due to the good SS. But the lower component 
G(r) of the wave function of an anti-neutron state deviates dramatically from that 
of its spin partner. In the right panel of Fig. 28, the differential relation of the 
lower components given in Eq. (129) is presented. This differential relation is 
satisfied remarkably well, which gives a further support to the SS in the anti- 


nucleon spectra in nuclei. 


3.6.3. SS in single-anti-nucleon spectra in local potentials 
In the past decade, there have been intensive investigations of the SS in single- 
particle or anti-particle spectra in various local potentials. For some special 


potentials, exact or approximate analytical solutions can be obtained, which 
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makes it very convenient to discuss not only the PSS but also the SS. The readers 
are referred to Section 2.2 for further details. 

The HO potential is widely used and discussed in nuclear physics. For the 
study of the PSS and/or the SS, the RHO potentials are of the most relevance. 
Several examples have been illustrated in Sections 2.2, 3.2, and 3.4. 

In Ref. [148], the eigenfunctions and eigenenergies for a Dirac Hamiltonian 
with equal scalar and vector HO potentials for spherical, axially deformed, and 
triaxially deformed shapes are derived. It has been shown that under the condition 
of equal scalar and vector potentials, the spectrum has a SS. In particular, for 
the spherical case, a higher symmetry analogous to the SU(3) symmetry of the 
non-relativistic HO was discussed. 

In Ref. [92], the generalized RHO Hamiltonian in 1+1 dimensions was solved. 
Both positive and negative quadratic potentials were considered and the bound- 
state solutions for particles and anti-particles were discussed. The main features 
of these bound states are the same as the ones of the generalized three-dimensional 
RHO bound states. The solutions found for zero pseudoscalar potential are related 
to the SS and PSS of the Dirac equation in 3 + 1 dimensions. It has been shown 
how the charge conjugation and y° chiral transformations relate the several spectra 
obtained and that for massless particles the SS- and PSS-related problems have the 
same spectrum but different spinor solutions. 

In Ref. [74], the SS and PSS in the spectra of nucleons and anti-nucleons 
in SPb are studied in scalar and vector Woods-Saxon potentials with different 
depths. Lisboa et al. examined the SO and PSO couplings for selected spin 
and pseudospin partners in both spectra. To assess the perturbative nature of the 


spin and pseudospin symmetries for particles and anti-particles, the SO and PSO 
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Figure 29: Left panel: Spin-orbit terms splittings for three neutrons and anti-neutrons spin partners 
plotted against the respective energy splittings AE. The thicker solid line represents the values for 
which AES® = AE. Right panel: Pseudospin-orbit terms splittings for three neutrons and anti- 
neutrons spin partners plotted against the respective energy splittings AE. The points labelled by 


‘X?’ mark the point of vector potential Vo = O for each pair of levels. Taken from Ref. [74]. 
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contributions are defined for the energy of the level with quantum numbers nk 


(particles), 


ji A’ l+k = 
ESO == Gl? ar f Grud” ° d > 
"= Í (Enx + 2M — Ay r | | r í 0 | | r r 


i X l-kK ee 
gon- f Fpl? 2d f Fpl? 2d ‘ 
7 i Ger iNee IMA) r Mim art y Ena rar 


(130) 


and with quantum numbers ñk (anti-particles), 


i 5 1+x ba 
ES? = ne Fr Ad f Fazl? ad š 
ik f eM Pal rar ; Fag! 1 dr 


co AN’ 1 E lee} 
EPSO — f EG rdr l { Gael? rdr. 
0 (Gx + A)(Gx + 2M + x) r 0 


(131) 


In the left panel of Fig. 29 are shown the SO splittings for three spin partners 
for both neutrons and anti-neutrons, i.e., the difference of the SO term defined 
in Eq. (130) for neutron spin partners and those in Eq. (131) for anti-neutron 
spin partners. In the right panel of Fig. 29 are shown the PSO splittings for three 
pseudospin partners, again for both neutrons and anti-neutrons. From these figures 
one can clearly see the different behaviors of spin and pseudospin in nuclei for 
both neutrons and anti-neutrons. There is a correlation between the values the SO 
splittings and the energy splittings for spin partners, the ratio AES°/AE being very 
close to 1 for anti-neutrons. This is a sign of the perturbative nature of SS in nuclei 
for both neutrons and anti-neutrons. The situation for the pseudospin partners is 
completely different. There is no correlation between the PSO splittings AEPS° 
and the energy splittings AZ, even for small values of AE. Thus it was concluded 


that the PSS in nuclei is not perturbative for both neutrons and anti-neutrons 
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[74]. The readers are referred to Section 4.1 for more discussions of the non- 


perturbative behaviors of PSS. 


3.6.4. SS in single-anti-nucleon spectra in non-local potentials 

In the above investigations, the SS in the (anti-)particle spectrum of local 
potentials has been discussed. In fact, the PSS and its origin as well as 
the importance of the Fock terms have also been investigated [45, 171, 175]. 
Although the PSS was still found to be a good approximation in the RHF theory, 
its mechanism becomes rather complicated by the presence of the non-local 
potentials, see Section 3.3 for details. 

In Ref. [198], the SS in the Dirac negative-energy spectrum and its origin were 
investigated in non-local potentials within the RHF theory. Taking the nucleus 
16O as an example, the SS in the negative-energy spectrum was found to be a 
good approximation and the dominant components of the Dirac wave functions 
for the spin doublets are nearly identical. In comparison with the relativistic 
Hartree approximation where the origin of SS lies in the equality of the scalar and 
vector potentials, in RHF the cancellation between the Hartree and Fock terms 
is responsible for the better SS properties and determines the subtle spin-orbit 
splitting. 

In the RHF theory, the radial Dirac equations (91) are coupled integro- 
differential ones due to the non-local Fock terms [182]. By introducing the 
effective local potentials Xç, Xr, Yg, and Yr defined in Eqs. (92), the integro- 
differential equations can be formally rewritten as equivalent differential ones 
as shown in Eq. (93). Finally, one is able to estimate the Hartree and Fock 
contributions to the SO splittings in the Dirac negative-energy spectrum by using 


Eqs. (96). 
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The denominator X(r) — E in Eqs. (96) contains a state-dependent potential 
Y¢(r). However, the quantity Y¢(r) is around a few MeV and is negligible in 
comparison with X(r) — E which is of the order of 1 GeV. Similar argument also 
holds for the time component of the vector potential V(r) which contains the 
rearrangement term from Fock channels [198]. 

Within the RMF framework, it has been pointed out that the strong centrifugal 
barrier and weak spin-orbit potential lead to the PSS in the single-nucleon 
spectrum [68] and the SS in the single-anti-nucleon spectrum [195]. In Ref. [198], 


this point has also been examined in the RHF scheme. 


160 neutron . 160 neutron 


Figure 30: (Color online) Centrifugal barriers Vcg and spin-orbit potentials Vso multiplied by 
the factor #F?/(= — E) for the spin doublets vlp (the left panel) and v4p (the right panel) in 
the negative-energy spectrum of '°O. The insets show the Hartree contributions of the spin-orbit 


potentials. Taken from Ref. [198] and modified to present notations. 


In Fig. 30 are shown the centrifugal barriers Vcg and the spin-orbit potentials 
Vso multiplied by the factor FF?/@Œ — E) for the neutron spin doublets 1p and 
4p, and their integrals over r are respectively proportional to their contributions to 
the single-particle energy. It is clearly shown that the contribution of Vcg is much 


larger than that of Vso. Therefore, similar reasons as in the RMF theory lead to 
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the SS in the negative-energy spectrum in the RHF theory, and the SO splitting is 
due to the different SO potentials Vso of the spin doublets. 

In the insets of Fig. 30 are given the Hartree contributions to the SO potentials. 
It was found that the contributions from the Fock terms to Vso are one order of 
magnitude larger than those from the Hartree terms. Therefore, the Fock terms 


must play important roles in the SO splittings. 
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Figure 31: (Color online) Hartree and Fock contributions to spin-orbit splitting in the negative- 
energy spectrum of '°O versus the average energy of the spin doublets. The vertical dashed line 


shows the continuum limit. Taken from Ref. [198] and modified to present notations. 


From Eq. (96), the contributions to the single-particle energies E from 
different channels can be estimated quantitatively. To examine the role of the Fock 
terms, the contributions from the Hartree and Fock channels to the SO splittings 
in the negative-energy spectrum of '°O versus the average energies of the spin 
doublets are shown in Fig. 31. It is found that the absolute contributions from both 
Hartree and Fock parts decrease monotonously with the average energy Ew. The 


contributions from the Hartree terms have an energy dependence similar to those 
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in the RMF scheme [195]. It is also found that the contributions of Vip are one 
order of magnitude larger than Vi, as shown in Fig. 30, but they are substantially 
cancelled by the other Fock contributions, i.e., the VE and V, terms [198]. In total, 
the contributions from the Fock terms have an opposite tendency and cancel with 
the Hartree ones, thus leading to better spin symmetry. The competition between 
the Hartree and Fock terms determines the sign of the spin-orbit splitting. 

In summary, the spin symmetry in the anti-nucleon spectra has been shown 
to be very well conserved in nuclei and has also been tested by examining 
the wave functions. The spin symmetry in single-(anti-)particle spectra in 
various potentials, including local and non-local potentials, has been investigated. 
Nowadays, the studies of the spin symmetry and the pseudospin symmetry are 
usually combined together and many interesting topics arise, e.g., the perturbative 
nature of these two symmetries and the contribution of Fock terms. For the spin 
symmetry in anti-nucleon spectra, the problem concerning the polarization effects 


of an anti-nucleon on this symmetry is still open. 


3.7. From nucleon spectra to hyperon spectra 


As it has been emphasized in Ref. [195], the annihilation probability of the 
anti-nucleon in the nucleus is very large and makes it very difficult to observe 
the small spin-orbit splitting of the anti-nucleon levels experimentally. Due to 
the additional strangeness degree of freedom, it is expected that the annihilation 
probability of an anti-hyperon in a normal nucleus is much smaller than that of an 
anti-nucleon. Therefore, it might be easier to observe the small spin-orbit splitting 
of the anti-hyperon levels, when the spin symmetry is approximately conserved. 
In Refs. [200-202], Song, Yao, and Meng investigated the spin symmetry in the 


single-A spectrum, as well as the corresponding polarization and tensor effects 
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with the relativistic mean-field theory. 
In this Section, the spin symmetry in the single-A spectrum will be introduced, 


and the polarization and tensor effects will be discussed. 


3.7.1. SS in single-anti-Lambda spectra 
In the RMF theory, the A hyperon is described as a Dirac spinor moving in the 


potentials generated by the meson fields [338, 339], 
{a - p + Val) + PIMA + Sa) VA) = QUA), (132) 


where M; is the rest mass of A and chosen as Mx = 1115.7 MeV, e; is the single- 
particle energy. As A is charge neutral and isoscalar, it couples only to the o and 
w mesons. As a consequence, the scalar Sx(r) and vector Vq(r) potentials are 


given by, 
SA) = 8A"), VAC) = 8.Av(r). (133) 


According to the charge conjugation transformation, the coupling constants for A 


are related to those for A by the following relations, 


Eci = Zoa > Ewi = E8 wA . (134) 


The Dirac equation for A (132) can be solved similarly as those for nucleons. 
In Ref. [200], the SS in the A spectrum in atomic nuclei was studied by taking 
160 as an example. With the mean-field and no-sea approximations, the coupled 
Dirac equations for nucleons and A together with the Klein-Gordon equations 
for mesons can be self-consistently solved. The effective interaction PK1 [290] 


is adopted for the nucleon part, and the coupling constants for A are chosen as 
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S+V (MeV) 


Figure 32: (Color online) Potential and spectrum of A in '°O. For each pair of spin doublets, the 
left level is with x < 0 and the right one with x > 0. The inset gives the potential and spectrum of 


A in !6O. Taken from Ref. [200]. 


é = 1, while g,4 = goa = 2/38cn, and 2,4 = —8on = —2/38un according to the 
SU(3) symmetry in naive quark model. 

The potential and single-A spectrum in '°O are plotted in Fig. 32, where for 
each pair of spin doublets, the left level is with x < 0 and the right one with x > 0. 
For comparison, the potential and single-A spectrum in !°O are given as well. As 
seen in Fig. 32, the single-A energies for each spin doublets are almost identical, 
and the energy differences between spin doublets eq j. — €,,;, in the A spectrum 
are around 0.09 ~ 0.17 MeV for the p states, which are much smaller than that in 
the A spectrum, 2.26 MeV. 

In order to see the splitting and its energy dependence more clearly, the 


reduced A SO splittings, 


AEso = (Ez j — Ex,j.)/2l + 1), (135) 


104 


0.20 


0.16 

= 0.12 
® 
= 

o 0.08 
m’ 
<j 

0.04 

0.00 

-600 -400 -200 0 
E (MeV) 


Figure 33: (Color online) Reduced spin-orbit splittings AEso for A and anti-neutrons in !ĉO as a 
function of the average energy E,,. For each / orbital, from left to right read the radial quantum 


numbers n = 1,2,... Taken from Ref. [200]. 
in '°O are plotted in Fig. 33 as a function of the average energies, 
Ew = (EX j. + Ez j, )/2 > (136) 


where Ex = €x — Mx excluding the rest mass of A. For comparison, the reduced 
SO splittings for anti-neutrons are also plotted. It was found that the AE so for the 
p States in the A spectrum are around 0.03 ~ 0.06 MeV, which is much smaller 
than those both in the A spectrum, 0.75 MeV, and in the anti-neutron spectrum, 
0.09 ~ 0.18 MeV. This indicates that the SS in the A spectrum is even better 
conserved than that in the anti-neutron spectrum shown in Ref. [195]. The main 
reason is that the SO coupling term for A is about 2/3 of that for anti-neutrons 
[200]. 

In Fig. 34, the radial wave functions F(r) and G(r) for several A spin doublets 


in '°O are plotted. Since the SO splittings in the single-A spectrum are so 
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Figure 34: Radial wave functions of A spin doublets in '°O. Taken from Ref. [200] and modified 


to present notations. 


small, the dominant components F(r) of the wave functions of spin doublets are 
almost identical, while the small components G(r) are quite different. The nodal 


relation (125) is satisfied. 


3.7.2. Polarization effects on SS 

In the above discussions, the polarization effects of A on the SS was neglected. 
For a A really inside !°O, i.e., the vo system, the mean fields including the scalar 
and vector ones are modified by the A [340]. 

In Ref. [201], the polarization effect on the SS due to the A in the '°O+A 
system has been studied with the reduction factor £ = 0.3 for A-meson couplings. 
The polarization effect from the valence A hyperon leads to a highly compression 
of nucleus with the central density up to 2 ~ 3 times of the normal saturation 
density, as seen in Fig. 35. As a result, the energy differences between the spin 


doublets in the A spectrum are around 0.10 ~ 0.73 MeV for the p; states, which 
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Figure 35: (Color online) Density distributions for A (solid line) and neutrons (dashed-dotted 
line) in the '°O+A system. For comparison, the density distribution for neutrons in !°O is plotted 


with the dashed line. Taken from Ref. [201]. 


is larger than the results without polarization effect, 0.09 ~ 0.17 MeV, but still 
much smaller than the SO splittings in the A spectrum, 2.26 MeV. The dominant 
components of the Dirac spinor for the A spin doublets are found to be near 


identical. 


3.7.3. Tensor effects on SS 

The tensor force has been discussed over decades. Recently, the tensor force 
was shown to have a distinct effect on the evolution of the nuclear shell structure 
and appropriate conservation of pseudospin symmetry, see Section 3.4. The 
importance of tensor effects on reducing the spin-orbit splitting of the A single- 
particle spectrum has been extensively discussed in the single-A hypernuclei 
[341-344]. Therefore, it is essential to examine further the SS of A in A-nucleus 


system with the presence of AAw-tensor coupling, 
a, 
Tx = -— iy: VV}, (137) 
2m A 
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with a = f.,x/g a in the RMF theory. This has been done in Ref. [202]. Note that 
the tensor coupling discussed here is also of the Lorentz type, as shown with T(r) 


in Eq. (102). 
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Figure 36: (Color online) Single-particle energies for A with (a = f,,x/g,, = —1) and without 
(æ = 0) tensor coupling in '*C+A. For comparison, the insets show the corresponding results for 


Ain m C. Taken from Ref. [202]. 


Taking the !?C+A system as the first example, the effects of tensor coupling 
on the SS in the single-A spectrum were studied in Ref. [202] by using the 
self-consistent RMF theory with the effective interaction PK1-Y1 [345, 346]. 
Figure 36 shows the single-particle spectrum for A in !?C+A. In order to illustrate 
the tensor effects on the SO splittings, the single-particle spectrum for A without 
(æ = 0) tensor coupling is also plotted. For comparison, the corresponding results 


for A in °C are given in the insets. It is shown clearly that the SO splittings of 
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each spin doublets for A are much smaller than those for A if the tensor coupling is 
not considered. However, the opposite is true after taking into account the tensor 
coupling (w = —1), i.e., the SO-splitting size becomes negligible for A states as 
found in Refs. [342, 343, 347-349], but quite noticeable for the A states. 

In Fig. 37 are shown the scalar potential S(r), the vector potential V(r), 
and S(r) — V(r), as well as the SO potentials for A in '?C+A calculated by 
the self-consistent RMF theory without and with the AAw-tensor coupling. For 
comparison, the corresponding results for A in (°C are given as well. As seen in 
Fig. 37, the vector potential of A changes its sign because of G-parity symmetry. 
The derivative of the difference between the vector and scalar potentials changes 
dramatically with the radial coordinate only for r < 1.5 fm, which leads to the 
central part of SO potential decreasing rapidly to zero at ~ 1.5 fm. However, for 
Ain aG, it is shown that the difference between the vector and scalar potentials is 
quite large. As the consequence, the corresponding central part of SO potential is 
much larger than that for A. 

Of particular interest is the onset of almost opposite phenomena after taking 
into account the tensor effects. The contribution from tensor coupling (“Tensor”) 
reduces the SO potential for A, but enhances that for A. These effects can be 
observed from the splittings of SO partner states, as seen in Fig. 36. 


In Table 5, the SO splittings AEso = E;. — E;, are given for A in °C and for 


j< 
A in *C+A calculated by the RMF theory without (œ = 0) and with (œ = —1) the 
tensor coupling. To show the tensor effects, in Ref. [202] the expectation values 
of the SO potentials with the dominant components in the Dirac spinors were 
divided into two parts, “Central” and “Tensor”, respectively. The difference of the 


expectation values of total SO potentials between the spin doublets AEso gives 
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Figure 37: (Color online) Comparison of the scalar, vector, and total potentials (upper panels) and 
spin-orbit potentials (lower panels) for A in ree and A in C+A calculated by the self-consistent 


RMF theory without (œ = 0, dashed line) and with (œ = —1, solid line) the tensor coupling. Taken 
from Ref. [202]. 
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Table 5: Spin-orbit splittings of A in {C and of A in '*C+A calculated by the self-consistent 


RMF theory without (a = 0) and with (a = —1) tensor coupling. In the calculations with 


tensor coupling, the expectation values of the SO potentials labelled with “Central”, “Tensor”, 


and “Total” in Fig. 37 are calculated with the dominant components in the Dirac spinors of spin 


doublets. Their differences are shown respectively in column AEso. All energies are in units of 


MeV. The data are taken from Ref. [202]. 
AEso 
Ape 
Central Tensor 
Bo ip 151 147 =1.20 
lp 0.64 0.64 1.85 
2p 0.33 0.32 1.03 
ld 0.48 0.50 2.87 


lf 028 0.30 3.18 


PCA 
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Total 
0.27 
2.49 
1.35 
3.37 
3.48 


A Eo! 


0.26 
2.49 
1.37 
3.37 
3.47 


mainly the observed SO splittings. 

As seen in Table 5, the SO splitting for the p; states of °C is 0.26 MeV, which 
is in agreement with the corresponding data 152 + 54 + 36 keV [350]. For A, the 
SO splittings of the 1p, 2p, ld, and 1f states with the tensor coupling are found 
to be 1.37 ~ 3.47 MeV, which is an order of magnitude larger than those without 
the tensor coupling, 0.28 ~ 0.64 MeV. 

It is noted that the SO splitting without the tensor coupling is almost the 
same as the “Central” part of the SO splitting in the calculations with the tensor 
coupling. This indicates that the tensor coupling has negligible contribution to 
the “Central” part of SO potential through the rearrangement of the mean fields. 
However, the additional contribution from the tensor coupling to the SO potential 
of A, corresponding to the “Tensor” term, dominates the final SO splittings in the 
calculations with the tensor coupling. Table 5 shows clearly that the “Tensor” part 
of the SO splitting almost cancels the “Central” part for the A states in )°C, but 
enhances that for the A states greatly in CC+A. 

In Fig. 38 are shown the radial wave functions for the pa states in °C and 
px States in '7C+A calculated by the RMF theory with and without the tensor 
coupling. It shows clearly that the tensor effect is significant on the dominant 
components of Dirac spinors, which recovers the SS on the wave functions of the 
Px spin doublets. For A, the SS is also well conserved from the calculations with 
and without the tensor coupling, because the SO potential of A (~ 1 MeV) is much 
smaller than the corresponding total potential Vā + Sã (~ 280 MeV). Therefore, 
the changes of SO potentials due to the tensor coupling has negligible influence 
on the final wave functions of the A states. 


The tensor effects on the SO splittings for A have been studied systematically 


112 


F,G(r) 


r (fm) 


Figure 38: (Color online) Radial wave functions for the pa states in rc (left panels) and px 
states in '7C+A (right panels). In each case, the top panel represents the results without tensor 
coupling (œ = 0) and the lower part displays the results with tensor coupling (@ = —1). Taken 


from Ref. [202] and modified to present notations. 
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Table 6: Spin-orbit splittings for different single-particle states of A in '©O+A, 4°Ca+A, and 
208Pb+A calculated by the RMF theory without (a = 0) and with (a = —1) tensor coupling. 
All energies are in units of MeV. The data are taken from Ref. [202]. 
ages —— 0o agent 
Central Tensor Total 

lp 0.39 0.40 148 1.88 1.88 

2p 0.23 0.22 0.89 1.11 1.12 

IO+A ld 0.29 0.30 2.11 241 241 

2d 0.16 0.16 0.95 1.11 1.12 

If 0.18 0.19 2.30 2.49 2.48 

lp 0.26 0.28 1.22 1.50 1.48 

2p 0.23 0.23 0.90 1.13 1.14 

“Ca+A ld 0.08 0.09 0.73 0.82 0.82 

2d 0.17 0.18 1.29 1.47 1.46 

1f 0.05 0.05 0.70 0.75 0.75 

lp 0.12 0.15 0.64 0.79 0.76 

2p 0.15 0.15 0.60 0.75 0.76 

8PH+A Id 0.00 0.00 0.05 0.05 0.05 

2d 0.03 0.03 0.27 0.30 0.30 

1f 0.00 0.01 0.06 0.07 0.06 
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for A-nucleus in different mass regions, including '°O+A, *°Ca+A, and 7°°Pb+A 
[202]. The corresponding results are given in Table 6. The tensor effects on the SO 
splittings for these three cases are similar as those for ?C+A. Specifically, the SO 
splittings of A in the calculations with the tensor coupling are found to be 1.12 ~ 
2.48 MeV in '©O+A, 0.75 ~ 1.48 MeV in *°Ca+A, and 0.05 ~ 0.76 MeV in 
208Pbh+A, which are an order of magnitude larger than those from the calculations 
without the tensor coupling, i.e., 0.16 ~ 0.39 MeV in '!°O+A, 0.05 ~ 0.26 MeV 
in *°Ca+A, and 0 ~ 0.15 MeV in 78Pb+A. Moreover, it was found that the SO 
splittings for A decrease with the mass number A, no matter the tensor coupling 
is considered or not. 

In summary, it has been shown that the spin symmetry in the A spectrum 
is even better conserved than that in the anti-nucleon spectrum. Even if the 
polarization and tensor effects, both enlarging the spin-orbit splittings, are 


included, the spin symmetry in the A spectrum remains well conserved. 


3.8. From spherical nuclei to deformed nuclei 


After it was proposed in spherical nuclei, the pseudospin symmetry in 
deformed nuclei was found to be an important physical concept also in axially 
deformed nuclei [13—16] and even in triaxially deformed nuclei [17, 18]. 

The pseudospin symmetry in deformed nuclei has been investigated exten- 
sively. On one hand, a lot of nuclear structure phenomena have been interpreted 
in connection with the pseudospin symmetry, such as nuclear superdeformed 
configurations [29, 30], identical bands [31-33], and pseudospin partner bands 
[35, 36]. On the other hand, in early years, much efforts were devoted to 
revealing connections between the normal spin-orbit representation and the 


“pseudo” spin-orbit one and to exploring the microscopic origin of the pseudospin 
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symmetry with deformed harmonic oscillator potentials [17, 18, 30, 55-58]. In 
particular, it was found that the pseudospin symmetry is conserved exactly for the 
Nilsson Hamiltonian with one-body orbit-orbit (vy) and spin-orbit (1,) interaction 
strengths satisfying the condition v; = 4v,. Moreover, this condition is found to 
be consistent with the relativistic mean-field results [30]. 

Following the success of understanding the pseudospin symmetry in spherical 
nuclei [65, 68], the study of the pseudospin symmetry within the relativistic 
framework was quickly extended to deformed systems [77, 78, 81, 82, 178, 179]. 
First, the pseudospin symmetry limit in deformed nuclei was discussed [78] 
and the evolution of the pseudospin symmetry with deformation was revealed 
[77]. Later, the pseudospin symmetry in realistic deformed nuclei was tested by 
examining the single-particle wave functions [82]. In Ref. [180], the influences 
from different fields of mesons on the pseudospin symmetry were investigated for 
deformed nuclei. Quite recently, the similarity renormalization group was used to 
study the pseudospin symmetry in axially deformed Dirac Hamiltonian [264]. 

In this Section, we will focus on progress on the pseudospin symmetry in 


deformed systems made within the relativistic framework. 


3.8.1. PSS in deformed nuclei 
In Ref. [78], two kinds of conditions were derived for the PSS in deformed 


nuclei. One is exact, 


ax ox 


—=0, —=0, 138 
Or, Oz me) 


and the other is approximate, 


1 0M_Q 


©? 
M_ Or, ry Pp 


2 
L 


(139) 


bd 


r 


where rı = yx +y?, M- = -M -È +e, and Q = Q + 1/2 corresponding to /.. 
The exact conditions correspond to the PSS limit but, similar to the spherical case, 
they are certainly not satisfied in realistic nuclei, while the approximate ones can 
be checked in realistic nuclei. In Refs. [81, 178], the axially deformed RMF theory 
[351, 352] was used to make the numerical check of the approximate conditions 
for the PSS in a well-deformed nucleus '**Sm. It was found that these conditions 
are well satisfied. 

The contributions from different fields of mesons and photons to the PSS 
has been investigated for spherical nuclei [177]. In Ref. [180], such study was 
extended to deformed nuclei and it was found that (i) the o- and w-fields are 
dominant in influencing the PSS, and (ii) the PSO splitting is mainly determined 


by the cancellation of these two fields. 


3.8.2. Evolution of PSS with deformation 

Since the PSS is observed in both spherical and deformed nuclei, how the 
PSS evolves with the shape of nuclei becomes an interesting problem. In 
Ref. [77], the broken PSS was investigated in axially deformed nuclei and the 
evolution of the PSS with deformation was studied by making the constrained 
RMF calculations [351, 352] with the effective interaction NL3 [289]. The quasi- 
degenerate pseudospin doublets were confirmed to exist near the Fermi surface 
for deformed nuclei. 

In Fig. 39, the single-particle states corresponding to pseudospin doublets in 
IDy are plotted against the deformation 62. The asymptotic Nilsson quantum 
numbers [N, n3, A, Q] are good for large values of the deformation parameter 62. 
The pseudospin doublets [Ñ, ñ, A,Q = A + 1/2] [55] are indicated by [N, ñ, A] 


with fT and |. For zero deformation (G2 = 0) the orbitals are indicated by the 
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Figure 39: Single-particle spectra in the deformed nucleus '**Dy as a function of the quadrupole 


deformation parameter 82. Asymptotic pseudospin quantum numbers are given and the pseudospin 
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partners are indicated by arrows T and |. Taken from Ref. [77]. 
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corresponding spherical states. In Fig. 39, it is found that the PSO splitting 
stays almost constant and does not vary with deformation for p2 sufficiently large. 
Furthermore, the energy difference between the | and 7 partners always remains 
positive except for [404], where there appears a crossing at 8 ~ 0.3. Similar 


crossing has been observed in Ref. [55]. 


3.8.3. PSS in deformed single-particle wave functions 

One consequence of the fact that the PSS is a relativistic symmetry of the 
Dirac Hamiltonian is that the relativistic wave functions of the corresponding 
pseudospin doublets satisfy certain relations. The relations (71) and (72) have 
been tested in spherical nuclei [77, 79, 80]. In Ref. [82], the Dirac wave functions 
in deformed nuclei were examined extensively. The conditions of the Dirac single- 
particle wave functions constrained by the PSS were derived and tested in the 
self-consistent RMF calculations. 

In summary, the pseudospin symmetry is connected to many experimental 
observations in deformed nuclei. Theoretically, much effort has been devoted 
to the study of the conditions of the pseudospin symmetry, the evolution of 
the pseudospin symmetry with deformation, and the pseudospin symmetry in 
single-particle wave functions. In a recent work [264], the pseudospin and spin 
symmetries of the Dirac Hamiltonian with axially deformed scalar and vector 
potentials were studied by using the similarity renormalization group theory, see 
Section 4.3.3 for details. There are still many open problems, e.g., the evolution 


of the pseudospin symmetry with the triaxial deformation. 
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4. Open Issues on PSS and SS 


4.1. Perturbative or not 


Since the pseudospin symmetry was recognized as a relativistic symmetry 
of the Dirac Hamiltonian [65], the perturbative nature of this symmetry has 
become a hot topic. The main concern is there are no bound states at the exact 
pseudospin symmetry limit, dX(r)/dr = 0 (17), thus the pseudospin symmetry is 
always broken in realistic nuclei. The non-perturbative behaviors of pseudospin 
symmetry have been carefully considered since the quantitative investigations 
by Marcos et al. in 2001 [73]. It is found that the small splitting of the two 
pseudospin partners cannot be justified by the smallness of the pseudospin-orbit 
term (18), but rather by a complicated counter balance between contributions 
from different terms due to the characteristics of the whole Dirac equation 
itself. Following Arima’s definition of a dynamical symmetry [70], Alberto et 
al. [71, 72] related such non-perturbative behaviors to the dynamical nature of the 
pseudospin symmetry. 

In a sense, such kind of complicated counter balance between contributions 
from different terms indicates the pseudospin-orbit term should not be consid- 
ered as the appropriate symmetry-breaking term [84]. This triggered recent 
investigations on the alternative PSS limits, from which a smooth transition to 
realistic nuclei can be performed [84, 255]. In particular, Typel [255] derived the 
pseudospin symmetry-breaking potential by using the supersymmetric quantum 
mechanics, and showed that the relativistic harmonic oscillator potential is the 
simplest case leading to a vanishing pseudospin symmetry-breaking potential. 
Therefore, the Dirac Hamiltonian with the relativistic harmonic oscillator poten- 


tial can be regarded as an alternative pseudospin symmetry limit, based on which 
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one hopes to understand the pseudospin symmetry breaking in realistic nuclei in 
a perturbative way. 

In Ref. [250], Liang et al. used the perturbation theory to investigate the 
spin and pseudospin symmetries of the Dirac Hamiltonian and their breaking 
in realistic nuclei. The perturbation corrections to the single-particle energies 
and wave functions were calculated order by order. In such a way, the link 
between the single-particle states in realistic nuclei and their counterparts in the 
symmetry limits is constructed explicitly. It is found that the energy splitting of 
the pseudospin doublets can be regarded as a result of small perturbation around 
the Dirac Hamiltonian with the relativistic harmonic oscillator potentials. 

In this Section, we will focus on these discussions on the non-perturbative or 


perturbative nature of the pseudospin symmetry. 


4.1.1. Non-perturbative behaviors of PSS 
In Ref. [73], the Schrédinger-like equation (14) for the lower component of 


the Dirac spinor was expressed as 


aps, +2M(S + V) + 2EV +(S?- V’) - E’ |F = 2MEF, 


mE s) D 


M \F r r2 


(140) 


where E = e— M excluding the rest mass of nucleon and M_(r) = E- V(r)-8S (r). 
It was proven in Ref. [167] that the term (M’/M_)(F’/F — K/r) is not singular 


even at the so-called singularity point rọ with M-(ro) = O, and thus, the integrals, 


e.g., 


[oe] M 
f T F’dr, (141) 
0 = 


121 


are finite. Similar to Eq. (99), one can evaluate the contributions of the terms 


containing F”, F’, S, V and/or E, Ī, and x to the single-particle energy E, 


E(F”) + E(F’) + E(S, V, E) + E(Ď + E(k) = E. (142) 


Table 7: Contributions from different terms in Eq. (142) to the single-particle energy E in Ca 
calculated by the RMF theory with the effective interaction NL3 [289]. All units are in MeV. The 
data are taken from Ref. [73]. 
mo State E(F”) E(F’) E(S,V,E) EQ Elk) E 
508.2 2sı2 18.53 —0.86 —45.16 12.34 -1.82 -16.96 
k=0 13.14 3.07 —49.86 11.66 0.00 -21.98 
ldz 14.00 1.99 —47.58 11.76 3.65 —16.17 
541.0 2sız 19.01 -0.82 —46.39 12.79 -1.54 -16.95 
k=0 13.68 3.29 —51.33 11.97 0.00 -22.40 
ldz2 14.93 2.14 —49.14 12.10 3.02 -16.95 
565.0 2sız2 19.26 -0.76 —47.10 13.12 -1.28 -16.77 
k=0 1404 3.47 —52.22 12.18 0.00 -22.52 
ld32 15.65 2.24 -50.15 12.31 2.38 —17.56 


In Table 7, the contributions from different terms in Eq. (142) to the single- 
particle energy E are shown by taking the 2s1;2 and 1d3;2 pseudospin doublets in 
“°Ca as examples. The results obtained by using three different values of the o- 
meson mass are compared. In order to illustrate the self-consistent effects caused 
by the x term, the results appearing in the intermediate line correspond to the 
solution of Eq. (142) by taking x = 0, i.e., neglecting the PSS-breaking term but 
keeping the PCB. 
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It is pointed out from these results that [73]: (1) For a pair of pseudospin 
doublets, the difference between the contributions of the x term |E(k,) — E(k;)| is 
not negligible in comparison with the contributions of their PCB, [E(,), E(i,)]. 
(ii) In all cases considered, the difference between the contributions to the 
single-particle energy coming from the different terms [[E(F7) — E(F;)|, IEP) — 
E(F;)|,...], in particular of the x term, is larger than the net energy splitting 
|E, — Epl itself, except for the PCB term. (iii) The x term cannot be considered at 
all as a perturbation, as it can be inferred from the two previous assertions. 

Point (ii) reveals the importance of the self-consistent effects due to the x term, 
which modify the small component F (r) of the Dirac spinor. These effects are not 
so small as it has been supposed. As a result, they yield the differences between 
the contributions to the single-particle energy from most of terms in Eq. (142) 
larger than the net splitting |E, — Epl itself. 

Point (iii) becomes more evident if one compares the results for the 2.1/2 state 
(k = —1) to the corresponding fictitious case of x = 0. Supposing the correction 
is perturbative, the effect of the x term on the energy of the 2s,/. state should be 
negative. However, the total energy E of the 2s; state indeed increases a few 
MeV when such a x term is included. That is to say, E,--; — E,=9 has an opposite 
sign than the case that the x term can be regarded as a small perturbation. 

Furthermore, the ordering of the pseudospin doublets strongly depends on the 
particular model considered, mainly through the value of m,, which has a strong 
influence on the smoothness of the nuclear surface. An exact degeneration of the 
pseudospin doublets can even be achieved with a particular value of m,. Thus, the 
accomplishment of the PSS does not necessarily require a small x term. 


Therefore, based on the Schrédinger-like equation for the lower component 
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of the Dirac spinor and the energy decomposition given in Eq. (142), the small 
splitting of the two pseudospin partners cannot be justified by the smallness of the 
k term, but rather by a complicated counter balance between contributions from 


different terms due to the characteristics of the whole Dirac equation itself [73]. 


4.1.2. PSS as a dynamical symmetry 

The non-perturbative behaviors of PSS was then related to the dynamical 
nature of PSS [71, 72], in the sense of Arima’s definition of a dynamical symmetry 
[70]: (i) a symmetry of the Hamiltonian which is not geometrical in nature, or (ii) 
an ordered breaking symmetry from dynamical reasons. 


By using the Woods-Saxon potentials in the Dirac equation (10), 


Xo Ao 
X) = Tame “i AW) = ae 


(143) 
the sensitiveness of the PSO splitting with the potential parameters Xo, Ao, R, 
and a has been investigated systematically [71]. It is found that the degeneracy 
of pseudospin doublets is very much dependent on the shape of the nuclear 
mean field potential. The energy splittings decrease as a increases and R or |Xo| 
decreases. Varying |Zo| and R simultaneously but keeping their product |Xo|R? 
fixed, the pseudospin doublet splittings remain almost constant. That is the reason 
why the PSS is better realized for neutrons than for protons in an isotopic chain 
towards neutron-rich nuclei. It is also found that the pseudospin levels can cross 
each other if a increases or R, |Xo| decreases enough. This stresses an aspect of the 
dynamical nature of PSS [72]. 

As a Step further, in Ref. [72] the contributions to the single-particle energies 


E are decomposed into 


2 
(=) + (Veso) + (Vparwin) + (©) = E, (144) 
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corresponding to the kinetic, PSO, Darwin, and mean-field potential terms defined 
in Eq. (109). This decomposition is quite close to that shown in Eq. (142), while 
( p’/(2M ) stands for the sum of E(F”) and E(Ď. The calculated results can be 
found in Table 3 for the cases without tensor potential. Taking the 1% pseudospin 
doublets as an example, it is clearly seen that the contribution to the PSO splittings 
of the PSO term is larger than the net splitting itself, and it is at least of the same 
order of magnitude but has the opposite sign compared to the kinetic and (£) 
terms. 

Therefore, it is concluded that the onset of the PSS in nuclei is dynamical, 
since it results mainly from cancellations of several terms contributing to the 
single-particle levels, instead of being a consequence of a small PSO coupling 


[71, 72]. 


4.1.3. RHO as a PSS limit 

In the above studies, either the non-perturbative or the dynamical behaviors of 
the PSS shows that the approximate PSS in realistic nuclei cannot be justified by 
the smallness of the PSO potential (18) in the Schrédinger-like equations (14) for 
the lower component F(r). That is to say, this term should not be considered as the 
appropriate symmetry-breaking term [84]. An explicit and quantitative connection 
between the ideal PSS limits and realistic nuclei was still missing. That triggered 
recent investigations on the alternative PSS limits, from which a smooth transition 
to realistic nuclei can be performed [84, 255]. 

In Ref. [84], Marcos et al. aimed at separating the single-particle Hamiltonian 


H as H = Ho + W, where Hp exhibits exact degeneracy of pseudospin doublets, 
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1.€., 


E, = Ey, (145) 


while W represents the PSS breaking term. Various kinds of potentials have 
been investigated, including the harmonic oscillator, Woods-Saxon, relativistic 
Nilsson, and self-consistent RMF potentials. Among these potentials, the RHO 
potential exhibits the exact PSS, in which X(r) is a harmonic oscillator and A(r) is 
a constant. Different from Eq. (71), at such a PSS limit, the lower components of 


a pair of pseudospin doublets are no longer equal to each other, i.e., 


Fr) + Fp(r). (146) 


By comparing their single-particle wave functions, it was pointed out that 
the transition from the RHO model, satisfying exact PSS, to a more realistic 
one as the Woods-Saxon or self-consistent RMF potential, with broken PSS, can 
be considered not far from being perturbative [84]. Meanwhile, because of the 
relation shown in Eq. (146), for a pair of pseudospin doublets, the conditions 
E, ~ E, and F,(r) x F;(r) are no longer dependent on each other. 

In Ref. [255], Typel derived the PSS breaking potential based on the Schrédinger- 
like equation for the upper component of the Dirac spinor, by using the SUSY 
quantum mechanics. The details will be presented in Section 4.3, but one of the 
essential conclusions is that the RHO potential is the simplest case leading to a 
vanishing PSS breaking potential in the SUSY framework. 

Therefore, the Dirac Hamiltonian with the RHO potential can be regarded as 
an alternative PSS limit, based on which one hopes for understanding the PSS 


breaking in realistic nuclei in a perturbative way. 
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4.1.4. Perturbative investigation of PSS 

In Ref. [250], the perturbation theory was used to investigate the symmetries of 
the Dirac Hamiltonian (2) and their breaking in realistic nuclei. The perturbation 
corrections to the single-particle energies and wave functions were calculated 
order-by-order. In such a way, the link between the single-particle states in 
realistic nuclei and their counterparts in the symmetry limits can be constructed 
explicitly. 

First of all, following the idea of Rayleigh-Schrédinger perturbation theory 
[353], the Dirac Hamiltonian H in Eq. (2) or (10) can be split as 


H=H).+W, (147) 


or equivalently 


H) =H-W, (148) 


where Hp leads to the exact SS or PSS, and W is identified as the corresponding 


symmetry-breaking potential. The condition 


Wink 


rae «l for m#k, (149) 


where Wmk = (Wm|WlW;.), determines whether W can be treated as a small 
perturbation and governs the convergence of the perturbation series [353]. 
For the spin and pseudospin symmetry limits shown in Eqs. (16) and (17), 


respectively, the Dirac Hamiltonian with exact symmetries reads 


d kK K 
M+rx -— +- M+X%) -— +- 
HSS = dr r HPSS — dr r (150) 
0 d K 2 0 d kK i 
— +- -M+^o — +- -M+A 
dr dr 
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whose eigenenergies are denoted as Eo, and the corresponding symmetry-breaking 


potentials are 


0 0 L-5, 0 
w5 = Wess = i , (151) 
0 A-A 0 0 


In contrast to adopting the Schrödinger-like equations in the above studies 
[65, 71-73], it is remarkable that all involved operators H, Ho, and W are 
Hermitian, and they do not contain any singularity. This allows us to perform 
the order-by-order perturbation calculations. In addition, only the W term 
corresponds to the symmetry-breaking potential within the present decomposition, 
thus the ambiguity caused by the strong cancellations among the different terms 
in the Schrédinger-like equations can also be avoided. Therefore, this method can 
provide a clear and quantitative way for investigating the perturbative nature of 
SS and PSS. For the case that the nature of the symmetry is perturbative, the link 
between the single-particle states in realistic nuclei and their counterparts in the 
symmetry limits can be constructed quantitatively. For the case that the nature of 
the symmetry is non-perturbative, the divergence of the perturbation series can be 
found explicitly. 

The neutrons in '3?Sn will be taken as examples in the following discussions. 
The corresponding potentials and single-particle calculated by the self-consistent 
RMF theory with the effective interaction PK1 [290] are shown in Figs. 3 and 4, 
respectively. 

For the SS case, the values of |Wmg/(Em — E,)| for the spin doublets k = 1f 
are plotted as functions of the energy differences E„ — Ez in Fig. 40, where 
the unperturbed eigenstates are chosen as those of H§* in Eq. (150), and the 


constant potential are chosen as -M + Ap = —350 MeV. It can be checked that 
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Figure 40: (Color online) Values of |Win /(Em — Ep)| versus the energy differences Em — Ex for the 


spin doublets k = 1f. The unperturbed eigenstates are chosen as those of HSS 


r5 and the single- 


particle states m include the states in the Dirac sea and Fermi sea. Taken from Ref. [250]. 


the convergence of the perturbation calculations is not sensitive to the value of Ao. 
For the completeness of the single-particle basis, the single-particle states m must 
include not only the states in the Fermi sea, but also those in the Dirac sea. It is 
seen that the values of |W,,./(Em — Ep)| decrease as a general tendency when the 
energy differences |E,, — E;,| increase. From the mathematical point of view, this 
property provides natural cut-offs of the single-particle states in the perturbation 
calculations. It is critical to find that the largest value of |W,nx/(Em — Ep)l is ~ 0.1. 
This indicates the criterion in Eq. (149) can be well fulfilled. 

The perturbation corrections to the single-particle energies of the spin doublets 
1f are then examined. In Fig. 41, by choosing the unperturbed eigenstates as 
those of HS, the single-particle energies obtained at the exact SS limit, and 
their counterparts obtained by the first-, second-, and third-order perturbation 
calculations, as well as those obtained by the self-consistent RMF theory, are 


shown from left to right. It can be clearly seen that the SO splitting is well 
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Figure 41: (Color online) Single-particle energies of spin doublets 1f obtained at the exact SS 
limit, and by the first-, second-, and third-order perturbation calculations, as well as those by the 


RMF theory. The unperturbed eigenstates are chosen as those of HS. Taken from Ref. [250]. 


reproduced by the second-order perturbation calculations. The perturbation 
corrections to the single-particle wave functions can be examined in the same 
way, and the same conclusions hold [250]. 

Similar perturbation calculations have been performed for the PSS case in 
Ref. [250]. In the PSS case, since there are no bound states in the pseudospin- 
symmetric Hamiltonian H5S, the perturbation calculations are only performed 
from H to H}°®, i.e., the unperturbed eigenstates are chosen as those of H and the 
perturbation is taken as -WPS in Eq. (151). It is critical to find that the largest 
value of |Wmg/(Em — E,)| is about 0.6 for the PSS case, compared to 0.1 for the 
SS case. That is because, although the potentials satisfy |A — Aol > |È — Xo|, one 


should never forget that different components of the Dirac spinor are involved: 
Wat = (Fal (A-A) Fx), Wink = (Gnl Œ - Zo) |Gx) , (152) 


where for the Fermi states the upper component G(r) ~ O(1), and the lower 
component F(r) ~ O(1/10). 


Therefore, from the perturbative point of view, the bridge can be constructed 
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to connect the Dirac Hamiltonian in realistic nuclei with the symmetry limit of 
constant A, but not constant 2. This indicates that the realistic system can be 
treated as a perturbation of the spin-symmetric Hamiltonian. This also confirms 
in an explicit way that the behavior of PSS is non-perturbative, if dx/dr = 0 in 
Eq. (17) is regarded as the symmetry limit. These conclusions are in agreement 
with those given in Ref. [84] by observing the single-particle wave-function 
behaviors, but now they are demonstrated in a quantitative way. 

However, it is pointed out in the previous Subsection that the energy splitting 
of the pseudospin doublets in realistic nuclei can be alternatively considered as 
the breaking of their degeneracy appearing in the Hamiltonian with the RHO 
potentials [84, 255]. In order to assess this statement in a perturbative way, the 


Dirac Hamiltonian H in Eq. (10) is split as 


H = HREOC + wRHO (153) 
with 
M + Xpo E + a 
msj g e T Tj (154) 
—+- -M+^ 
dr r 
and 
2-2 0 
wrHo (155) 
0 A- ^o 


where Zy0(r) = co + C27” has the form of a harmonic oscillator. Here, HRO leads 
to the energy degeneracy of the whole major shell, and W®"° is identified as the 
corresponding symmetry-breaking potential. The constants -M+Ay = —350 MeV 


and M + co = 865 MeV are chosen in HB"°. As discussed before, the perturbative 
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Figure 42: (Color online) Upper panel: Same as Fig. 40, but for the case of the RHO potentials. 
Lower panel: Same as Fig. 41, but for all single-particle states in the pf major shell. The 


unperturbed eigenstates are chosen as those of Ae. Taken from Ref. [250]. 
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properties are not sensitive to these two constants. Meanwhile, the coefficient c2 
is chosen as 1.00 MeV/fm? to minimize the perturbations to the pf states. 

In the upper panel of Fig. 42, the values of |WRHO Emn- E,)| for the 
pseudospin doublets k = 1d are shown as functions of the energy differences 
En — Ex. It is found that the general patterns shown in this panel are the same as 
those in Fig. 40, and the largest perturbation correction is ~ 0.16. This indicates 
that the criterion in Eq. (149) is fulfilled, even though not as well as in the SS case. 

In the lower panel of Fig. 42, the perturbation corrections to the single- 
particle energies of the states in the pf major shell are shown. As shown in 
this panel, both SO and PSO splittings are well reproduced by the third-order 
perturbation calculations. Thus, the link between the pf states in realistic nuclei 
and their counterparts in the symmetry limit with RHO potential can be explicitly 
established. Furthermore, it is found that the single-particle wave functions of 
H can also be reproduced by the second-order perturbation calculations starting 
from HR". 

Therefore, the quantitative connection between the Dirac Hamiltonian in 
realistic nuclei and that with RHO potentials is constructed by using perturbation 
theory. This indicates that the energy splitting of the pseudospin doublets can be 
regarded as a result of small perturbation around the Dirac Hamiltonian with RHO 
potentials, where the exact degeneracy of the pseudospin doublets appears [250]. 

In summary, whether or not the pseudospin symmetry breaking behaves 
perturbatively depends on whether an appropriate symmetry limit is chosen and 
an appropriate symmetry-breaking term is identified. In Refs. [71-75], the non- 
perturbative behaviors of PSS were shown from various aspects, if one takes the 


condition dX/dr = 0 in Eq. (17) as the symmetry limit. This has also been 
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confirmed with perturbation calculations in Ref. [250]. However, by using the 
supersymmetric quantum mechanics, it is shown that a Dirac Hamiltonian with 
the relativistic harmonic oscillator potential is the simplest example that holds 
the exact pseudospin symmetry [255]. Based on this symmetry limit (154), the 
perturbation corrections to the single-particle energies and wave functions have 
been investigated order by order [250]. In such a way, the link between the single- 
particle states in realistic nuclei and their counterparts in the symmetry limits 
is constructed explicitly. It is found that the energy splitting of the pseudospin 
doublets can be regarded as a result of small perturbation around this symmetry 
limit. Further discussions on perturbation calculations with supersymmetric 


quantum mechanics will be continued in Section 4.3. 


4.2. Intruder states 


For the nuclear single-particle spectra in the normal spin-orbit scheme, the 
s orbitals with / = O correspond to the spin singlets, and all other orbitals with 
l > 0 have their spin partners. In contrast, in the pseudospin-orbit scheme, one 
observes that not all states with / > 0 have their own pseudospin partners. More 
precisely, one finds the § orbitals with I = 0 corresponding to the pseudospin 
singlets, and most of the / > 0 orbitals showing up by pairs as pseudospin doublets, 
but not for the 15 1/2, 13/2, 1ds/2, 1f7/2, 189/2 . . . states. These states are in general 
of j = /+1/2 and have no internal nodes in neither the upper nor lower component 
of the Dirac spinor, ng = np = 0. Beyond the 17/2, these particular states are 
usually called the intruder states, which intrude from the the major shell above to 
the shell below, and form the conventional nuclear magic numbers 28, 50, 82, etc. 
For convenience, we call all these states, i.e., Ls1/2, 1p3/2, as well as 1d5/⁄2, 1 fz/2, 


and 1go9/2, etc., the pseudospin-unpaired states. 
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Such a special difference between the spin and pseudospin schemes has been 
one of the longstanding puzzles in the studies of pseudospin symmetry. The origin 
and the meaning of this special feature are still not fully understood. In this 


Section, we will discuss this puzzle from several different and interesting aspects. 


4.2.1. Bound states in non-confining potentials 
Leviatan and Ginocchio [272] have pointed out that, for the single-particle 
bound states in the Fermi sea, the numbers of internal nodes of the upper and 


lower components of their Dirac spinor, ng and np, obey 


nr =ng for x<0O, Nr =ngt1 for x>O0. (156) 


This was proven analytically for the case that the spherical scalar S(r) and vector 
V(r) potentials are local and S (r), V(r) — 0 as r — ov. In this case, the eigenstates 
in the Fermi sea are bound by the X(r) = V(r) + S(r) potential. The demonstration 
is shown in detail in Section 3.2.1. 

In addition, it is generally believed that, for a pair of pseudospin doublets, the 
lower components F (r) of their wave functions are very similar to each other [65], 
as illustrated in Fig. 8. Even in the alternative PSS limit, the Dirac Hamiltonian 
with RHO potentials discussed in Section 4.1.3, it is also true that the numbers 
of internal nodes of their F(r) components are equal. Namely, for any pair of 


pseudospin doublets with x, < 0 and x, = 1 — Ka > 0, there holds 


Np(Ka) = NF(Kp) « (157) 


Combining Eqs. (156) and (157), one naturally obtains 


Ng(Ka) = Ng(kp) + 1. (158) 
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This explains the reason why the nj, orbitals pair with the (n — 1)j< ones in the 
pseudospin scheme. In addition, the eigenstates 1j, with ng = 0 cannot find their 
pseudospin partners simply because ng cannot be —1. 

In such a way, the special feature concerning the intruder states can be 
understood. However, many interesting aspects about this puzzle are still open, 


which will be discussed in the following parts. 


4.2.2. Bound states in confining potentials 

All states with / > 0 have their own pseudospin partners, if the single-particle 
eigenstates in the Fermi sea are bound by the A(r) = V(r) — S (r) potential instead 
of X(r) [88, 273]. It is the case when lim,_,., A(r) > 2M + X(r) shown in Eq. (85) 
is satisfied [88]. 

On one hand, Alberto, de Castro, and Malheiro [273] have proven that the 
number of internal nodes of the upper and lower components, ng and nr, now 


obeys 
nr =ng-1 for x<0O, nr =ng for x>O0, (159) 


for the local and spherical scalar S(r) and vector V(r) potentials. On the 
other hand, the relation (157) is held as pseudospin doublets have similar lower 
components F(r). As a result, one still holds ng(k,) = ng(kp) + 1 in Eq. (158) for 
all pairs of pseudospin doublets with x, < 0 and x, = 1 — xk, > 0. However, it is 
interesting and important to note that now there exist no eigenstates with quantum 
numbers ng = 0 and j = l + 1/2, because np cannot be —1 in the first expression 
of Eq. (159). 

In such a way, all states with / > 0 have their own pseudospin partners. The 


corresponding typical examples are shown in Fig. 14. However, one can comment 
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that in the present case not all states with / > 0 have their own spin partners, i.e., 
there are some spin-unpaired states which do not have spin partners. This may 
be even more serious because the nuclear shell structure and magic numbers are 


substantially changed. 


4.2.3. Zeros of Jost function: Bound and resonant states 

In Section 3.5, it has been shown that the PSS in the single-particle resonant 
states can be investigated by studying the asymptotic behavior of the nucleon 
Dirac wave functions [247]. In particular, by examining the zeros of the Jost 
function in Eq. (118) corresponding to the small component of the radial Dirac 
wave function, one can investigate the PSS not only for bound states but also for 
resonant states. Here, the problem concerning the intruder states will be discussed. 

In Ref. [249], the Jost functions (118) corresponding to the small component 
of the radial Dirac wave function in square-well potentials (117) with C = 
—66 MeV and D = 650 MeV were investigated in detail. The Jost functions 
are plotted as a function of the binding energy E = e — M for several bound pairs 
of pseudospin partners in Fig. 43. For each pseudo-orbital angular momentum /, 
there exist two x’s, one with a positive value x = /+1 and the other with a negative 
value x = —/, respectively. From Fig. 43, one can clearly find that the number of 
zeros for Jost functions with negative x is always one more than that with positive 
k. For example, for Ī = 1 there are two zeros for the d3 j2 States but three for the 
$1/2 States. This means that there are two bound states for d3/. orbitals and three 
for sı;2 orbitals. Therefore, there is always one bound state with x < 0 which does 
not have a partner. This state is simply an intruder state. 

In Fig. 44 are shown the calculated energies of bound and resonant states for 


the / = 2 pseudospin doublets in the square-well potentials as functions of the 
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Figure 43: (Color online) Jost functions J,(£) (in arbitrary unit) on the real E axis for several pairs 
of pseudospin partners in square-well potentials (117) with C = —66 MeV and D = 650 MeV. 
The results for the j = 7+ 1/2 orbitals are denoted as solid and dashed curves, and their zero 


points representing the bound states are denoted as black and red dots, respectively. Taken from 
Ref. [249]. 
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Figure 44: (Color online) Energies of bound and resonant states for p32 and fs;2 with [T= 2in 
square-well potentials (117) as functions of the potential depth C. The results for the j = 7+ 1/2 
orbitals are denoted as solid and dashed curves, respectively. All the levels are paired off except 
the lowest one. The bottom of the potential E = C and the bound state threshold E = 0 are shown 


as dotted lines. Taken from Ref. [249]. 
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potential depth C. As the potential depth varies from 0 to -70 MeV, the energies 
are always paired off except the lowest one which is an intruder states. 

This is an alternative way to show the origin of the appearance of intruder 
states: The lowest zero of the Jost functions (118) with negative x is always 
isolated while the others are paired off with those of the Jost functions with 


positive x [249]. 


4.2.4. A continuous transformation between SS and PSS 

Recently, Desplanques and Marcos [85] re-examined the puzzle of the intruder 
states in a different way, by setting up two single-particle spectra satisfying 
the SS and PSS, and between these two limits constructing a continuous and 
unitary transformation that commutes with the kinetic-energy operator. This was 
performed in the non-relativistic framework but involving non-local potentials. 


The starting point is the potentials 
Vi = Vis Vo = (0 - PVrX@ - P), (160) 


where V(r) and V2(r) are local and spin independent. The operator (o - p) is 
unitary, but it makes the potential V> non-local. It is clear that [V;, o] = 0, and 
consequently the potential V; satisfies the spin symmetry. Analogously, [V2, (© - 
p)o(o - p)] = 0, since (0 - p)(o- p) = 1 and [V2, 0] = 0, and consequently V> 
satisfies the pseudospin symmetry. Therefore, the Hamiltonian H, can be written 
as 


2, 2 
A, = a +t, (op 2 + va) (o-), (161) 


implying that the spectra for the potentials V. and V2(r) are the same. This 
indicates that the pseudospin-symmetric Hamiltonian A has the same spectra as 


those of the spin-symmetric Hamiltonian Hz = p*/(2M) + V2(r). 
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In order to construct a continuous and unitary transformation from the SS to 


PSS limits, a superposition of the above two potentials is made as [85] 
V = (1 -= XV + xV = (1- Vi (r) + (0 - P)V2(rX0 - P), (162) 


where x is supposed to vary from 0 to 1. 


Woods-Saxon model 
Spin symmetry Pseudospin symmetry 


E [MeV] 


Figure 45: (Color online) Single-particle spectra at the SS and PSS limits and the continuous 
transition in between. The potentials are chosen as the Woods-Saxon type with parameters 


appropriate for a medium-size nucleus like *°Ca. Taken from Ref. [85]. 


In Fig. 45, the single-particle spectra at the SS and PSS limits as well as the 
continuous transition in between are shown. In these results, the potentials are 


chosen as the Woods-Saxon form, 


Vo 


EEFE (163) 


Vi(r) = V(r) = 


with Vo = 40 MeV, R = 5 fm, a = 0.65 fm, and x is varied from —0.2 to 1. 
First of all, it is more or less as expected but still very important that a full 


continuity between the spectra at the SS and PSS limits is obtained. Focusing on 
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the single-particle spectrum around x = —0.2, it agrees with the expectation for a 
medium-size nucleus like *°Ca. In particular, the quasi-degeneracies between the 
standard pseudospin doublets can be seen, such as (24/2, 1d3/2) and (2p3/2, 1 f5/2), 
but these combinations do not correspond to those at the hypothetical PSS limit 
with x = 1 discussed here. Looking at how the spectrum changes when going 
from x = —0.2 to 1, a striking feature is observed. The energies of the standard 
pseudospin doublets, (21/2, 1d3/2) or (2p3/2, 1fs/2), tend to get apart from each 
other, whereas those for states (1 51/2, 1d3/2) and (1p3/2, 1 fs/2) get closer and closer 
and eventually coincide at x = 1. These results suggest a pattern quite different 
from the one generally proposed. In such a pattern, not only all states with / > 0 
have their own spin partners at the SS limit, but also all states with / > 0 have their 


own pseudospin partners at the PSS limit [85]. 
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Figure 46: Wave functions R(r) [f (r) in the original figures] of the 15/2 (the left panel) and 21/2 
(the right panel) states from the SS to the PSS limits given by the Woods-Saxon model. Taken 
from Ref. [85]. 


As a step further, the corresponding single-particle wave functions are studied. 
The wave functions R(r) of the 1s1/2 and 2.1/2 states are shown in Fig. 46 from the 


SS to PSS limit. While it confirms that the transition shown in the wave functions 


142 


is continuous, it shows an interesting phenomenon that an extra node progressively 
appears in their wave functions when x varies from 0 to 1. This is indeed the 
case for all the j = / + 1/2 orbitals, whereas no extra node appears for all the 
j = 1-1/2 orbitals [85]. Unfortunately, the analytical demonstrations about the 
nodal structure shown in Section 3.2 cannot be straightforwardly applied, since 
the non-local potentials are involved here. Therefore, in a sense, it is not clear 
whether these states should be labelled as 15/2, 251/2 or as 21/2, 351/2 at the PSS 
limit, according to the number of internal nodes. In Ref. [85], it is argued that 
these extra nodes are not the usual ones, because they occur outside the range 
of the potential. The contribution to the norm beyond the extra node amounts to 
about 3% in the largest case, compared with the usual contribution which amounts 
up to several 10%. 
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Figure 47: Left panel: Wave functions of the 151/2, 251/2, 1p1/2, 1p3/2, and 1d3/2 states in the 
momentum space R(p) [f(p) in the original figure] calculated at the SS limit x = 0. Right 
panel: The corresponding wave functions in the coordinate space &(r) [f (r) in the original figure] 


calculated by Eq. (164). Taken from Ref. [85]. 


To investigate the relations in the wave functions of pseudospin doublets, one 
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gets benefit from the facts that the wave functions of spin doublets are the same at 
the SS limit and the operator (ø - p) for the continuous and unitary transformation 
is local in the momentum space. Thus, the effect of this operator will not affect 
the angular-momentum-independent part of the wave function denoted as R(p). 
In the left panel of Fig. 47, the wave functions of the 1351/2, 251/2, lpi, 
1p3/2, and ldz; states in momentum space R(p) are compared. These results 
are calculated at the SS limit x = 0. Since V,(r) = V(r) is set, the wave functions 
of the 15 1/2 and 1d3/2 states at the PSS limit x = 1 are nothing but those of the 
Ip states at the SS limit. In addition, to make a fair comparison between these 
functions in the coordinate space, the same spherical Bessel function is used for 


the transformation from the p- to r-space, i.e., [85] 


2 
a0 = fèr f aoho: (164) 


The corresponding results are shown in the right panel of Fig. 47. Note that they 
are different from R(r) shown in Fig. 46. 

From these results, it was pointed out that [85] the wave functions of 1s1/2 
and 1d3;2 do not show much qualitative difference. Furthermore, it is interesting 
that they indeed evolve towards each other when x goes from 0 to 1. In contrast, 
the wave functions of the standard pseudospin doublets 2s1;2 and 1d3;2 evidence a 
striking difference for the whole range of x. 

In short, a new pattern of the SS and PSS limits was set up in Ref. [85]. 
There exists a continuous and unitary transformation between these two symmetry 
limits, while the non-local potentials are involved. In this pattern, all states with 
l > 0 have their own spin partners at the SS limit and all states with / > 0 have 
their own pseudospin partners at the PSS limit. By observing the evolutions in the 


single-particle energies and wave functions, it is tempting to assign the pseudospin 


144 


doublets as, e.g., (151/2, 1d3/2), instead of the standard ones, e.g., (281/2, 1d3/2). It 
represents a departure from the standard PSS. To clarify its validity, analytical 


investigations similar to those shown in Section 3.2 are highly demanded. 


4.2.5. Supersymmetric transformation 

Motivated by the special feature of the intruder states in the pseudospin- 
orbit scheme, the SUSY quantum mechanics was used to investigate the spin 
and pseudospin symmetries [254—256, 270, 271]. It is found that by employing 
both exact and broken SUSY, the phenomenon that all states with / > 0 have 
their own pseudospin partners except for the intruder states can be interpreted 
within a unified scheme. Furthermore, the “striking” difference between the wave 
functions of standard pseudospin doublets mentioned above can be understood, 
as they are indeed almost identical by performing the SUSY transformation 
[255, 270] as shown in Figs. 50 and 53 below. It will be the main task in the 
next Section to illustrate the key ideas of the related investigations. 

In summary, the spin-unpaired (pseudospin-unpaired) states, i.e., the states 
with Z > 0 (/ > 0) which do not have spin (pseudospin) partners, are listed in 


Table 8 for the five different schemes discussed in this Section. 


4.3. SUSY and SRG 


In the above Sections, we have discussed many special features of pseudospin 
symmetry, for example, the singularities in the pseudospin-orbit potential (18), the 
different nodal structures for the pseudospin doublets, the perturbative nature, the 
puzzle of intruder states, and so on. These special features strongly motivated 
recent studies on the pseudospin symmetry [254—256, 270, 271] by using the 


supersymmetric quantum mechanics [252, 253], instead of the Schrödinger- 
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Table 8: Spin-unpaired (S-unpaired) and pseudospin-unpaired (PS-unpaired) states within 


different schemes. The symbol “—” means that there is no S-unpaired or PS-unpaired states. 


S-unpaired states  PS-unpaired states 


Non-confining potentials [272] — n=1&joist 1/2 
Confining potentials [88, 273] n=1&j=!1-1/2 — 
Jost functions [247] — n=1&j=14+1/2 


Non-local unitary transformation [85] — — 


SUSY transformation [255, 270] — n=1& j=14+1/2 


like equations for the lower component of the Dirac spinor. By employing 
both exact and broken supersymmetries, the phenomenon that all states with 
i > 0 have their own pseudospin partners except for the intruder states can be 
interpreted within a unified scheme. In addition, the wave functions of each pair of 
pseudospin doublets become almost identical by performing the supersymmetric 
transformation. 

Typel [255] investigated the pseudospin symmetry by using the supersym- 
metric quantum mechanics, based on the Schrédinger-like equation (12) for 
the upper component of the Dirac spinor. Different from the Schrédinger-like 
equation (14) for the lower component, the equation for the upper component has 
no singularities. As a result, a regular pseudospin symmetry-breaking term can be 
obtained. Unfortunately, the effective Hamiltonian involved in this scheme is not 
Hermitian, which prevents us from being able to perform quantitative perturbation 


calculations. 
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Recent works by Guo and coauthors [262-264] bridged the perturbation cal- 
culations and the supersymmetric descriptions by using the similarity renormal- 
ization group. With the similarity renormalization group, the Dirac Hamiltonian 
is transformed into a diagonal form, expanding in a series of 1/M. The effective 
Hamiltonian in the Schrédinger-like equation thus obtained is Hermitian, which 
makes the perturbation calculations possible. 

By combining the similarity renormalization group, supersymmetric quantum 
mechanics, and perturbation theory, Liang and coauthors [270, 271] pointed out 
that it is now promising to understand the pseudospin symmetry and its breaking 
mechanism in a fully quantitative way. 

Alternatively, Leviatan [254] established the supersymmetric scheme directly 
based on the first-order differential Dirac Hamiltonian by using the intertwining 
relation. As the supersymmetric scheme is established directly on the Dirac 
Hamiltonian, the higher-order terms in the effective Hamiltonian transformed 
by the similarity renormalization group can be avoided. However, it is still an 
open problem how to identify the perturbative pseudospin symmetry limit and the 
corresponding symmetry-breaking term within such a scheme. 

In this Section, we will start with a brief introduction of supersymmetric 
quantum mechanics and its application to the Schrödinger and Schrédinger-like 
equations. We will then introduce the similarity renormalization group. Detailed 
discussions will be mainly focused on the perturbation calculations with similarity 
renormalization group and supersymmetric quantum mechanics. Finally, the 
application of supersymmetric quantum mechanics to the Dirac equations will 


be illustrated. 
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4.3.1. Supersymmetric quantum mechanics 
It has been shown that every second-order differential Hamiltonian can 
be factorized in a product of two Hermitian conjugate first-order differential 


operators, 1.e., [354] 
H,=B'B , (165) 


with B- = [B*]". Its SUSY partner Hamiltonian can thus be constructed by [252, 
253] 


H =B B. (166) 


The Hermitian operators 


0 B 0 -iBt 
Qı = Q = ; (167) 
B 0 iB 0 
are the so-called supercharges with respect to the involution 
1 0 
T=T'= : (168) 
0 -1 


because {Q,,T} = {Q2, T} = 0. The extended SUSY Hamiltonian Hs is the square 


of these Hermitian supercharges, 


H, 0 
0 M 


Hs = 0, =O; = (169) 


The supercharges Q1, Q2 and the extended Hamiltonian H;, together with the 


commutators, 


[Hs, Q1] = [Hs,Q2] = 0, (170) 
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and the anti-commutator, 


{Q;,Q2} =0, (171) 


form the most simple example of a supersymmetric algebra. 
Since the extended Hamiltonian Hs is the square of the supercharges, all 


eigenvalues E;(n) of the eigenvalue equation 
Hs Ys(n) = Es(n)¥s(n) (172) 
are non-negative, and the two-component wave functions read 


y(n) | a 


p(n) 


where y(n) and y2(n) are the eigenfunctions of H, and H3, respectively. For 


Ps (n) = | 


each eigenstate with E;(n) > 0, it is an eigenstate for both H; and H3, and the 


corresponding eigenfunctions satisfy 


B- Bt 
= = — 174 
y(n) T, mo y(n) T; Da (174) 


with the normalization factor 1/ YEs). 

The property of SUSY can be either exact (also called unbroken) or broken 
[252]. The SUSY is exact when the eigenvalue equation (172) has a zero energy 
eigenstate Es(0) = 0. In this case, as a usual convention, the Hamiltonian H; 
has an additional eigenstate at zero energy that does not appear in its partner 


Hamiltonian H2, because 
By (0)=0, W(0) = 0, (175) 


i.e., the trivial eigenfunction of H, identically equals zero. The SUSY is broken 


when the eigenvalue equation (172) does not have any zero energy eigenstate. 
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exact SUSY broken SUSY 


Figure 48: (Color online) Schematic patterns of the exact and broken supersymmetries. Taken 


from Ref. [270]. 


In this case, the partner Hamiltonians H, and H, have identical spectra. The 
schematic patterns of both cases are illustrated in Fig. 48. 
For detailed derivations and various examples of SUSY quantum mechanics, 


the readers are referred to Ref. [252] for a review and Ref. [253] for a textbook. 


4.3.2. SUSY for Schrédinger-like equations 

A second-order differential Hamiltonian (165) is the starting pointing of the 
SUSY quantum mechanics presented above, so this scheme is straightforwardly 
applicable to the Schrödinger or Schr6édinger-like equations. Based on the 
Schrédinger-like equation (12) for the upper component of the Dirac spinor, Typel 
[255] investigated the PSS and derived a regular symmetry-breaking term. 

In Ref. [255], the effects of tensor interactions have also been taken into 
account. The corresponding radial Dirac equations read [cf. Eq. (102)] 


| d K ) 
— + - - T(r)| G(r) = MOF), 
dr r 


-4 p ro) Fr) = MG Ar), (176) 
dr r 


where T(r) is the tensor potential and M.(r) = 6, + M — V(r) + S(r). In the 
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Schr6dinger-like equation for the upper component of the Dirac spinor, 
AG(K) Gn = €nxGne » (177) 


the effective Hamiltonian reads 


If æ 1 M’ 
z d pi -arfar (248-7) 
M, r 


E E n E 
a) dr? r? r M, \dr 


+(M+Z). (178) 


The main task in the SUSY description is to construct the operators B% and 


Bz. The particular form of the Hamiltonian in Eq. (178) suggests the ansatz 


1 l d 
B; = (+l, 179 
B: = |ou i E VM) : M,(r) lo ae A a 


where the superpotentials Q,(r) are the functions of r to be determined. Then, the 


SUSY partner Hamiltonians read 


a M, d 
A(x) = BEB, = T2- Oga" 7 (2+5] , 


1 a Ms d M" 3(M’ ig 
H = B- Bt oe ee Paes ace yey +4 Eoo + 
a) = BB, = T jor Q- Mdr 2M, 4M 


| . (180) 


In order to identify explicitly the «(«+ 1) structure and the tensor T(r) terms shown 


in Eq. (178), the reduced superpotentials g,(r) are introduced as [255] 
K 
GO) = Qr) eT): (181) 


The Hamiltonians H, and H, can be further expressed as 


1 d Kkk+1) , K j 
A(x) = —| - — +4 +2qe- — 4. — 2q, T +T 
i(k) val T7 z de Veh Ged 
M’, d k 
- 27 +T’ + + o+2-7)I, 
a dr r ) 
1 Pi Kk-1) 4 K 5 
Hy(k) = —| - — + G+ 2Qe- +4, -2q T +T 
2(K) ral dr 2 qk Ik", dk q 


ee preni 
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In general, the effective Hamiltonian Hg in the Schrédinger-like equation 


(178) differs from the SUSY Hamiltonian H, in Eq. (182) by a constant, i.e., 


A¢(k) = Ai(k) + e(k), (183) 


where e(x) is the so-called energy shift [252]. Thus, the reduced superpotentials 
q,(r) satisfy the first-order differential equation, 


M! (K) 
M, (K) 


q+ (= OF 4 Ja - q. = —M,(k)N(k) -T' . (184) 


Note that N(x) = e(xk) — M — X(r) depends on the energy shift but M, (K) = Enk + 
M — A(r) depends on the single-particle energy. For regular S(r), V(r), and T(r) 


potentials, the boundary condition for the reduced superpotentials reads 
4x0) = 0. (185) 


At small radius, g,(r) behaves asymptotically as a linear function of r, 


M,(k)N(k) + T’ , 


; 186 
1-2k Leni 


limq,(r) = 
r>0 
and at large radius, g,(7) becomes a constant, 


lim g(r) = ¥(M + €x)(M — e(k)), (187) 


if the S ,V, and T potentials vanish there. 

It is very important to examine the asymptotic behaviors of the full super- 
potentials Q,(r), because they determine the type of SUSY [252]. If there is a 
change of sign in Q,(r) when comparing the limits r — 0 with r — ov, the exact 
SUSY follows, and thus there exists a single non-degenerate state at zero energy. 
In contrast, if the Q,(r) does not change the sign between these limits, the SUSY 


is broken, and thus all eigenstates are doubly degenerate with positive energy. 
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In the present cases, Q,(r) are always positive at r — oo. In contrast, 
Q,(r) for r — 0 is determined by the angular-momentum term x/r. One finds 
lim,_.9 O,(r) = —co for x < 0, and lim,_,9 O,(r) = +00 for x > 0. In other words, 
the SUSY is exact for all cases of x < 0, whereas SUSY is broken for all cases 
of x > 0. This is crucial to understand the puzzle of intruder states discussed in 
Section 4.2. 

The x-dependent energy shifts e(x) can be determined in the following ways: 


(i) For the case of x, < 0, the SUSY is exact, and it requires 
€(Ka) = Elk - (188) 


(ii) For the case of x, > 0, the SUSY is broken, and thus the corresponding energy 
shift can be, in principle, any number which makes the whole set of H, eigenstates 
positive. In practice, the energy shifts are determined by assuming that the PSO 
potentials vanish as r — 0. This vanishing behavior is similar to that of the 
usual surface-peaked SO potentials. Therefore, lim,_,9 qx, (r) = lim,—o0 qx,(7) with 
Ka + Kp = 1 is satisfied for pseudospin doublets. Considering M,(x,) and M,(k;) 
are almost identical as €, ~ &, and neglecting the contribution of T’(0), the energy 


shifts read [255] 
e(kp) = 2 (M + ©)|,-0 — e(Ka) . (189) 


Finally, one can evaluate the corresponding PSS-breaking terms. The Hamil- 
tonians H>(k,) + e(K,) and H>(k,) + e(k,) for the pseudospin doublets are almost 


identical, and their difference is given by the potential [255] 


Ka — Kp 


as 2 d Ka Kp 1 
WSS = [Hy (kz) +e(Ka)]—[Ho (kp) +e(Ks)] = a ie a 


, (190) 
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where the difference between M,(xk,) and M,(x,) is neglected. It is essential that 
this symmetry-breaking potential is a regular function for all r, in contrast to that 
shown in Eq. (18). One of the simplest cases that such a symmetry-breaking 


potential vanishes is nothing but the relativistic harmonic oscillator potentials as 


S(r) = Vin = Hr and T(r)=0. (191) 
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Figure 49: (Color online) Reduced superpotentials g,(r) (the upper panel) and full superpotentials 
Q,(r) (the lower panel) for the sı;2 and d3/2 orbitals. Taken from Ref. [255], and the values of x 


for s1/2 and d3/2 orbitals were shown with | and —2 in the original figures, respectively. 


In Fig. 49, a typical example of the reduced q,(r) and full Q,(r) superpotentials 
is shown for the p pseudospin doublets, in which the tensor potential is neglected 


and scalar and vector potentials are chosen as the Woods-Saxon type, S(r) = 
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So/[1 + exp((r — R)/a)] and V(r) = Vo/[1 + exp((r — R)/a)] with So = —450 MeV, 
Vo = 370 MeV, R = 3.8 fm, and a = 0.65 fm. The single-particle energies 
of 1p pseudospin doublets thus obtain are E>,,,, = —15.604 MeV and Eia = 
—15.424 MeV, and the corresponding energy shifts are e(—1)—M = —56.849 MeV 
and e(2)— M = —102.680 MeV. It is shown that the reduced superpotentials 
q,(r) are almost identical for radii below 5 fm. At large r they become constant 
and approach different values for r — oo as predicted by Eq. (187). The full 
superpotentials Q(r) contain the angular-momentum contribution and, hence, they 
diverge for r — 0. The change of sign for Q_;(r) indicates the exact SUSY with 
a single state at zero energy for its H,. In contrast, no change of sign for Q2(r) 
corresponds to the broken SUSY. 

In Fig. 50, the single-particle wave functions for the 1p pseudospin doublets 
are shown. The upper panel shows the wave functions of the original Hamiltonian 
Hg in Eq. (178). No doubt they have different numbers of nodes and look quite 
different. However, as shown in the lower panel, the eigenfunctions of the SUSY 
partner Hamiltonian in Eq. (182) are remarkably similar to each other. They 
are nothing but the original wave functions with the transformation shown in 
Eq. (174). 

Furthermore, it is proven in the SUSY quantum mechanics that the number of 
nodes in the radial wave functions y⁄2(r) is one less than those in their counterparts 
y(r) when the SUSY is exact, while the number of nodes in y2(r) is the same 
as those in their counterparts Wı(r) when the SUSY is broken [252]. This 
indicates that the nodal relation between the pseudospin doublets [272] discussed 


in Section 3.2.1, 


D 


ñ=n-1 for «x<QO, =n for x>0O, (192) 
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Figure 50: (Color online) Upper panel: Single-particle wave functions of the original Hamiltonian 
for the 2sı;2 and 1d3/2 pseudospin doublets. Lower panel: Those of the supersymmetric partner 


Hamiltonian. Taken from Ref. [255] and modified to present notations. 
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is one of the intrinsic properties of the SUSY quantum mechanics [270]. 

The quasi-degeneracy in the single-particle energies, the similarity in the 
wave functions, and the smallness of the symmetry-breaking term (190), all of 
these facts imply that one should be able to understand the PSS breaking in a 
quantitative and perturbative way. Unfortunately, both the effective Hamiltonian 
Hg (178) and its SUSY partner (182) are not Hermitian, since the upper 
component wave functions, as the solutions of the Schrédinger-like equation, 
are not orthogonal to each other. This prevents us from being able to perform 


quantitative perturbation calculations. 


4.3.3. Similarity renormalization group 

Recent works in Ref. [262—264] bridged the perturbation calculations and the 
SUSY descriptions by using the SRG. 

The idea of SRG [265-269] is to drive the Hamiltonian toward a band- 
diagonal form via the flow equation and unitary transformations that suppress 
off-diagonal matrix elements. In recent years, the SRG is also widely used in 
nuclear effective field theory for decoupling the low-energy physics from the high- 
energy physics [355-357]. The SRG-evolved two-nucleon and three-nucleon 
interactions have been used for the ab initio calculations of nuclear ground states 
[358, 359], low-lying spectra [360], nucleon scatterings [361], and so on. The in- 
medium SRG has been developed for the medium-mass open-shell nuclei [362— 
365], which decouples the physics of valence nucleons from the full Hilbert space, 
enabling exact diagonalizations in the valence space that are impossible in the full 
problem where all nucleons are active. Recent reviews on the related topics can 
be found in, e.g., Refs. [366, 367]. 


For the Dirac Hamiltonian shown in Eq. (2), it can be transformed with the 
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SRG into a diagonal form and expanded in a series of 1/M. It is crucial that the 
effective Hamiltonian in the Schrédinger-like equation thus obtained is Hermitian. 
This makes the perturbation calculations possible. 

Fist of all, according to the commutation and anti-commutation relations with 
respect to the 6 matrix, the Dirac Hamiltonian in Eq. (2) is separated into the 


diagonal £ and off-diagonal o parts, H = € + o, which satisfy 
[e,8] =0 and {0,6} =0. (193) 


In order to obtain the equivalent Schrédinger-like equation for nucleons, the main 
task is to decouple the eigenvalue equations for the upper and lower components 
of the Dirac spinor. A possible way is to make the off-diagonal part of the Dirac 
Hamiltonian vanish with a proper unitary transformation. 

According to the SRG [265-269], the Hamiltonian H is transformed by a 


unitary operator U(/) as 
H(D = U(DHUŻ (D, (194) 


with H(/) = e(/) + o(D, H(0) = H, and a flow parameter l. Then, the so-called 
Hamiltonian flow equation can be obtained by taking the differential of the above 


equation, 
d 
ET (D = O, HO], (195) 


with the anti-Hermitian generator 


_ du) 
= dl 


As pointed out in Ref. [268], one of the proper choices of 7(/) for letting the off- 


n(1) ud, (196) 


diagonal part o(/) = 0 as l — o reads 
m1) = [BM, HD]. (197) 
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Then, the diagonal part of the Dirac Hamiltonian e(/) at the / — oo limit can be 
derived analytically in a series of 1/M [262], 
gœ)  &3(00) E E4(00) Ge 

M M? M? 


_ Li a ae, 
= BM + (BS + V) + Bla py + lla PES +V]. a p] 


1 
= CPT ed ae “py a {a ' P, [læ ' P, (BS + V)| (BS + V)I} 


elco) = Melo) + €;(0C0) + 


- 2l: piss + VP) + (198) 


In such a way, the eigenvalue equations for the upper and lower components 
of the Dirac spinor are decoupled. The equivalent Schrédinger-like equations 
for nucleons with Hermitian effective Hamiltonians can be obtained. The 
corresponding details can be found in Refs. [262, 263] for the spherical case and 
Ref. [264] for the axially deformed case. 

For example, for the spherical case, the effective Hamiltonian for the nucleons 


in the Fermi sea up to the (1/M°)-th order reads [262] 


2 
1 
He=M+E4 EE- (sph s’ 


d k ï r D 
2M 2M? dr} r4M? 8M? 
S 5 d\ KSN (XYP-2'A +452" pi 
a LG) da L A A T aA LL 
t Im (s Pr z) r 2M? 16M? 8M? 
(199) 
with the operator p% = -d°/(dr°) + x(k + 1)/r?. This Hamiltonian can be 


decomposed into five Hermitian components: the non-relativistic term H,,, the 


spin-orbit term H,, the dynamical term H4, the relativistic modification of kinetic 
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energy Hz, and the Darwin term H,, as 


H =M+2+ ËE, 
2M 
k ^ KSA’ 


Pr 
Hasti 
E 8M3’ 
w” (XP - 2A’ 4 45D” 
maar =e e ered at 
3m2 16M? on 


Since all these terms are Hermitian, one can calculate the contribution of each 
term to the single-particle energies, which is very helpful to disclose the origin of 
relativistic symmetries. 

In Ref. [264], Guo et al. have made a new exploration of the PSS in 
deformed nuclei by carrying out such a decomposition for axially deformed Dirac 
Hamiltonian with Woods-Saxon-like vector and scalar potentials. The single- 
particle states for “Dy were obtained and the SO and PSO splittings were 
analyzed in detail. 

In Fig. 51, the variation of the energy splitting between the spin (pseudospin) 
doublets with the quadrupole deformation parameter 82 is shown. It can be seen 
that the SO splitting is mainly from the contribution of the SO term H, when the 
system is prolate deformed. However, on the oblate side, the situation becomes 
more complicated: The non-relativistic, the SO coupling, and the dynamical terms 
all influence very much the SO splitting. The different behaviors between the 
SO energy splitting on the prolate and oblate sides are quite interesting. This 
difference was attributed to the presence of configuration mixing on the oblate 


side. This should be examined in detail in the future. 
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Figure 51: (Color online) Comparison of the contributions of all the terms in Hp to the spin (left 
panels) and pseudospin (right panels) energy splitting and their correlations with the deformation 
parameter 62 for four pairs of spin (pseudospin) doublets. Here “nonrela’”, “dynam”, “spinorb”, 
“relakin”, and “Darwin” denote the non-relativistic part, the dynamical term, the spin-orbit term, 


the relativistic modification of kinetic energy, and the spin (pseudospin) Darwin term, respectively. 


As a guide to the eyes, the total energy splittings are marked as “total”. Taken from Ref. [264]. 
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Compared with the SO splitting, the change of the PSO splitting with the 
deformation is more complicated, as one can see in the right panel of Fig. 51. The 
energy difference between a pair of pseudospin partners is mainly determined by 
three parts, the non-relativistic term H,, the spin-orbit term H., and the dynamical 
term H4. Four pairs of pseudospin partners are shown here and their energies are 
in the range of 0 ~ —25 MeV. In this energy range and in the whole shown range 
of £5, the contribution of H, to the PSO splitting is always positive, which means 
its contribution to the pseudospin anti-aligned state is larger. The opposite is that 
of H., i.e., its contribution to the pseudospin aligned state is larger. Therefore, the 
SO term always plays a role in favor of the PSS. The contribution of the dynamical 
term H; to the PSO splitting varies with the deformation and changes its sign when 


the pseudospin doublets approach the threshold. 


4.3.4. Perturbation with SRG and SUSY 

Gathering all pieces presented above, it is promising to understand the PSS 
and its breaking mechanism in a fully quantitative way by combining the SRG, 
SUSY quantum mechanics, and perturbation theory [270, 271]. 

The Dirac equation can be transformed into a diagonal form in a series of 
1/M by using the SRG as shown in Eq. (199), and its lowest-order approximation 
corresponds to a Schrödinger equation. In Ref. [270], taking the lowest-order 
approximation as an example, the idea for applying the SUSY quantum mechanics 
to trace the origin of the PSS was illustrated and the PSS breaking mechanism was 
explored quantitatively by the perturbation theory. 


By assuming the spherical symmetry, the radial Schrödinger equation is cast 
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in the form of 
A(k)R(r) = ER(r) (201) 


with the single-particle Hamiltonian 


d K(k + 1) 


yorsa eee 
©) = sage + Me 


+ V(r), (202) 


and the single-particle wave functions 


Rar) lq 


r JaMa 


Yalt) = (f). (203) 


Here V(r) is the non-relativistic central potential standing for the sum of the scalar 
and vector potentials, Ł(r), in Eq. (199). 

It is clear that H conserves the explicit SS for the spin doublets a and b with 
Ka+k, = —1, which leads to the same centrifugal barrier k(k+1)/(2Mr7). Similarly, 
in order to investigate the origin of the PSS and its breaking mechanism, it is 
crucial to identify the corresponding term proportional to (7+ 1) = x(k — 1), 
which leads to the same x(k — 1) values for the pseudospin doublets a and b with 
Ka + Kp = 1. The SUSY quantum mechanics is one of promising approaches for 
identifying such x(x — 1) structure. 

Following the similar procedures shown in Section 4.3.2, one can start with a 


couple of Hermitian conjugate first-order operators 


d 1 1 d 
By = jou = 2| Sae B; = —— jow es (204) 
dr] V2M V2M dr 
and the reduced superpotentials 
K 
q«(r) = Q,(r) — a (205) 
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and end up with the SUSY partner Hamiltonians 


1 d Kk+1) 
_ ptp- — | _ 
eal a ai dr? r? 


1 | d Kk-1) 


H. = BBY = — |-— 
lk) = BB 2M | dr r? 


2 2K , 
ge Ue Ge > 


: (206) 


2 2K , 
+ dk + —qk + qk 
r 


It is important to note that these Hamiltonians are Hermitian, but those in 
Eqs. (182) are not [270]. 
The reduced superpotentials qg,(r) satisfy the first-order differential equation 


[270], 


na Ér) + Hat) — g (r)| + e(k) = V(r), (207) 


with the asymptotic behaviors 


lim adr) = ¥-2Me(xk) and lim QA’) = —. (208) 


The energy shifts are determined in the same way as that shown in Eqs. (188) and 


(189), i.e, 
C(Ka) = Ex, and = e(kp) = 2 V| =o — C(Ka), (209) 


for the cases of k < 0 and x, > 0, respectively. 
Before the numerical calculations, it is interesting to seek a possible exact PSS 
limit analytically. The SUSY partner Hamiltonian reads 


d K(k — 1) 


ES, «gE Gig, 210 
Mar © ome * Ye”) oy) 


A(k) = Hy(k) + e(k) = 
with 


Vr) = V(r) + q'(r)/M. (211) 
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By definition the exact PSS limit holds Enk, = Eq—1)x, With Ka < 0 and ky +x, = 1, 
which indicates H>(k,) + e(Ka) = H2(k,) + e(k,). By combining Eqs. (206) and 


(207), as well as the boundary condition g,(0) = 0, one can readily have 
Ik (1) = qo (r) = Mar , (212) 


with a known constant w, = (e(k,)—e(K,))/(Kp—Ka). As the reduced superpotentials 
qr) are linear functions of r, the central potential V(r) in H has the form [cf. 


Eq. (191)] 


27? + V(0). (213) 


M 
Vyo(r) = y PE 


The corresponding PSS limit is nothing but the well known case with harmonic 
oscillator potentials, which leads to the energy degeneracy of the whole major 
shell. 

In the following calculations, the mass of nucleon is taken as M = 939.0 MeV, 
and the central potential V(r) is chosen as a Woods-Saxon form 


Vo 


YO) = Ty aria’ 


(214) 


with the parameters V) = —63.297 MeV, R = 6.278 fm, and a = 0.615 fm, which 
corresponds to the neutron potential provided in Ref. [368] by taking N = 82 
and Z = 50. This potential is illustrated as the solid line in Fig. 55 below. In 
this Section, we use a tilde to denote the operators, potentials, and wave functions 
belonging to H. 

In Fig. 52, the reduced PSO splittings (EF. — E;,)/ (27+ 1) versus their average 
single-particle energies Ea = (E;. + E;,)/2 are plotted, where j< (j>) denotes the 
states with j = Ī— 1/2 (j = Ï+ 1/2). It is seen that the amplitudes of the reduced 


PSO splittings are less than 1 MeV. Moreover, as a general tendency, the splittings 
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Figure 52: (Color online) Reduced PSO splittings (E; — E;,)/ (27+ 1) versus their average single- 
particle energies (E; + E;,)/2. Taken from Ref. [270]. 


become smaller with the increasing single-particle energies, which is in agreement 
with the results shown in Figs. 5 and 10 by the self-consistent relativistic mean- 
field and relativistic continuum Hartree-Bogoliubov theories. It is very interesting 
to investigate the physical mechanism for such energy-dependent behavior. This 
also helps to figure out whether or not the PSS is an accidental symmetry [84]. 

In the left panel of Fig. 53, the single-particle radial wave functions R,,;(r) 
of H are shown by taking the pseudospin doublets 2/ and 1f as examples. It is 
clear that the wave functions of the spin doublets are identical since there is no 
spin-orbit term in H. In contrast, the wave functions of the pseudospin doublets 
are very different from each other. This leads to difficulties in analyzing the origin 
of the PSS and its breaking. 

Prior to the quantitative analysis by using the perturbation theory in Ref. [250], 
the investigation of PSO splittings was usually done by decomposing the contri- 


butions by terms, where each contribution is calculated with the corresponding 


166 


nlj 


R (fm) 


Figure 53: (Color online) Single-particle wave functions R,,;(r) of H (the left panel) and R 


A (the right panel) for the 3.51 /2, 2d3/2, 2ds/2, and 1g7/2 states. Taken from Ref. [270]. 
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(r) of 


Table 9: Contributions from the kinetic term (kin), centrifugal barrier (CB), and central potential 


(cen) to the single-particle energies E and the corresponding PSO splittings AEpso for the 


pseudospin doublets 2 and 1f. All units are in MeV. The data are taken from Ref. [270]. 


State 
3514/2 
2d3/2 
AEpso 
2ds/2 
187/2 
AEpso 


Exin 
28.953 
16.845 
12.109 
16.845 

6.197 
10.648 


Ecg 
0.000 
11.758 
—11.758 
11.758 
20.483 
—8.725 
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Ecen 
—50.545 
—51.746 

1.201 
—51.746 
—54.188 

2.442 


E 
=21.591 
—23.143 

1.552 
—23.143 
—27.508 

4.365 


operator O; by 
E; = i R*(r)O;R(r)dr . (215) 


Within the representation of H shown in Eq. (202), the operators of the kinetic 
term, centrifugal barrier, and central potential read —d?/(2Mdr7), K(k + 1)/(2Mr’), 
and V(r), respectively. In Table 9, the contributions from these terms to the 
single-particle energies E as well as the corresponding PSO splittings AEpso are 
shown for the pseudospin doublets 2 and 1f. It is not surprising that, within 
this representation, the contributions to AEpso come from all channels, while they 
substantially cancel to each other in a sophisticated way. 

As discussed in Section 4.1, the phenomenon of such strong cancellations 
among different terms was usually associated with the dynamical [71, 72] and 
even the non-perturbative [73-75] nature of PSS. However, such connection is 
mystified and sometimes even misleading. Indeed, by using the perturbation 
calculations based on the PSS limit shown in Eq. (213), we demonstrated that the 
nature of PSS discussing here is perturbative [270], the corresponding calculated 
results are quite similar to those shown in Fig. 42. 

Much more important is now the origin of the PSS and its breaking mechanism 
can be studied in an explicit way within the representation of the SUSY partner 
Hamiltonian H [270]. 

To obtain the SUSY partner Hamiltonian A in Eq. (210), one solves the first- 
order differential equation (207) for the reduced superpotentials g,(r) with the 
boundary condition g,(0) = 0. The corresponding q,(r) are shown in Fig. 54 in 
units of MeV/c. The q,(r) are x-dependent since the left-hand side of Eq. (207) 
contains a x-dependent term. However, it should be emphasized that g¢,(r) does 


not depend on the radial quantum number n for a given x. One will discover that 
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Figure 54: (Color online) Reduced superpotentials g,(r) for the p and f orbitals. Taken from 
Ref. [270]. 


such an n-independent property is essential for understanding the general pattern 
of AEpso versus Ea shown in Fig. 52. 

Then the x-dependent central potentials V,(r) in H can be obtained, and their 
asymptotic behaviors read 


2(e(k) — V) 


lim V(r) = V + 
a? (1— 2”) 


and lim V,(r) = 0. (216) 


It is important that these potentials are regular and converge at both r — 0 andr —> 
co, In the upper panels of Fig. 55, these central potentials V,(r) are shown, while 
the Woods-Saxon potential V(r) in H is also shown for comparison. For all x, the 
potentials V,(r) approximately remain a Woods-Saxon shape, and they become 
shallower than the original potential V(r). By comparing two upper panels, it is 
seen that the amplitude of the difference between V,(r) for a pair of pseudospin 
partners increases with the difference of their quantum numbers |k, — Kol. 


After getting the central potentials V,(r), one is ready to calculate the single- 
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Figure 55: (Color online) Upper panels: x-dependent central potentials V,(r) in A as a function of 
r for the (a) p and (b) f orbitals, while the Woods-Saxon potential in H is shown for comparison. 
Lower panels: The corresponding single-particle energies obtained with H and H. Taken from 


Ref. [270]. 
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particle energies and wave functions of the SUSY partner Hamiltonians A(x). In 
the lower panels of Fig. 55, the discrete single-particle energies obtained with 
H are compared with those obtained with H. It is clear that the eigenstates of 
Hamiltonians H and H are identical, except for the lowest eigenstates with x < 0 
in H, which are the so-called intruder states. In other words, the fact that the 
intruder states have no pseudospin partners can be interpreted as a natural result 
of the exact SUSY for x < 0 and broken SUSY for x > 0. By holding the one-to- 
one mapping relation in the two sets of spectra, the origin of PSS, which is deeply 
hidden in H, can be now traced by employing its SUSY partner Hamiltonian 7. 


The single-particle radial wave functions R,7,(r) of H for the 2p and 1f 


ñlj 
pseudospin doublets are shown in the right panel of Fig. 53. The nodal relation 
in Eq. (192) can be seen by comparing the wave functions shown in the left and 
right panels. In fact, not only are the numbers of nodes equal, but also the wave 
functions of pseudospin doublets are almost identical to each other. Therefore, 
within this representation, the quasi-degeneracy of pseudospin doublets is closely 
related to the similarity of their wave functions, and vice versa [270]. 

The same strategy as done in Table 9 is then adopted to investigate the PSO 
splittings, but now within the representation of H shown in Eq. (210) instead. 
The corresponding operators include the kinetic term —d*/(2Mdr’), the PCB 
K(k — 1)/(2Mr’), and the central potential V,(r). The corresponding results for the 
pseudospin doublets 2f and 1 f are listed in Table 10. It can be seen that for each 
pair of pseudospin doublets the energy contributions from the PSS-conserving 
terms, i.e., the kinetic and PCB, are very similar. The PSO splittings AEpso 
are mainly contributed by the difference in the central potentials AF..,, which 


is due to the slight x-dependence of V,(r) as shown in Fig. 55. In other words, the 
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Table 10: Contributions from kinetic term (kin), pseudo-centrifugal barrier (PCB), and central 

potential (cen) to the single-particle energies E and the corresponding PSO splittings AEpso for 

the pseudospin doublets 2 and 1f. All units are in MeV. The data are taken from Ref. [270]. 
State Ekin Epcg Ecen E 


2Đı;2 16.602 6.723 -44.916 -21.591 
2Ď32 17.331 6.857 -47.332 -23.143 
AEpso —0.729 -0.134 2.415 1.552 
lf 5.710 16.286 -—45.139 -23.143 
lpn 6.293 16.591 -50.392 -27.508 
AEpso —0.584 -0.305 5.253 4.365 


sophisticated cancellations among different terms in H can be clearly understood 
by using a proper decomposition with the help of the SUSY quantum mechanics 
[270]. 


In order to perform the perturbation calculations, the Hamiltonian A is split as 
H = AS. + WS, (217) 


where HSS and WPS are the corresponding PSS-conserving and PSS-breaking 
terms, respectively. By requiring that W?SS should be proportional to x, which is 


similar to the case of the SO term in the conventional scheme, one has 


: 1 œ Kk-1)| ~- 
an = 7M -5 z + Vpss(r), 
WSS = KVpso(r). (218) 


In such a way, the PSS-conserving Vess(r) and PSS-breaking Veso(r) potentials 
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can be uniquely determined as [270] 


Vess(r) = 


= MV gag Fgo(p) = LEO HO 9 


Ka — Kp M Ka- 


Vass (MeV) 
Toso (MeV) 


Figure 56: (Color online) PSS-conserving potentials Vpss(r) (the left panel) and PSS-breaking 
potentials Vpso(r) (the right panel) for the p and F orbitals. Taken from Ref. [270]. 

In Fig. 56, the Vess(r) and Vpso(r) potentials are illustrated by taking the p 
and f orbitals as examples. On one hand, it can be seen that the PSS-conserving 
potentials Vpss(r) remain an approximate Woods-Saxon shape, and they are K- 
dependent to a small extent. On the other hand, the PSS-breaking potentials 
Vpso(r) show several special features [270]: (i) These PSS-breaking potentials 
are regular functions of r, in particular, they vanish at r — oo. (ii) It can be seen 
that the amplitudes of Vpso are around 1 MeV, which directly lead to the reduced 
PSO splittings AEpso S 1 MeV as shown in Fig. 52. (iii) More importantly, the 
PSO potentials Vpso(r) are negative at small radius but positive at large radius with 
anode at the surface region, which is totally different from the usual SO potentials 
with a surface-peaked shape. The particular shape of the PSO potentials Vpso(r) 
can explain well the variations of the PSO splitting with the single-particle energy. 


More details for point (iii) are as follows: First of all, Veso(r) do not depend 
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on the radial quantum number n. Meanwhile, the single-particle wave functions 
R(r) extend to larger distances with higher energies. Thus, the matrix element 
(R|VpsolR) is negative when the wave function is centralized in the inner part. 
As the wave function becomes more extended, the positive part of Vpso(r) 
compensates for the negative value of the matrix element. In this way, the PSO 
splittings AEpso decrease while the radial quantum numbers ñ increase. The 
splittings can even reverse for the resonance states, where the outer part of the 


PSO potentials plays the major role. 


“i PSS , PSS 
A =A; 


+W 


E (MeV) 


Figure 57: (Color online) Single-particle energies obtained at the exact PSS limit HPSS, and their 
counterparts obtained by the first-, second-, and third-order perturbation calculations with Wess | 


as well as those obtained with H. Taken from Ref. [270]. 


Finally, the perturbation calculations are performed based on the pseudospin 
symmetric Hamiltonian }SS with the perturbation W?SS. In Fig. 57, the single- 
particle energies obtained at the PSS limit H®SS, and their counterparts obtained 
by the first-, second-, and third-order perturbation calculations, as well as 
those obtained with H are shown from left to right. It can be seen that the 


pseudospin doublets are exactly degenerate at the PSS limit YSS, For the present 
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decomposition, the largest perturbation correction |WP8S/(E,, — Ep)| is less than 


0.03 [270]. This indicates that the criterion in Eq. (149) is satisfied very well. As 
shown in Fig. 57, not only the PSO splittings but also the energy degeneracy 
of the spin doublets are excellently reproduced by the first-order perturbation 
calculations. 

In such an explicit and quantitative way, the PSO splittings AEpso can be 
directly understood by the PSS-breaking term WSS within the representation of 
the SUSY partner Hamiltonian #. Furthermore, this symmetry-breaking term can 
be treated as a very small perturbation on the exact PSS limit H?SS. Therefore, the 
PSS discussing here is of pertubative nature [270]. 

Recently, it has been demonstrated in Ref. [271] that the perturbative nature 
of PSS maintains even when a substantial SO potential is included. The SO term 
shows both indirect and direct effects on the PSS-breaking potentials Vpso(r). The 
indirect effect due to the changes of the reduced superpotentials q,(r) is rather 
small. In contrast, the direct effect, i.e., the SO potential itself appearing in both 


H and H, reduces the PSO splittings AEpso systematically and substantially. 


4.3.5. SUSY for Dirac equations 

The above discussions are based on the second-order differential Hamiltonian, 
i.e., the factorizable Hamiltonian. Alternatively, Leviatan established the SUSY 
schemes directly based on the first-order differential Dirac Hamiltonian by using 
the intertwining relation [254]. 

In usual applications of SUSY quantum mechanics, one starts from a factoriz- 
able Hamiltonian H4, then identifies the pair of Hermitian conjugate operators B+ 
and B- in Eq. (165), and eventually obtains the SUSY partner Hamiltonian H in 
Eq. (166). 
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This procedure can be regarded in a different way. Assuming one holds 
BH, = MWB, (220) 


which is the so-called intertwining relation between H; and H3 [369], this 
intertwining relation ensures that, if yı(n) is an eigenstate of Hı, then also 
yaln) x B-w,(n) (174) is an eigenstate of Hy with the same energy Es (n), unless 
B-w,(n) vanishes (175) or produces an unphysical state, e.g., non-normalizable. 
In other words, the SUSY schemes shown in Fig. 48 can be set up as long as the 
intertwining relation is satisfied, but Hamiltonians H, and H; are not necessarily 
factorizable. 

For example, one can insist that both SUSY partner Hamiltonians H(x,) and 
H(k,) be the Dirac Hamiltonian of the form prescribed in Eq. (10), and looks for 


possible solutions of By that satisfy 


waxy 222 Was. <2 

Ba ow a dr r |= A a dr r B7. (221) 
L +2 -M+A() +2 -M+A(Ð) 
dr dr r 


By considering a matrical Darboux transformation operator, 
2 d 
B; = Per) + OA7), (222) 
dr 


where P, and Q, are 2 x 2 matrices, and assuming certain forms in the functions 
[P.(r)]i; and [Q,(r)];;, Leviatan found three different kinds of solutions in 
Ref. [254]. They correspond to three different kinds of symmetry limits: (i) 
Coulomb symmetry limit, (ii) spin symmetry limit, and (iii) pseudospin symmetry 
limit. The schematic patterns in these symmetry limits are illustrated in Fig. 58. 
In the Coulomb symmetry limit, the orbitals with x, + Kp = O form the 


degenerate doublets, e.g., (2512, 1pi/2), (23/2, 1d3/2), and the scalar and vector 
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Figure 58: Schematic patterns in (a) SUSY quantum mechanics and in the (b) Coulomb, (c) 
pseudospin, and (d) spin symmetry limits of the Dirac Hamiltonian. Taken from Ref. [254], and 


the notation of the single-particle spectra starts from n = 0. 
potentials are in the shape of 
S(r) =as/r and V(r) = ay/r, (223) 


respectively. The corresponding transformation operator reads [254] 


d z Ey P Ma, _ as d ay 
-Iid r Ka Ka dr r 

B, = a;sd ay d se Ma l’ (224) 
Ka dr r dr r Ka 


where €, = Ka + Œs @+/Ka and Q, = Qg + Qy. 
The SS limit corresponds to that shown in Eq. (16), i.e., V(r)— S(r) = Ao, and 
the transformation operator reads [254] 


2M + X(r) — Ao —-—+— 


a g k ana C (225) 
—+— 0 
dr r 
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with Ka + Kp = —l. 
The PSS limit corresponds to that shown in Eq. (17), i.e., V(r) + S(r) = Xo, 


and the transformation operator reads [254] 


0 4 + Kb 
Belg g dr r , (226) 
—+— -2M+A(r)- Xo 
dr r 


with kK, + Kp = 1. 

The higher-order terms in the effective Hamiltonian (199) transformed by 
SRG can be avoided, if the SUSY scheme is established directly on the Dirac 
Hamiltonian by using the intertwining relation (220). However, the present studies 
have not included the perturbative PSS limit shown in Eq. (154). Furthermore, it 
is still an open problem about how to identify the corresponding PSS-breaking 
term W?SS in such a scheme. 

In summary, it was pointed out in Refs. [270, 271] that it is promising 
to understand the pseudospin symmetry and its breaking mechanism in a fully 
quantitative way, by combining the similarity renormalization group [262-264], 
the supersymmetric quantum mechanics [254—256], and the perturbation theory 
[250, 251]. It is shown that while the spin-symmetry-conserving term appears 
in the single-particle Hamiltonian H, the pseudospin-symmetry-conserving term 
appears naturally in its supersymmetric partner Hamiltonian H. The eigenstates 
of Hamiltonians H and É are exactly one-to-one identical except for the so-called 
intruder states. In such a way, the origin of pseudospin symmetry deeply hidden 
in H can be traced in its supersymmetric partner Hamiltonian H. 

By using the similarity renormalization group, the Dirac Hamiltonian can 
be transformed into a diagonal form and expanded in a series of 1/M. So 


far, the cases corresponding to the lowest-order approximation and the lowest- 
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order spin-orbit potential have been investigated in Refs. [270, 271], respectively. 
The perturbative nature of pseudospin symmetry has been demonstrated by the 
perturbation calculations for both cases without and with a substantial spin-orbit 
potential. A general tendency that the pseudospin-orbit splittings become smaller 
with increasing single-particle energies can also be interpreted in an explicit way. 
The corresponding studies with the effective mass term and higher-order terms are 
expected in the near future. 

Alternatively, the supersymmetric scheme can be directly established based 
on the first-order differential Dirac Hamiltonian by using the intertwining relation 
[254]. In such a way, the higher-order terms in the effective Hamiltonian 
transformed by the similarity renormalization group can be avoided. Meanwhile, 
it is still an open problem how to identify in such a scheme the perturbative 
pseudospin symmetry limit, e.g., a Dirac Hamiltonian with relativistic harmonic 


oscillator potential, as well as the corresponding symmetry-breaking term. 


5. Summary and Perspectives 


After the independent observation of the near degeneracy between pairs of 
single-particle states (n,/,j = 1+ 1/2) and (n — 1,1 + 2, j = 1+ 3/2), i.e., the 
so-called pseudospin symmetry, by Hecht and Adler [11] and Arima, Harvey, 
and Shimizu [12] in 1969, it raised a fascinating question whether such near 
degeneracy is accidental (a degeneracy not explained by an obvious symmetry) or 
due to symmetry breaking (more descriptively hidden symmetry). The pioneering 
work by examining the ratio between the strengths of the spin-orbit and orbit-orbit 
potentials by Bohr, Hamamoto, and Mottelson in 1982 [55] opened the door to 


explore the hidden symmetry in understanding the origin of pseudospin symmetry. 
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With the success of the relativistic mean-field theory [59-64], it is quite 
encouraging to find that the special ratio between the strengths of the spin-orbit 
and orbit-orbit interactions can be reproduced by the relativistic mean-field theory 
[30]. The reveal of the pseudo-orbital angular momentum as the orbital angular 
momentum of the lower component of the Dirac spinor by Ginocchio in 1997 
[65] proved to be an unexpected success of the relativistic mean-field theory. It 
is then followed by lots of exciting discoveries. For examples, the vanishing of 
the derivative for the sum of the scalar and vector potentials, i.e., dx(r)/dr = 0, 
can lead to the exact pseudospin symmetry [68], which means that the pseudospin 
symmetry becomes much better for exotic nuclei with a highly diffused potential 
[69]. While developing the relativistic mean-field theory in a Dirac Woods-Saxon 
basis [187], Zhou, Meng, and Ring [195] discovered the spin symmetry in the 
anti-nucleon spectra. A rigorous verification of the pseudospin symmetry in the 
single-particle resonant states is given by Lu, Zhao, and Zhou [247] by examining 
the Jost functions corresponding to small components of the radial Dirac wave 
functions. By using the similarity renormalization group, Guo and coauthors 
[262, 263] have made a new exploration of the pseudospin symmetry, including 
the axially deformed systems [264]. Combining the supersymmetric quantum 
mechanics, perturbation theory, and similarity renormalization group, Liang and 
coauthors [270, 271] have provided a promising way to understand the origin of 
pseudospin symmetry and its breaking fully quantitatively. 

All the pioneering works and exciting discoveries have triggered tremendous 
enthusiasms in understanding the pseudospin symmetry in various physical 
systems as well as exploring their origin. Although the present paper is intended 


to provide a comprehensive overview of the related works done recently in the 
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exploration of the pseudospin and spin symmetries, it is unavoidable to have 
overlooked some of them. 

The historical development and recent progress have been summarized in the 
introduction by trying to exhaust the existing literatures. In Section 2, the general 
features for the Dirac equation and its corresponding Schrédinger-like equations 
were discussed. The pseudospin symmetry and spin symmetry in various systems 
and potentials have been reviewed in Section 3. The discussions ranged from 
stable to exotic nuclei, from non-confining to confining potentials, from local to 
non-local potentials, from central to tensor potentials, from bound to resonant 
states, from nucleon to anti-nucleon spectra, from nucleon to hyperon spectra, 
and from spherical to deformed nuclei. Extensive discussions have been devoted 
in Section 4 to the open issues on the pseudospin and spin symmetries as well, 
including the perturbative nature, the puzzle of intruder states, and the connection 
with the supersymmetric quantum mechanics and similarity renormalization 
group. 

The current Review is focused on the theoretical exploration of the hidden 
pseudospin and spin symmetries and their origins in atomic nuclei. There 
might be questions on the unique experimental signal of these symmetries. 
While introducing some themes in the study of very deformed rotating nuclei, 
Mottelson [370] preluded the link between the pseudospin symmetry and newly 
observed experiments. In particular, a lot of phenomena in nuclear structure have 
been successfully interpreted directly or implicitly by the pseudospin symmetry, 
including nuclear superdeformed configurations [29, 30], identical bands [31-33], 
quantized alignment [34], and pseudospin partner bands [35, 36]. In addition, 


the relevance of pseudospin symmetry in the structure of halo nuclei [45] and 
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superheavy nuclei [46, 47] was pointed out. The fingerprint for pseudospin 
doublet bands and their difference with the chiral doublet bands [371] have been 
briefly listed in Ref. [372]. 

The experimental verification of the pseudospin symmetry can also be done 
from the single-particle energies. The observed single-particle energies obtained 
by pickup or knockout reactions contain polarization effects, except for the 
cases where the spectroscopic factors of the single-particle states are close to 
one [303, 373]. Therefore, one should take into account the particle-vibration 
coupling and polarization effects [374, 375]. Litvinova and Afanasjev [374] have 
studied the impact of particle-vibration coupling and polarization effects due to 
deformation and time-odd mean fields on the single-particle spectra systematically 
in doubly magic nuclei from >°Ni up to superheavy ones. It has been shown that 
the particle-vibration coupling substantially improves the description of splitting 
energies in the pseudospin doublets. However, there are still cases where the 
observed pseudospin-orbit splittings are poorly reproduced. 

Finally, the exploration for the deep reasons on the approximate pseudospin 
symmetry observed in realistic nuclei is also a challenging problem. There are 
already some discussions on the physics behind the near equality of the vector 
and scalar potentials in the Dirac Hamiltonian [376-378] and more investigations 
are expected. In particular, it is highly necessary to search for the links with 
more fundamental models on strong interactions or realistic interactions between 


nucleons. 
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Appendix A. 


During the preparation of the present Review, we found the symbols and 
notations used in different papers are indeed quite similar. However, subtle 


differences, like a sign, a factor, etc., are sometimes crucial for understanding the 
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investigations. Therefore, some of the symbols and notations have been changed 
from the original papers and unified through the present Review by choosing 
the common conventions. We also recommend these unified conventions for the 
future studies. 

In this Appendix, we list the key symbols and notations as well as the 


abbreviations used in this Review in an alphabetic way. 


List of key symbols and notations 


A atomic mass number 

Bt, B- pair of Hermitian conjugate operators in SUSY 

E single-particle energy excluding rest mass, E = € — M 
ER single-A energy excluding rest mass, Ex = €} — M 
Ew average single-particle energy, Ea = (Ej. + E;,)/2 

E; contribution from each potential to 


the total single-particle energy E 


F(r) lower component of Dirac spinor in spherical case 
FAP) Fer) = 1 *F(r) 

G(r) upper component of Dirac spinor in spherical case 
GA) GAT) = WGA) 

H,H single-particle Hamiltonian and its SUSY partner 
Ho symmetry-conserving Hamiltonian 

H system Hamiltonian density 

Tk Jost function 

L Lagrangian density 

M mass of nucleon 
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[N, n3, A, Q] 


Qx 
R(r) 


S(r) 
T(r) 
V(r) 
Vece 
Vpso 
W 
wess 
X,Y 


XG, Xr, YG, Vr 


I 
Z 


effective mass, M, = M-A+e 

effective mass, M_ = -M -È + € 

neutron number 

asymptotic Nilsson quantum numbers 

SUSY superpotential 

non-relativistic single-particle radial wave function 
in the spherical case 

scaler potential, negative value means attractive 
tensor potential 

vector potential, positive value means repulsive 
(pseudo-)centrifugal barrier 

(pseudo)spin-orbit potential 
symmetry-breaking term 

PSS-breaking term in SUSY scheme 

non-local Fock potentials 

localized Fock potentials 

spherical harmonics spinor 


proton number 


labels of single-particle orbitals 

nucleon creation and annihilation operators 
SUSY energy shift 

meson tensor coupling strengths 

meson coupling strengths 


run over 1, 2,3 
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total angular momentum of single-particle states 

j = LF 1/2 for spin doublets or 

j = 1¥ 1/2 for pseudospin doublets 

angular momenta of upper component G(r) and 

lower component F (r) 

mass of mesons 

main quantum number of single-particle states, 

starting from n = 1 

number of internal nodes of G(r) and F(r) 

momentum operator 

SUSY reduced superpotential 

radial and angular coordinates in spherical coordinate system 
radial and height coordinates in cylindrical coordinate system 
spin operator s = o` /2 and pseudospin operator 


Š = (0 - p)s(o - P) 


width of a resonant state 

A(r) = V(r) — S (r), positive value means repulsive 

for the case of constant V — S 

reduced spin-orbit splitting, AEso = (Ej. — Ej,)/(21 + 1) 
reduced pseudospin-orbit splitting, AEpso = (Ej. — E;,)/(2/+ 1) 
X(r) = S(r) + V(r), negative value means attrative 
for the case of constant S + V 

X(r) of the form of a harmonic oscillator 


system ground-state wave function 
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l->) 


5 
a, B, y, y, o” 


HV 


Pcs Ps» PvP. 


o 
> 
T, T3 


Wr) 
W(X), WA) 


List of abbreviations 


physical vacuum 


Dirac matrices 

set of quantum numbers, œ = (a, Ma) = (Na, la, Ja, Ma) 
single-particle energy including rest mass M 
single-anti-nucleon and single-A energies including rest mass 
quantum number, x = +(j + 1/2) for j = L+ 1/2 states 
run over 0, 1, 2,3 

charge and scalar densities, isoscalar and isovector 
baryonic densities 

Pauli matrices 

isospin operator and its third component 
single-particle wave function 


single-anti-nucleon and single-A wave functions 


ACCC analytical continuation in coupling constant 
CB centrifugal barrier 

CDFT covariant density functional theory 

CSM complex scaling method 

HO harmonic oscillator 

NU Nikiforov-Uvarov 

PCB pseudo-centrifugal barrier 

PSO pseudospin-orbit 

PSS pseudospin symmetry 
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RCHB relativistic continuum Hartree-Bogoliubov 


RHF relativistic Hartree-Fock 
RHO relativistic harmonic oscillator 
RMF relativistic mean-field 
SO spin-orbit 
SRG similarity renormalization group 
SS spin symmetry 
SUSY supersymmetry or supersymmetric 
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